#include <stdio.h>

void PideDatos(int *Dim, float Sist[][102]);
void EscribeDatos(int Dim, float Sist[][102]);

void ResuelveGauss(int Dim, float Sist[][102]);


int main(void)
{
    int C,Dimension;
    float Sistema[101][102];
    PideDatos(&Dimension,Sistema);
    printf("\n\n\nEl SISTEMA introducido es el siguiente: \n\n");
    EscribeDatos(Dimension,Sistema);
    ResuelveGauss(Dimension,Sistema);
    printf("\n\n\nLas soluciones son:\n");
    for(C=1;C<=Dimension;C++) printf("\n X%d=%f\n",C,Sistema[C][Dimension+1]);
    
    scanf("");
    return(0);
}


void PideDatos(int *Dim,float Sist[][102])
{
    int A,B;
    printf("\n\n ||RESUELVE SISTEMAS LINEALES DETERMINADOS POR GAUSS||");
    printf("\n\n\n Introduce el numero de incognitas:(menor que 100)");
    scanf("%d",&*Dim);
    printf("\n\n PASE A INTRODUCIR CADA COMPONENTE DEL SISTEMA (A|B):");
    printf("\n\n MATRIZ A:\n");
    for(A=1;A<=*Dim;A++) for(B=1;B<=*Dim;B++){
        printf("\n Termino A(%d,%d):",A,B); scanf("%f",&Sist[A][B]);}
    printf("\n\n\n VECTOR B:\n");
    for(A=1;A<=*Dim;A++){
        printf("\n Termino B(%d):",A);scanf("%f",&Sist[A][*Dim+1]);
    }}

void EscribeDatos(int Dim, float Sist[][102])
{
    int A,B;
    for(A=1;A<=Dim;A++){
        for(B=1;B<=(Dim+1);B++){
            printf("%7.2f",Sist[A][B]);
            if(B==Dim) printf("   |");}
        printf("\n");
    }}

void ResuelveGauss(int Dim, float Sist[][102])
{
    int NoCero,Col,C1,C2,A;
    float Pivote,V1;
    
    for(Col=1;Col<=Dim;Col++){
        NoCero=0;A=Col;
        while(NoCero==0){
            if(Sist[A][Col]!=0){
                NoCero=1;}
            else A++;}
        Pivote=Sist[A][Col];
        for(C1=1;C1<=(Dim+1);C1++){
            V1=Sist[A][C1];
            Sist[A][C1]=Sist[Col][C1];
            Sist[Col][C1]=V1/Pivote;}
        for(C2=Col+1;C2<=Dim;C2++){
            V1=Sist[C2][Col];
            for(C1=Col;C1<=(Dim+1);C1++){
                Sist[C2][C1]=Sist[C2][C1]-V1*Sist[Col][C1];}
    }}
    
    for(Col=Dim;Col>=1;Col--) for(C1=(Col-1);C1>=1;C1--){
        Sist[C1][Dim+1]=Sist[C1][Dim+1]-Sist[C1][Col]*Sist[Col][Dim+1];
        Sist[C1][Col]=0;
    }
}

#include <stdio.h>
#include <stdlib.h>

main()
{
	
	int correc=0;
	float a11, a12, a13, b1, a21, a22, a23, b2, a31, a32, a33, b3, d, d2, d3, D, x, x2, x3, xx, X1, y, y2, y3, yy, X2, z, z2, z3, zz, X3;	
		
			while(correc==0)
			{
				printf("Metodo de Cramer con tres ecuaciones\n");
				printf("\nPrimera ecuacion");
				printf("\nInserta el valor de x1: ");
				scanf("%f",&a11);
				printf("\nInserta el valor de x2: ");
				scanf("%f",&a12);
				printf("\nInserta el valor de x3: ");
				scanf("%f",&a13);
				printf("\nInserta el b1: ");
				scanf("%f",&b1);
	
				printf("\nSegunda ecuacion");
				printf("\nInserta el valor de x1: ");
				scanf("%f",&a21);
				printf("\nInserta el valor de x2: ");
				scanf("%f",&a22);
				printf("\nInserta el valor de x3: ");
				scanf("%f",&a23);
				printf("\nInserta el b2: ");
				scanf("%f",&b2);
    
    			printf("\nTercera ecuacion");
				printf("\nInserta el valor de x1: ");
				scanf("%f",&a31);
				printf("\nInserta el valor de x2: ");
				scanf("%f",&a32);
				printf("\nInserta el valor de x3: ");
				scanf("%f",&a33);
				printf("\nInserta el b3: ");
				scanf("%f",&b3);
			
				printf("\nResultado:");
				printf("\n\n%.fx1 ",a11); printf("%.fx2 ",a12); printf("%.fx3 =",a13); printf("%.f",b1);
				printf("\n%.fx1 ",a21); printf("%.fx2 ",a22); printf("%.fx3 =",a23); printf("%.f",b2);
				printf("\n%.fx1 ",a31); printf("%.fx2 ",a32); printf("%.fx3 =",a33); printf("%.f",b3);
	
				printf("\n\nSi sus datos son correctos presione 1, si no presione 0 para ingresarlos nuevamente: ");
				scanf("%d",&correc);
				system("cls");	
		    }
			//Resolviendo Delta
			d=(a22*a33 - a23*a32)*a11; 
			d2=(a21*a33 - a23*a31)*-a12;
			d3=(a21*a32 - a22*a31)*a13;
			D=d+d2+d3;
			printf("Resultado de Delta");
			printf("\nD= " "%.f",D);
			 
			//Resolviendo X1
		    x=(a22*a33 - a23*a32)*b1;
			x2=(b2*a33 - a23*b3)*-a12;
			x3=(b2*a32 - a22*b3)*a13;
			xx=x+x2+x3;
			X1=xx/D;
			printf("\n\nResultado de X1");
			printf("\nX1= " "%.f",X1);
			
			//Resolviendo X2
		    y=(b2*a33 - a23*b3)*a11;
			y2=(a21*a33 - a23*a31)*-b1;
			y3=(a21*b3 - b2*a31)*a13;
			yy=y+y2+y3;
			X2=yy/D;
			printf("\n\nResultado de X2");
			printf("\nX2= " "%.f",X2);
			
			//Resolviendo X3
		    z=(a22*b3 - b2*a32)*a11;
			z2=(a21*b3 - b2*a31)*-a12;
			z3=(a21*a32 - a22*a31)*b1; 
			zz=z+z2+z3;
			X3=zz/D;
			printf("\n\nResultado de X3");
			printf("\nX3= " "%.f",X3);	
			
		return 0;
}