import com.imsl.math.*;

public class OdeRungeKuttaEx1 {
    public static void main(String args[]) throws com.imsl.IMSLException {
        OdeRungeKutta.Function fcn = new OdeRungeKutta.Function() {
            public void f(double t, double y[], double yprime[]) {
                yprime[0] = 2. * y[0] * (1-y[1]);
                yprime[1] = -y[1] * (1-y[0]);
            }
        };
        
        double y[] = {1,3};
        OdeRungeKutta q = new OdeRungeKutta(fcn);
        int	nsteps = 10;
        for (int k = 0;  k < nsteps;  k++) {
            q.solve(k, k+1, y);
        }
        System.out.println("Result = {"+y[0]+","+y[1]+"}");
    }
}