$0=K+1      # $0,$1 := ((K+1)*N),(((K+1)*N)/(K*K+K+1))
$1=K*K
$1=$1+$0
$0=$0*N
$1=$0/$1
$2=N*N      # $2,$3 := (N*N),(K*N);
$3=K*N
while$0>$1  # while $0>$1:
  $0=$1     #   $0,$1 := $1,(($1*$1+N*N)/(2*$1+K*N));
  $4=$1*$1
  $1=2*$1
  $4=$4+$2
  $1=$1+$3
  $1=$4/$1
end         # end
$0=$0+N     # return ($0+N);
return$0