A continuación les dejo un algoritmo para resolver sistema de ecuaciones con 3 incógnitas, me hacen falta algunos ajustes de estética pero por ahora lo dejo asi, ahi les dejo los demas a su criterio. Estoy utilizando el metodo de cramer.
#include <iostream>
using namespace std;
int main()
{
float x,y,z,re,x2,y2,z2,re2,x3,y3,z3,re3,deter,v;
cout<<" SISTEMA DE 3 ECUACIONES CON 3 INCOGNITAS ";
cout<<"Primera Ecuacion"<<endl;
cout<<"Inserta volor de X: ";
cin>>x;
cout<<"inserta valor de Y: ";
cin>>y;
cout<<"inserta valor de Z: ";
cin>>z;
cout<<"inserta el resultado: ";
cin>>re;
cout<<endl<<endl;
cout<<"Segunda ecuacion Ecuacion"<<endl;
cout<<"Inserta volor de X: ";
cin>>x2;
cout<<"inserta valor de Y: ";
cin>>y2;
cout<<"inserta valor de Z: ";
cin>>z2;
cout<<"inserta el resultado: ";
cin>>re2;
cout<<endl<<endl;
cout<<"Tercera Ecuacion"<<endl;
cout<<"Inserta volor de X: ";
cin>>x3;
cout<<"inserta valor de Y: ";
cin>>y3;
cout<<"inserta valor de Z: ";
cin>>z3;
cout<<"inserta el resultado: ";
cin>>re3;
cout<<endl<<endl;
deter= x*y2*z3 +x2*y3*z +x3*y*z2 -z*y2*x3 -z2*y3*x -z3*y*x2 ;
cout<<"n determinate es: "<<deter<<endl;;
v= (re*y2*z3 +re2*y3*z +re3*y*z2 -z*y2*re3 -z2*y3*re -z3*y*re2)/deter ;
cout<<"n El valor de X es: "<<v<<endl;
v= (x*re2*z3 +x2*re3*z +x3*re*z2 -z*re2*x3 -z2*re3*x -z3*re*x2 )/deter ;
cout<<"n El valor de Y es: "<<v<<endl;
v= (x*y2*re3 +x2*y3*re +x3*y*re2 -re*y2*x3 -re2*y3*x -re3*y*x2)/deter ;
cout<<"n El valor de Z es: "<<v<<endl;
return 0;
}