//
proc main;
	echo edit('%.15e',sqrt(%1+0));
	echo u_sqrt(%1+0);
//	dec(20,15) x = u_sqrt(%1+0);
//	echo x;
	return $ERROR;
end proc;

function u_sqrt(b);
	x0=b;
//echo x0;
	if (x0 > 0) then
		if (x0>=4) then x0/=2; endif;
echo x0;
		while (d=((y=x0*x0)-b)/(x0+x0)) do
//echo d x0;
			x0-=d;
echo x0;
		end do;
		if (y>b) then x0--; endif;
	end if;
	return x0;
end func;
