/* This is the header file for lu_decomp.c*/

#ifndef lu_decomp_h
#define lu_decomp_h

#include <default.h>
#include <matrix.h>

#define TINY 1.0e-40


void ludcmp(double **a, int n, int *indx, double *d);
void lubksb(double **a, int n, int *indx, double b[]);

double **mat_inverse(double **a, int n);
double mat_det(double **a, int n);


#endif