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