ptofijo

function [xs,nit,err] = puntofijo(gx, x0)
    xs = [];
    nit = [];
    err = [];
    syms x;

    g = gx+x;
    gxprima = diff(g);
    iteraciones = 0;
    if(abs(subs(gxprima, x0)) > 1)
        disp('La función diverge con el método de Punto Fijo.');
        return ;
    else
        error = 1;
        while (abs(error) > 0.00000001) %tolerancia 10e-8
            x1 = subs(g, x0);
            xs = [xs; x1]; %agregar cada raíz a un vector de raíces
            error = abs(x1-x0);
            err = [err; error]; %agregar cada error a un vector de errores
            iteraciones = iteraciones + 1;
            nit = [nit; iteraciones]; %agregar cada iteracion a un vector de iteraciones
        end
    end
end