/* this is the header file to solve ordinary differential equations */

#ifndef odeint_h
#define odeint_h

#include <default.h>
#include <nrutil.h>
#include <complex.h>
#include <matrix.h>

void odeint(double ystart[],int nvar,double x1,double x2,double eps,double h1,
	double hmin, int *nok, int *nbad,
	void (*derivs)(double, double [], double []),
	void (*rkqs)(double [],double [],int,double *,double,double,double [],
	double *, double *, void (*)(double, double [], double [])));
void hypdrv(double s, double yy[], double dyyds[]);
void bsstep(double y[],double dydx[],int nv,double *xx,double htry,double eps,
	    double yscal[], double *hdid, double *hnext,
	    void (*derivs)(double, double [], double []));
void mmid(double y[],double dydx[],int nvar,double xs,double htot,int nstep,
	  double yout[], void (*derivs)(double, double[], double[]));
void pzextr(int iest, double xest, double yest[], double yz[], double dy[], 
	    int nv);




#endif