Polinomio característico en c++

[jugo86]

Estimados:

 

Estoy haciendo un programa en c++ que verifica si una matriz de orden 3 es diagonalizable, para lo cual, en primer lugar, desarrollé una función que calcule el determinante de dicha matriz por la regla de Sarrus. Mi inconveniente es crear el polinomio característico a trvés la función que hice. Si sirve de algo, dejo la función sarrus que hice:

 

int sarrus(int m[i][i])
{
    int a, b, maux[i][j], diag1, diag2, diag3,diag4, diag5, diag6, s1, s2, det;
    diag1=diag2=diag3=diag4=diag5=diag6=1;    
    for(a=0;a<i;a++)
    {
        for(b=0;b<i;b++)
            maux[a][b]=m[a][b];
        for(b=i;b<j;b++)
            maux[a][b]= m[a][b-i];
    }    
    
    for(a=0;a<i;a++)
        {
            for (b=0;b<j;b++)
            {
                switch(b-a){
                
                    case 0:
                        diag1*=maux[a][b];
                        break;
                    case 1:
                        diag2*=maux[a][b];
                        break;
                    case 2:
                        diag3*=maux[a][b];
                        break;
                    }
                
                switch(a+b){
                
                    case 2:
                        diag4*=maux[a][b];                        
                        break;
                    case 3:
                        diag5*=maux[a][b];                        
                        break;
                    case 4:
                        diag6*=maux[a][b];
                        break;
                    } 
            }
        }
    s1=diag1+diag2+diag3;
    s2=diag4+diag5+diag6;
    det=s1-s2;
    return det;        
}

 

Espero que me puedan ayudar, (y que la función que hice sirva para algo), desde ya muchas gracias.