#include <math.h>
#include <stdio.h>


#define MNP    200
#define ALPHA  1.0
#define BETA   0.5
#define GAMMA  2.0
#define LW     5
#define ROOT2  1.414214


extern int M,NVPP,N, NSF, NPF;
extern int       *h, *l, *plist, *splst;
extern int       np, maxiter, niter;
extern double    *next, *mean, *error, *maxerr, *step, **simp;
extern double    **data;
extern FILE      *fpdata;

extern double    f();

extern double rwork1, rwork5, rwork6, rwork7;
extern double *atol, *rtol, t, tout; 
extern int iwork1, iwork2, iwork5, iwork6, iwork7, iwork8, iwork9;
extern int itol, itask, istate, iopt, jt;