zad3

%%
i=100;
a=-1;
b=2;
x=linspace(a,b,i);

h=0.2;
x0=0.3;
e=10^(-3);

x1=x0-(funkcja1(x0,h))/(funkcja2(x0,h));
x2=x1;

for j=1:i
    x3=x2-(funkcja1(x2,h)/funkcja2(x2,h));
    if abs(x3-x2)<=e
        break;
    else
        x2=x3;
    end
end


figure(3);
plot(x,funkcja(x));hold on;
plot(x0,funkcja(x3),'X');


function [funkcja] = funkcja(x)
funkcja=-(1./((x-0.3).^2+0.01))-(1./(((x-0.9).^2)+0.04))+6;
end

function [funkcja1] = funkcja1(x,h)
funkcja1=(funkcja(x+h)-funkcja(x-h))/(2*h);
end

function [funkcja2] = funkcja2(x,h)
funkcja2=(funkcja(x+h)-2*funkcja(x)+funkcja(x-h))/(h^2);
end