// crout1.h  solve simultaneous equation  A x = y  by Crout method
//           A and y destroyed

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