xrange=[0 0.5];
y0=1;
ode23(@myproblem,xrange,y0)