siecz

function [wyn] = met_siecznych(g,a,b)
f = @(x)eval(g);
fb = f(b);
fa=f(a);
ep = 0.001;
if(f(a)*f(b)>0)
   disp('punkty  nie spelniaja zalozen')
   return;
end
while(abs(a-b)>0)
    c = a - fa*(a-b)/(fa-fb);
    fc = f(c);
    if(abs(fc)<ep)
        break;
    end
    a = b;
    fa = fb;
    b = c;
    fb = fc;
end
wyn = [b,c];
end