
function [r, zpar] = root(x, zpar)
    if zpar(x) == x
        r = x;
    else
        zpar(x) = root(zpar(x), zpar);
        r = zpar(x);
    end
end
