From: Richard F. <fa...@gm...> - 2023-07-19 15:49:03
|
> > >> >> f : sin(x)^4 * cos(x)^2; >> f1 : expand(exponentialize(f)); >> f2 : subst([x=1/%i, %e=E], f1); >> f3 : rat(f2, E); >> >> converts f to a rational function on E. What do I need to use instead >> of the "rat" in the last step to convert f to a _Laurent polynomial_ >> on E instead of a _rational function_ on E? >> >>> >>>> powerdisp:true; ratexpand(f3) produces this display (it looks better in 2-D)... 1/16+1/(64*E^6)-1/(32*E^4)-1/(64*E^2)-E^2/64-E^4/32+E^6/64. It takes the E^0 term out front, which may not fit the pattern you are looking for. or just powerdisp:true; ratexpand(ratsubst(E, exp(%i*x),exponentialize(f))); I am not aware, offhand, of any advantage of this representation over the original form. RJF |