Anonymous

#include <stdlib.h>#include <string.h>#include "gauss.h"voidgauss_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];    }}