void cdgefa(double **a,int n, int ipvt[],int *info);
void cdgesl(double **a,int n,int ipvt[],double b[],int job);
void cdgeco(double **a,int n, int *ipvt, double *rcond,double *z);
void cdgedi(double **a,int n,int ipvt[],double det[],double work[],int job);
void cdpofa(double **a,int n,int *info);
void cdposl(double **a,int n,double *b);