void dgeev_ ( char *jobvl, char *jobvr, int *n, double *A, int *lda, 
  double *wr, double *wi, double *vl, int *ldvl, double *vr,  int *ldvr, 
  double *work, int *lwork, int *info );

double complex *eigs ( int n, double *A );