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