/* ---------------------------------------- integrate.h -------------- */

#ifndef __INTEGRATE_H__
#define __INTEGRATE_H__

/* Define parameters for integration.
 * 
 * 
 */

void integrate(double zstart, double ystart[], int nv,
	       int *nout, double *zout, double **yout);
void jacobn(double z, double *y, double *dfdz, double **fjac, int nv);
void derivs(double z, double y[], double fvec[]);
void stifbs(double y[], double dydx[], int nv, double *xx, double htry,
	    double eps, double yscal[], double *hdid, double *hnext,
	    void (*drvs)(double , double [], double []));
void odeint(double ystart[], int nvar, double x1, double x2, double eps,
	    double h1, double hmin, int *nok, int *nbad,
 	    void (*drvs)(double, double [], double []));
void simpr(double y[], double dydx[], double dfdx[], double **dfdy,
	   int n, double xs, double htot, int nstep, double yout[],
	   void (*drvs)(double, double [], double []));
void pzextr(int iest, double xest, double yest[], double yz[], double dy[],
	    int nv);
void rzextr(int iest, double xest, double yest[], double yz[], double dy[],
            int nv); 
void mmid(double y[], double dydx[], int nvar, double xs, double htot,
	  int nstep, double yout[], void (*derivs)(double, double[], double[]));

void stiff(double y[], double dydx[], int n, double *x, double htry, 
	   double eps, double yscal[], double *hdid, double *hnext,
	   void (*derivs)(double, double [], double []));

#endif /* !__INTEGRATE_H__ */

/* ------- end ---------------------------- integrate.h  ------------- */