y[i]=bess.jn(i,x1);
yexp[i]=bess.jn(i,x2);
}
Output out = new Output(20);
rhs_StepperDopr5 d = new rhs_StepperDopr5();
StepperDopr5 s = new StepperDopr5();
Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
ode.integrate();
for (i=0;i<nvar;i++) {
yout[i]=out.ysave[i][out.count-1];