10 void eigen(
double *a, 
double *ev, 
int n)
 
   14     dp = (
double *)calloc(n, 
sizeof(
double));
 
void trnm(double *a, int n)
 
void eigen(double *a, double *eval, int n)
 
int qrevec(double *eval, double *evec, double *dp, int n)
 
void housev(double *a, double *d, double *ud, int n)