/*

Copyright 1999-2002 David Joiner and the Shodor Education Foundation, Inc.

Shodor general use license

This code is provided as is, and is offered with no guarantee.

This code may be used freely in any non commercial product. Any works
using this code must publicly give credit for this work to the Shodor
Education Foundation, and should provide in that credit a link to the
location of the code on the Computational Science Education Reference
Desk (http://www.shodor.org/cserd) as well as to the Shodor Education
Foundation (http://www.shodor.org) and to any other authors listed in
the copyright statement found within the code.

Copyright statements within the code may not be modified.

This license does not grant permission for use in any commercial product.

*/

#ifndef _DIFFEQ
#define _DIFFEQ

#include <math.h>
#include <stdlib.h>

class Diffeq{
  public:
   int neq;
   double *x;
   double time;

  public:
   double * der;
   double * x_old;
   double * der_old;
   double * k1;
   double * k2;
   double * k3;
   double * k4;

  public:
   Diffeq(int j);
   ~Diffeq();
   void updateEuler(double, void no_name(int,double,double*,double*));
   void updateIEuler(double, void no_name(int,double,double*,double*));
   void updateRKutta4(double, void no_name(int,double,double*,double*));
};

#endif