System.out.println("Testing StepperSie");
y[0]=1.0; y[1]=1.0; y[2]=0.0;
Output out=new Output(20);
rhs_StepperSie d = new rhs_StepperSie();
StepperSie s= new StepperSie();
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];