xy' =3y−2x−2(xy−x^2)^1/2