Untitled

S=20;
x0=2;
radq=babilonianRoots(S,x0)
disp(['Il valore della radice quadrata è', num2str(radq)])

function y=babilonianRoots(S,x0)
x=x0;
xnew=(x+S/x)/2;
counter=0;

while abs(xnew-x)>0.01
    x=xnew;
    xnew=(x+S/x)/2;
    counter=counter+1;
end
y=x
end