ODE

package odeadapt;

import Jama.*;

/** 
 * class to model ODEs <br>
 * x' = f(x, t), x(t0) = x0<br>
 * in d dimensions
 */
public abstract class ODE {

  /** initial time */
  protected double t0;
  /** initial values at time t0 */
  protected Matrix x0; 

  /**
   * right hand side, giving the ODE<br>
   *  x' = f(x, t)
   */
  public abstract Matrix f(Matrix x, double t);

  /**
   * construct an ODE with given initial conditions
   */
  public ODE(Matrix x0, double t0) {
    this.t0 = t0;
    this.x0 = x0;
  }
}