function y = Newton_1d_modificado(der1f, der2f, der3f, x)
x1 = x;
x2 = x;
for iter = 1:10
    delta = der1f(x1)*der2f(x1)/(der2f(x1)^2 - der1f(x1)*der3f(x1));
    x2 = x1 - delta;
    disp([iter x2]);
    if abs(delta) < 1e-6 break;
    else
        x1 =x2;
    end;
end;
y = x2;