double fastcrawl(double func(double*),int m,int *nk,double *p0,double *dp,
                 double *p,int geo,int finalcall);
double slowcrawl(double func(double*),int m,int *nk,double *p0,double *dp,
                 double *p,int geo,int finalcall);