Implementar el método de krammer no es dificil, es pura algebra matricial.
Primero genera el código que te calcula el determinante de una mátriz.
Que tipo de sistemas de ecuaciones vas a resolver?
Un sistema de 3 variables o un sistema de n variables?
Si tienes un sistema de n variables, entonces necesitas matrices dinámicas, y una función que te calcula el determinante de una matriz.
Por cual metodo quieres resolver el determinante?
Yo tengo un snippet de como encontrar el determinante de una mátriz por el método montante, lo hice en mi tiempo libre hace ya como 6 años.
Voy a asumir, que sabes que tienes que llenar una matriz y declararla, es un codigo poco eficiente, pero almenos cumple con el trabajo.
Código C++:
Ver original/* Seleccionamos Pivote */
pivote = montante[0][0];
cout<<endl<<"Pivote seleccionado: "<<pivote<<endl<<endl;
for (i=0;i<=1;i++){
for(j=0;j<=1;j++){
matriz[i][j] = 0;
}
}
for(i=1;i<=2;i++){
for(j=1;j<=2;j++){
matriz[i-1][j-1] = pivote*montante[i][j] - montante[0][j]*montante[i][0];
}
}
/* Imprimimos la nueva matriz */
j=0;
cout<<"Paso 2. Reducimos la matriz de orden 3 a orden 2 por el Metodo Montante"<<endl<<endl;
do{
for(i=0;i<=1;i++){
cout<<matriz[i][j]<<"\t";
}
cout<<endl;
j++;
}while(j<=1);
cout<<"--------------"<<endl<<" "<<pivote<<"\t Dividimos la matriz entre el pivote anterior"<<endl;
cout<<endl<<endl;
/*Calculamos determinante de una matriz de orden 2
(Diagonal principal - Diagonal Secundaria) / Pivote anterior */
cout<<"[("<<matriz[0][0]<<")("<<matriz[1][1]<<") - ("<<matriz[0][1]<<")("<<matriz[1][0]<<")]"<<endl<<"------------------"<<endl<<" "<<pivote<<endl;
det = (matriz[0][0]*matriz[1][1] - matriz[0][1]*matriz[1][0])/pivote;
cout<<endl<<"El valor del determinante es: "<<det<<endl<<endl;