#include  <stdio.h>
#include <string.h>
#include <unistd.h>
#include <math.h>

/* linear algebra */
void pdinv(double *cinv, double *coeff, int n) ;

/* numer recipes p 97 */
int choldc (double *a, int n, double p[]);
void cholsl (double *a, int n, double p[], double b[], double x[]);
void cholesky(double *cf, double *a, int n) ;
void pmat(double *mat, int n)   ;