
x = roots([1 1 1 0 -1])
y = x.^2

x0 = [1/4;1/4];

[sol,it] = newtonrd('Fxy','dFxy',x0,1e-10,30);
if it > -1
    display('numero de iteraciones para obtener precision deseada')
    it
    display('aproximacion a cero de f')
    sol
else
    display('metodo no obtuvo precision deseada en maximo de iteraciones dado')
end
