styczne

function [wyn] = met_stycz(g,x0)
%('x^2+x-1',-1) % Miejsce zerowe: -1.6180 oraz 0.6180
syms x;
ep = 0.001
f = @(x)eval(g);
fp = diff(eval(g));

for i = 1:10
       f1 = subs(fp,x0);
       x1 = x0 - f(x0)/f1;
       if(abs(f(x0)-f(x1))<ep)
           break
       end
       x0 = x1;
end
wyn = eval(x0);
end