yy[1]=y[0].im();
yy[2]=y[1].re();
yy[3]=y[1].im();
Hypderiv d = new Hypderiv(a,b,c,z0,dz);
Output out = new Output();
StepperBS s = new StepperBS();
Odeint ode =new Odeint(yy,0.0,1.0,atol,rtol,0.1,0.0,out,d,s);
ode.integrate();
y[0]=new Complex(yy[0],yy[1]);
return y[0];
}