function [] = sqrts(x)

for i = 7:25,
    down=sqsqrt(x,i);
    up=sqrtsq(x,i);
    fprintf('%d, %1.15g, %1.15g\n', i, down, up)
    %fprintf('%d, %1.15g, %1.15g\n', i, (down-x)/x, (up-x)/x)
end

