int ludcmp(float *m, int n, int indx[], float *d);
int lubksb(float *m, int n, int indx[], float b[]);