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