y'=( sin(y)^2 - x^2)/(x sin(2y)), y(1)=Pi/3