Programma C++ - Gauss

#include <stdio.h>

#include <conio.h>

#include <stdlib.h>

#include <math.h>

const int SIZE=51;

typedef float matrix [SIZE][SIZE];

typedef float vector [SIZE];

void errore(int);

void presentazione ();

void leggi_dim (int&);

float modulo(float);

void leggi_mat (int,matrix,vector);

void Gauss (int,matrix,vector);

int find_max (int, matrix, int, int);

void change_lines (int,int,int,matrix,vector);

void solve (int,matrix,vector,vector);

void print_sol (int,vector);

main ()

{ presentazione ();

int n; matrix a; vector b,sol;

leggi_dim (n);

leggi_mat (n,a,b);

Gauss (n,a,b);

solve (n,a,b,sol);

print_sol (n,sol);

system ("PAUSE");

return 0;

}

void presentazione ()

{ printf ("* * * Questo programma esegue l'algoritmo dell'eliminazione di Gauss

* * *\n");

printf ("\n\n* * * * su una matrice non singolare di dimensione nxn * * *

*\n\n\n");

system ("PAUSE");

return;

}

void leggi_dim (int& dim)

{ printf ("\nQual e` la dimensione della matrice dei coefficineti associata al

sistema?");

printf("\nInserici la dimensione: ");

scanf ("%d",&dim);

while (dim<=0 || dim>=SIZE)

{ printf ("\nAttenzione!");

printf("\nLa dimensione inserita non e' permessa.\n Inserisci dinuovo la

dimensione: ");

scanf ("%d",&dim);

}

return;

}

void leggi_mat (int dim,matrix a,vector b)

{ int i,j;

printf ("\nInserisci i coefficienti della matrice:");

for (i=1;i<=dim;i++)