#include <stdlib.h>
#include <string.h>
#include "gauss.h"

void
gauss_elimination( double *Matriz, double *Vetor, double *A,int t )
{
    int i,j,k;
    double m,s=0,X[t*t],Y[t];
    memcpy(X,Matriz,sizeof(double)*t*t);
    memcpy(Y,Vetor,sizeof(double)*t);
    //Escalonamento do sistema
    for(k=0;k<t-1;k++){
        for(i=k+1;i<t;i++){
            m=X[(i*t)+k]/X[(k*t)+k];
            for(j=0;j<t;j++){
                X[(i*t)+j]=X[(i*t)+j]-(m*X[(k*t)+j]);
            }
            Y[i]=Y[i]-(m*Y[k]);
        }
    }
    //Fim do escalonamento
    //Resposta Final
    A[t-1]=Y[t-1]/X[((t-1)*t)+t-1];
    s=0;
    for(k=t-2;k>=0;k--){
        s=0;
        i=0;
        for(j=t-1;j>=k+1;j--){
            i=i+1;
            s=s+(X[(k*t)+i+k]*A[i+k]);
        }
        A[k]=(Y[k]-s)/X[(k*t)+k];
    }
}