dsolve(diff(diff(y(x),x),x)*y(x)-diff(y(x),x)^2+(tan(x)+cot(x))*y(x)*diff(y(x),x)+(cos(x)^2-n^2*cot(x)^2)*y(x)^2*ln(y(x))=0,y(x))