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