codigo en c++q implemta el gauus yordan

lo hice bajo el interprete netbenas no tiene todas las validaciones necesarias pero trata de cumplir su prometido

// 

// File:   newmain.cc

// Author: felix

//

// Created on 13 de julio de 2008, 11:14 PM

//

#include

#include

double a1,b1,c1,r1,a2,b2,c2,r2,a3,b3,c3,r3;

void calcular () {

////////obtener la ecuasion

    cout<<"escribir la ecuasion 1 \n";

cin>>a1; cin>>b1; cin>>c1;cin>>r1;

cout<<"escribir la ecuasion 2 \n";

cin>>a2; cin>>b2; cin>>c2;cin>>r2;

cout<<"escribir la ecuasion 3 \n";

cin>>a3; cin>>b3; cin>>c3;cin>>r3;    

////////resolver la ecuasion

///uno en la primera fila

b1=b1/a1;

c1=c1/a1;

r1=r1/a1;

a1=a1/a1;

//primer cero de la columna 1

b2=((a2 * -1) * b1) + b2;

c2=((a2 * -1) * c1) + c2;

r2=((a2 * -1) * r1) + r2;

a2=0;

///segundo cero columna 1

b3=((a3 * -1) * b1) + b3;

c3=((a3 * -1) * c1) + c3;

r3=((a3 * -1) * r1) + r3;

a3=0;

//uno en la fila 2

c2=c2/b2;

r2=r2/b2;

b2=b2/b2;

///primer cero en la columna 2

c1=((b1 * -1) * c2) + c1;

r1=((b1 * -1) * r2) + r1;

b1=0;

///segundo cero en la columna 2

c3=((b3 * -1) * c2) + c3;

r3=((b3 * -1) * r2) + r3;

b3=0;

///uno en la fila #3

r3=r3/c3;

c3=c3/c3;

////primer cero en la columna 3

r1=((c1 * -1) * r3) + r1;

c1=0;

//segundo cero en la columna 3

r2=((c2 * -1) * r3) + r2;

c2=0;

//////////

cout<<"\n \n \n";

cout<<"r1 = "<

cout<<"\n r2 = "<

cout<<"\n r3 = "<

}

int ayuda (){

 int opc;   

 do{

     cout<<" si es la primera ves q usas esta aplicasion te convendria acceder a la ayuda, siendo de otra manera puedes empesar a trabajar sin ningun problema  \n \n";

    cout<<"ayuda (1)     ";

    cout<<"calcular(2)   ";

    cout<<"salir (3)   \n    ";

   cout<<"opcion  = "; cin>>opc;

 }while(opc>3 || opc<1);

 

 if  (opc==1){

                    

                     cout<<"este programa toma las costantes de la ecuasion de  3*3 y trata de deducir el valor de las variables, asi, si la ecausion es del tipo  \n";

                     cout<<"5a + 6b - 7c = 8  \n";

                    cout<<"1a - 2b + 3c= 4   \n ";

                    cout<<"-7a + 8b -9c =10 \n";

                    cout<<"5 seria a1, 1 seria a2 y -7 seria a3 asi sucesivamente al igual q 8 seria r1 hasta 10 que seria r3 ";

                    cout<<"ademas teniendo en cuenta que cuando pide la ecausion 1 lo que digitaremso seran las constantes de a1, b1, c1 y r1";                                         

                    }

 

  if  (opc==3){

                         return (EXIT_SUCCESS); 

                         

                      }

 

 if (opc==2){

                      calcular();

     

                    }

   

}

int main(int argc, char** argv) {

    

    

    cout<<"sistema gauss-jordan \n \n \n";

    

    ayuda();

    

    

     return (EXIT_SUCCESS);

}