y'=sqrt(y^4/2-y^2+3/2), y(0)=1