Yes I have an example m function below. I run it in freemat, since I don't have license of Matlab

function [k,c] = bisection(f, a, b, Tol)

%

% bisection(f, a, b, Tol)

%

% Input:

% f - function given as a string

% a - lower bound

% b - upper bound

%

% Output:

% k - number of iterations performed

% c - root of the function

%

% Example:

% [k,c] = bisection( 'x^6-x-1', 1, 2, 0.00001 )

format long;

N= 1e4;

if nargin == 3

Tol = 1e-4;

end

f = inline(f);

if (f(a) * f(b) > 0) || (a >= b)

error('INCORRECT INPUT f(a)*f(b)>0 or a>=b');

end

k = 1;

x(k) = (a + b) / 2;

while ((k = Tol)

if f(x(k)) == 0

error([ 'bisection: condition f(' num2str(x(k)) ...

')~=0 didn''t apply' ]);

end

if (f(x(k)) * f(a)) < 0

b = x(k);

else

a = x(k);

end

k = k + 1;

x(k) = (a + b) / 2;

c=x(k);

end

Oh sorry, it is better use the attachment.

On MATLAB File Exchange you can find few implementations:

http://www.mathworks.com/matlabcentral/fileexchange/index?utf8=%E2%9C%93&term=bisection

اه

you may use the example in applied numerical method with matlab (Steven Chapra)

function [root,fx,ea,iter]=bisect(func,xl,xu,es,maxit,varargin)

% bisect: root location zeroes

% [root,fx,ea,iter]=bisect(func,xl,xu,es,maxit,p1,p2,...):

% uses bisection method to find the root of func

% input:

% func = name of function

% xl, xu = lower and upper guesses

% es = desired relative error (default = 0.0001%)

% maxit = maximum allowable iterations (default = 50)

% p1,p2,... = additional parameters used by func

% output:

% root = real root

% fx = function value at root

% ea = approximate relative error (%)

% iter = number of iterations

if nargin0,error('no sign change'),end

if nargin

% Simple and Easy way to under stand of Bisection Method in MATLAB

a=input('Enter function with right hand side zero:','s');

f=inline(a);

xl=input('Enter the first value of guess interval:') ;

xu=input('Enter the end value of guess interval:');

tol=input('Enter the allowed error:');

if f(xu)*f(xl)