Interval Halving (Bisection) Method

>>bisec(1,2,50)

function bisec(xa,xb,n)
fx='x^3+x^2-3*x-3';
x=xa; fa=eval(fx);
x=xb; fb=eval(fx);
eps=1.0e-6;
for i=1:n
    xc=(xa+xb)/2; x=xc; fc=eval(fx);
    X=[i,xa,xb,xc,fc];
    disp(X)
    if abs(fc) < eps
        return;
    else
    end
    if fc*fa <0
        xb=xc;
    else
        xa=xc;
    end
end



Newton's Method

>>newton(1,50)

function newton(x0,n)
fx='x^3+x^2-3*x-3'; dfdx='3*x^2+2*x-3';
eps=1.0e-6;
for i=1:n
    x=x0; f0=eval(fx); df0=eval(dfdx);
    X=[i,x0,f0,df0];
    disp(X)
    x0=x0-f0/df0
    if abs(f0) < eps
        return;
    else
    end
end