dsolve(diff(diff(y(x),x),x)-(1+2*tan(x)^2)*y(x)=0,y(x))