Re: complex line integral in octave

[Stephen Montgomery-Smith]

On 04/16/2013 08:35 AM, Urs Hackstein wrote:
> There was another typo, since I guess you meant
> f = (((((a6.*z.+a5).*z.+a4).*z.+a3).*z.+a2).*z.+a1).*z.+a0;
> fp = ((((6*a6.*z.+5*a5).*z.+4*a4).*z.+3*a3).*z.+2*a2).*z.+a1;
> 
> g = z.*(fp.+7.*z.^6)./(f.+z.^7) .*(b-a)

I wasn't sure what you were trying to do.

This is the way it is usually done.  Apparently it is more efficient in
how many arithmetic operations it uses.

f = ((((((z.+a6).*z.+a5).*z.+a4).*z.+a3).*z.+a2).*z.+a1).*z.+a0;
fp = (((((7*z.+6*a6).*z.+5*a5).*z.+4*a4).*z.+3*a3).*z.+2*a2).*z.+a1;

g = z.*fp./f .*(b-a)

http://en.wikipedia.org/wiki/Horner%27s_method