Sunday, 11 March 2007

ca.analysis and odes - How to invert this series?

The inverse series doesn't have that form. If z=eQ(1+Q1)z=eQ(1+Q1) then



Q=logz+log(1+Q1)=logz+(logz)1+O((logz)2)quadmboxaszto0+.Q=logz+log(1+Q1)=logz+(logz)1+O((logz)2)quadmboxaszto0+.



The general form should be



Q=logz+sumi>0bi(logz)i.Q=logz+sumi>0bi(logz)i.



You might be able to coerce this into the Lagrange inversion form, but I don't see how right now. I would just generalize the solution above:



Write
Q=logz+logleft(1+sumaiQiright)=logz+sumfrac(1)kkleft(sumaiQiright)k.Q=logz+logleft(1+sumaiQiright)=logz+sumfrac(1)kkleft(sumaiQiright)k.



Expanding this will give you a formula of the form
Q=logz+sumNi=1ciQi+O(QN1)quad()Q=logz+sumNi=1ciQi+O(QN1)quad()
for any NN you like. If you already know that Q=logz+sumN1i=0bi(logz)i+O((logz)N)Q=logz+sumN1i=0bi(logz)i+O((logz)N), then plug your known values into ()() to deduce the value of bNbN.

No comments:

Post a Comment