Topic 10.6: Müller's Method (Maple)

Contents Previous Chapter Start of Chapter Previous Topic Introduction Notes Theory HOWTO Examples Engineering Error Questions Matlab Maple Next Topic Next Chapter

Implementing Müller's method in Maple is not that difficult:

To be completed...

> eps_step := 1e-5;
> eps_abs := 1e-5;
> p := x -> x^4 + 2*x^3 + 3*x^2 + 4*x - 7;
> x = <0, -0.1, -0.2>;
> y = map( p, x );
> for i from 1 to 100 do
     V = vander( x - x(2) );
     c = V \ y;
     x = ;

     if abs( x[2] - x[3] ) < eps_step and abs( y[3] ) < eps_abs  then
	break;
     elif ( i == 100 )
        error ( 'Muller\'s method failed to converge' );
     end if;
 end do;
>> x(3)

Copyright ©2005 by Douglas Wilhelm Harder. All rights reserved.