System.out.println("Testing StepperDopr853");
Bessjy bess = new Bessjy();
for (i=0;i<nvar;i++) {
y[i]=bess.jn(i,x1);
yexp[i]=bess.jn(i,x2);
}
Output out=new Output(20);
rhs_StepperDopr853 d = new rhs_StepperDopr853();
StepperDopr853 s = new StepperDopr853();
Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);