dsolve(diff(diff(y(x),x),x) = -(-a*cos(x)^2*sin(x)^2-m*(m-1)*sin(x)^2-n*(n-1)*cos(x)^2)/cos(x)^2/sin(x)^2*y(x),y(x))