#define DIM	3

void  ludcmp(double **, int n, int *, double *);
void lubksb(double **, int n, int *, double []);
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);