Tugas Polinomial karakteristik matriks 4x4 dengan program bahasa C++

tugas aljabar linear .

baru belajar pake C++ jadi masi belum tau bagaimana cari akar2 polinomial di C++

kalau ada yang tau mohon di share ya .

ini Program yang aku buat

//Program Aljabar Linear

//Determinan matriks 4x4

#include<iostream>

using namespace std;

int main()

{

    int brs1=4, kol1=4;

    float A[brs1][kol1];

    float CA00, CA11, CA22, CA33, CA01, CA02, CA03, DA11, DA22, DA33;

    float DetA, tr, C, D;

    int i,j;

    char off;

    

    //masukkan A (Matriks kofisien)

    cout << "Masukkan Matriks A :" << '\n';

    cout<< '\n'; 

    for (i=0; i<brs1; i++){

        for (j=0; j<kol1; j++){

            cout<<"baris " << i+1 << " kolom " << j+1 << ": ";

            cin>>A[i][j];

        }

}     

   cout<< '\n';  

    //Menghitung nilai koefisien a b c d e persamaan Polinomial 

       {for (i=0; i<brs1; i++)

        {for (j=0; j<kol1; j++)

               CA00 = (A[1][1]*A[2][2]*A[3][3])+(A[1][2]*A[2][3]*A[3][1])+(A[1][3]*A[2][1]*A[3][2]) -

                       (A[3][1]*A[2][2]*A[1][3])-(A[3][2]*A[2][3]*A[1][1])-(A[3][3]*A[2][1]*A[1][2]);

               CA01 = (A[1][0]*A[2][2]*A[3][3])+(A[1][2]*A[2][3]*A[3][0])+(A[1][3]*A[2][0]*A[3][2]) -

                       (A[3][0]*A[2][2]*A[1][3])-(A[3][2]*A[2][3]*A[1][0])-(A[3][3]*A[2][0]*A[1][2]);  

               CA02 = (A[1][0]*A[2][1]*A[3][3])+(A[1][1]*A[2][3]*A[3][0])+(A[1][3]*A[2][0]*A[3][1]) -

                       (A[2][1]*A[1][3])-(A[3][1]*A[2][3]*A[1][0])-(A[3][3]*A[2][0]*A[1][1]);

               CA03 = (A[1][0]*A[2][1]*A[3][2])+(A[1][1]*A[2][2]*A[3][0])+(A[1][2]*A[2][0]*A[3][1]) -

                       (A[3][0]*A[2][1]*A[1][2])-(A[3][1]*A[2][2]*A[1][0])-(A[3][2]*A[2][0]*A[1][1]); 

               CA11 = (A[0][0]*A[2][2]*A[3][3])+(A[0][2]*A[2][3]*A[3][0])+(A[0][3]*A[2][0]*A[3][2]) -

                       (A[0][3]*A[2][2]*A[3][0])-(A[0][0]*A[2][3]*A[3][2])-(A[0][2]*A[2][0]*A[3][3]);

               CA22 = (A[0][0]*A[1][1]*A[3][3])+(A[0][1]*A[1][3]*A[3][0])+(A[0][3]*A[1][0]*A[3][1]) -

                       (A[0][3]*A[1][1]*A[3][0])-(A[0][0]*A[1][3]*A[3][1])-(A[0][1]*A[1][0]*A[3][3]);

               CA33 = (A[0][0]*A[1][1]*A[2][2])+(A[0][1]*A[1][2]*A[2][0])+(A[0][2]*A[1][0]*A[2][1]) -

                       (A[0][2]*A[1][1]*A[2][0])-(A[0][0]*A[1][2]*A[2][1])-(A[0][1]*A[1][0]*A[2][2]);                

               DA11 = (A[2][2]*A[3][3])-(A[2][3]*A[3][2]);

               DA22 = (A[0][0]*A[3][3])-(A[0][3]*A[3][0]);

               DA33 = (A[0][0]*A[1][1])-(A[0][1]*A[1][0]);  

    

    //Polinomial Karakteristiknya         

         

    tr = A[0][0]+A[1][1]+A[2][2]+A[3][3];    

    C = CA00+CA11+CA22+CA33;

    D = DA11+DA22+DA33;

    DetA = (CA00*A[0][0]*1)+(CA01*A[0][1]*(-1))+(CA02*A[0][2]*1)+

           (CA03*A[0][3]*(-1)); 

           

   cout<< '\n';

   cout<< "Determinan Matriks A = "<<DetA<<'\n';

   cout<< '\n';  

   

    //Hasil

    cout<< "Polinomial Karakteristiknya = t^4 - "<<tr<<"t^3 + "<<C<<"t^2 - "<<D<<"t + "<<DetA<<'\n';

    

    cout<< '\n';

   cout<< '\n';  

   cout<<"== PROGRAM SELESAI =="<<endl<<endl;

   cin>>off;

 }}}