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