Maple Program for Bisection Method

 

Bisection Method

restart;

 with(Student[NumericalAnalysis]):

f := x^3 + 4*x^2 - 10;

Bisection(f, x = [1, 2], tolerance = 10^(-4));

1.365112304

Bisection(f, x = [1, 2], tolerance = 10^(-4), output = sequence);

[1., 2.], [1., 1.500000000], [1.250000000, 1.500000000], [1.250000000, 1.375000000], [1.312500000, 1.375000000], [1.343750000, 1.375000000], [1.359375000, 1.375000000], [1.359375000, 1.367187500], [1.363281250, 1.367187500], [1.363281250, 1.365234375], [1.364257812, 1.365234375]

Bisection(f, x = [1, 2], tolerance = 10^(-2), stoppingcriterion = absolute);

1.367187500

Bisection(f, x = [1, 2], output = animation, tolerance = 10^(-3), stoppingcriterion = function_value);


Bisection(f, x = [1, 2], output = plot, tolerance = 10^(-3), maxiterations = 15, stoppingcriterion = relative);


Bisection(f, x = [1, 2], output = animation, tolerance = 10^(-3), maxiterations = 15, stoppingcriterion = relative);

Bisection(f, x = [1, 2], output = information, tolerance = 10^(-9), maxiterations = 15, stoppingcriterion = function_value);


 

Comments

Post a Comment

Popular posts from this blog

Maple Program for Secant Method

Image
Working Rules for Secant Method in Maple 2020  The Secant command numerically approximates the roots of an algebraic function, f(x) , using a technique similar to Newton's method but without the need to evaluate the derivative of f(x) . Given an expression f and an initial approximate a, the Secant command computes a sequence  " p[k],  k=0..n ", of approximations to a root of f(x)=0 , where "n" is the number of iterations taken to reach a stopping criterion. The Secant command is a shortcut for calling the Roots command with the method=secant option The criterion that the approximations must meet before discontinuing the iterations. The following describes each criterion: [> restart; [> with(Student[NumericalAnalysis]); [> f := x^3 + 4*x^2 - 10; [> Secant(f, x = [1, 2], tolerance = 10^(-2));                           1.365211903 [> Secant(f, x = [1, 2], tolerance = 10^(-2), stoppingcriterion =...
Read more