void  jacobi(double **a, int n, double *d, double **v, int *nrot);

double *vector(long size);
void free_vector(double *v);
double **matrix(long rows, long cols);
void free_matrix(double **m);

//int fileopen(void);