callfunc sqrt{
while (@x*@x < @IST){
    set @x, (@x+(@IST/@x))/2;
}
return(@x);
}