First, yes, the characteristic equation is [itex]r^2+ 1= 0[/itex] which has roots i and -i. That does NOT mean [itex]y(x)= C_1e^x cos(x)+ C_2e^x sin(x)[/itex]. That would correspond to roots 1+ i and 1- i. The general solution to the associated homogeneous equation is [itex]y(x)= C_1 cos(x)+ C_2 sin(x)[/itex]
Since "sec(x)" is NOT one of the functions that forms solutions for homogeneous linear equations with constant coefficients (polynomials, exponentials, sine and cosine, and combinations of those), there will be no simple way to use "undetermined coefficients".
But you can use "variation of parameters". Look for a solution of the form y(x)= u(x)cos(x)+ v(x)sin(x) where u(x) and v(x) are unknown functions. Notice that we don't need any special knowledge here. Given any solution, y, we could create many functions "u" and "v" that would work. From y= ucos(x)+ vsin(x), y'= u' cos(x)- u sin(x)+ v' sin(x)+ v cos(x). Because there are, as I said, many different function u and v we can simplify our "search" by also requiring that u' cos(x)+ v' sin(x)= 0. With that condition, y'= -u sin(x)+ v cos(x) and then y"= -u' sin(x)- u cos(x)+ v' cos(x)- v sin(x).
Putting those into y"+ y= sec3(x) we have (-u' sin(x)- u cos(x)+ v' cos(x)- v sin(x))+ (u cos(x)+ v sin(x))= -u' sin(x)+ v' cos(x)= sec3(x).
That, -u' sin(x)+ v'cos(x)= sec3(x), and u' cos(x)+ v' sin(x)= 0 give two equations that can be treated as ordinary algebraic equations to be solved for u' and v'. Then integrate to find u and v.