fx=f1a*(e+c2 e2+c3 e3+c4 e4); f1x=D[fx,e]; |
ey=e-Series[fx/f1x,
{
e,0,8
}
]; |
fy=f1a*(ey+c2 ey2+c3 ey3+c4 ey4); |
t=Series[fy/fx,
{
e,0,8
}
]; |
ez=ey-Series[1/(1−2t)*fy/f1x,
{
e,0,8
}
]; |
fz=f1a*(ez+c2 ez2+c3 ez3); |
s=Series[fz/fy,
{
e,0,8
}
]; |
v=Series[fz/fx,
{
e,0,8
}
]; |
gt=t0+t10*t+t20*t2/2+t30*t3/6; |
gs=s0+s10*s+s20*s2/2; |
gv=v0+v10*v+v20*v2/2; |
f1z=f1x*gt*gs*gv; |
e1=ez-Series [fz/f1z,
{
e,0,8
}
]//Simplify |
C1: “Out[a4](-1+t0 s0 v0)/(t0 s0 v0)” |
t0=1; s0=1; v0=1; a5=Coefficient[e1,e5] //Simplify |
C2: “Outa5(2 + t10)” |
t10=−2; a6=Coefficient[e1,e6]//Simplify |
C3: “Outa6 (1+s10) + (4 + t20 + 2v10)” |
s10=−1; t20=-2; a7=Coefficient[e1,e7]//Simplify |
C4: “Outa7 v10t30v10” |
v10=−2; t30=0; a8=Coefficient[e1,e8]//Simplify |
C5: “Outa8” |