central case (b=0)

First poly(x)-a0 is set in foo1. The constant term a0 can be removed from poly in the following argument, but only be added in logexppoly later.

[Graphics:../Images/index_gr_146.gif]
[Graphics:../Images/index_gr_147.gif]

Make rule "rulepoly" to substitute "poly" for foo1.

[Graphics:../Images/index_gr_148.gif]

Considering

[Graphics:../Images/index_gr_149.gif]
[Graphics:../Images/index_gr_150.gif]

we can calculate [Graphics:../Images/index_gr_151.gif] as follows, and stored in foo4 up to [Graphics:../Images/index_gr_152.gif] terms.

[Graphics:../Images/index_gr_153.gif]
[Graphics:../Images/index_gr_154.gif]
[Graphics:../Images/index_gr_155.gif]

Then, we take the expectation of foo4, calculating [Graphics:../Images/index_gr_156.gif], and stored in foo5 below.

[Graphics:../Images/index_gr_157.gif]
[Graphics:../Images/index_gr_158.gif]
[Graphics:../Images/index_gr_159.gif]
[Graphics:../Images/index_gr_160.gif]
[Graphics:../Images/index_gr_161.gif]
[Graphics:../Images/index_gr_162.gif]

Considering

[Graphics:../Images/index_gr_163.gif]
[Graphics:../Images/index_gr_164.gif]

we take the log of foo5, and stored in foo8 up to [Graphics:../Images/index_gr_165.gif] terms below.

[Graphics:../Images/index_gr_166.gif]
[Graphics:../Images/index_gr_167.gif]
[Graphics:../Images/index_gr_168.gif]
[Graphics:../Images/index_gr_169.gif]

This is logeexppoly for b=0.

[Graphics:../Images/index_gr_170.gif]
[Graphics:../Images/index_gr_171.gif]

The following expression may be easier to read for us, but violating the summation convention rule of subscripts.

[Graphics:../Images/index_gr_172.gif]
[Graphics:../Images/index_gr_173.gif]


Converted by Mathematica      July 21, 2003