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