dsolve(diff(diff(y(x),x),x) = -a*(n-1)*sin(2*a*x)/cos(a*x)^2*diff(y(x),x)-n*a^2*((n-1)*sin(a*x)^2+cos(a*x)^2)/cos(a*x)^2*y(x),y(x))