%newtona
h=0.01;
pochodna3=@(x) ((funkcja(x+h)-funkcja(x))/(2*h));
pochodna22=@(x) ((funkcja(x+2*h)-2*funkcja(x)+funkcja(x-2*h))/(4*h^2));
a=0;
b = 6;
liczbaIteracji=0;
x0 = b;
while(1)
    liczbaIteracji=liczbaIteracji+1;
    xn=(x0-(pochodna3(x0)/pochodna22(x0)));
    if(abs(xn-x0)<epsilon)
        break;
    end
    x0 = xn;
end
newtona=xn
liczbaIteracji