// Gaussian Elemination
// 
// Copyright by mpC team, 1999
//
//    http://www.ispras.ru/~cbr
//    mpc@ispras.ru
//
// Version with fixed matrix size

#include<stdio.h>

#define N 5
#define VALUE1 10.
#define VALUE2 VALUE1*N

double A[N][N];

void init_matrix() {
  int i,j;
  
  for(i=0;i<N;i++)
    for(j=0;j<N;j++)
      A[i][j]=10.;

  for(i=0;i<N;i++) 
     A[i][i]=10. * N;
}


void print_matrix()
{
  int i,j;
  for(i=0;i<N;i++) {
    for(j=0;j<N;j++) printf("%lf\t",A[i][j]);
    printf("\n");
  }
  
}

main()
{
  int i,j;

  init_matrix();
  printf("\n Source Matrix :\n");
  print_matrix();
  for(i=0;i<N-1;i++) {
    for(j=i+1;j<N;j++) {
        double t;
	t=A[j][i]/A[i][i];
        if(A[j][i]!=0) A[j][i:N-1]-=t*A[i][i:N-1];
    }
  }
    

    printf("\n Matrix after LU transformation: \n");
    print_matrix();
    
}