y0=.025;
[t,y]=ode23('f',[0 2],y0);
y0=.25;
[t1,y1]=ode23('f',[0 2],y0);
plot(t,y,'r',t1,y1,'c')