babylon sqr

disp(babylonsqr(64, 10))

function guess=babylonsqr(x, k)
    guess = 1;
     for i=1:k
        guess = (guess + (x/guess))/2;
     end
end