// crout.h  solve simultaneous equation  A x = y  by Crout method
//           A and y destroyed
//           N is dimensioned size of matrix, n is number of equations

void crout(int n, int N, double A[N][N], double y[], double x[]);