x1=0.0;x2=1.0;
rhs_StepperStoerm d = new rhs_StepperStoerm();
StepperStoerm s1= new StepperStoerm();
Odeint ode1 = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s1);
ode1.integrate();
yexp[0]=bess.airy_ai(1.0);
yexp[1]=bess.airy_bi(1.0);
for (i=0;i<nvar;i++) {
yout[i]=out.ysave[i][out.count-1];
System.out.printf("%f %f\n", yout[i],yexp[i]);