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