%Secant uses the secant method to approximate the zero of an inline
%function with an interval given

f=inline('816*x^3-3835*x^2+6000*x-3125');

a=1;
b=2;

while abs(b-a) > eps*abs(b)
    c=a;
    a=b;
    
    b=b+(b-c)/(f(c)/f(b)-1);