/* 
   
   Copyright (C) Bernard Piette 
                 University of Durham (UK) 
                 Department of Mathematical Sciences
   
   This program is free software; you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation; either version 2 of the License, or
   (at your option) any later version.
   
   You can modify or use any part of the source code but only under the  
   condition that this Copyright notice is kept in the resulting source file.
*/
#ifndef PDE_RK4_H
#define PDE_RK4_H
#include <pde.h>

class pde_RK4 : public pde { 
public:
  pde_RK4(double xmin,double xmax,int nx);

  ~pde_RK4();
  virtual double integrate_1step();    
  void sub_RK4_1step(double *f,double *ef,double *mf,double *nef,
                     double Cnef,double Cmf);
protected:
  double *tmp1_;        // Same size as Fcts_
  double *tmp2_;        // Same size as Fcts_
  double *tmp3_;        // Same size as Fcts_
};
#endif