% Un algoritmo para aproximar raices cuadradas (Ejmeplo 1.2.3)

clear 
clc

n = 7; % quiero calcular \sqrt{n}
y = 2; % estimacion inicial

lista = [y]; % quiero ir guardando las aproximaciones que vamos generando

tol = 1e-12;

dif = inf;

while abs(dif) > tol
    y = 0.5 * (y + n/y);
    dif = y - lista(end);
    lista = [lista; y];
end

y