Source: methods/rk38.h

/*  rk38.h - 
    This file is part of Kalamaris
    Copyright (C) 2000  Antonio Larrosa Jimenez
    Kalamaris' homepage : http://www.arrakis.es/~rlarrosa/kalamaris.html

    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.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program; if not, write to the Free Software
    Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 

    Send comments and bug fixes to Antonio Larrosa <larrosa@kde.org>

 ***************************************************************************/
#include "mapmethod.h"

class MethodRungeKutta38 : public MapMethod
{
protected:
  TMatrix *method(const TMatrix &f,const T &x, const TMatrix &y,
              QDict <Map> &vars, unsigned int k, unsigned int n,const T &h);

private:
  void calcH0(const TMatrix &f, const TMatrix &y, QDict <Map> &vars,
      unsigned int k, unsigned int n,const T &h, TMatrix &H, 
      const T &x0,const T &,const T &, TMatrix &fx0, TMatrix &, TMatrix &);
  void calcH1(const TMatrix &f, const TMatrix &y, QDict <Map> &vars,
      unsigned int k, unsigned int n,const T &h, TMatrix &H, 
      const T &,const T &x1,const T &, TMatrix &fx0, TMatrix &fx1, TMatrix &);
  void calcH2(const TMatrix &f, const TMatrix &y, QDict <Map> &vars,
      unsigned int k, unsigned int n,const T &h, TMatrix &H, 
      const T &,const T &,const T &x2, TMatrix &fx0, TMatrix &fx1, TMatrix &fx2);
  
public:
  MethodRungeKutta38(void);
  MethodRungeKutta38(const MethodRungeKutta38 &m);
  virtual ~MethodRungeKutta38();
  
  virtual Map *copy(void) const;
  
  Map *eval(const QDict <Map> &vars);

  virtual QString string(void) const;
  virtual QString help(void) const;

};