dsolve(diff(diff(y(x),x),x)-diff(y(x),x)*cot(x)+y(x)*sin(x)^2=0,y(x))