Dear Wolfgang

yes I tried the ODE examples of the book, they all work fine.But they are all for second order ODE's. I have a few first order ODE's I would like to splve with ODE/ODEsolver.  Based on the book Open Source Physics (a User Guide with Examples) and the examples csm book I found out how to use ODe/ODEsolver for first order ODE's. Please have a look on the code below. May be you have some comments for improvements.

Thanks for your help

Martin

import org.opensourcephysics.numerics.ODE;
public class NucDecODE implements ODE{
double [] state =new double[3];
double lambda=0.01;
public  NucDecODE(double x){
state[0]=x;
state[1]=-lambda*state[0];
state[2]=0;
}
public double[] getState(){
return state;
}
public void getRate(double[] state,double[] rate){
rate[0]=state[1];
rate[1]=0;
rate[2]=1;
}
}
public class NucDecODEApp {
public NucDecODEApp(double x0) {
}
public static void main(String[] args) {
double dt=0.1;
double x0=100.0;
NucDecODE ode = new NucDecODE(x0);
ODESolver ode_solver=new Euler(ode);
ode_solver.initialize(dt);
while(ode.state[0]>90.0){
ode_solver.step();
}
System.out.println("final time= " +ode.state[2]);
System.out.println("final x =" +ode.state[0]);
}
}