#define DIM	3

void  ludcmp(double **a, int n, int *indx, double *d);
double *vector(long size);
void free_vector(double *v);
double **matrix(long rows, long cols);
void free_matrix(double **m);
int *ivector(long nh);
void free_ivector(int *v);