void amoeba(double **p, double y[], int ndim, double ftol,double (*funk)(double []), int *nfunk);
double amotry(double **p, double y[], double psum[], int ndim,double (*funk)(double []), int ihi, double fac);

double *vector(long size);
void free_vector(double *v);
double **matrix(long rsize, long csize);
void free_matrix(double **m);