\ SQR 99-01-31 KBB

( Square root using [ a + b ]^2 method   KBB         Sept 22 1990 )

: (SQR) ( LO HI RESULT ... LO HI RESULT )
   8 0 DO
      >R D2* D2*
      R> 2* DUP >R
      2* 1+
      OVER SWAP -
      DUP 0< IF
         DROP R>
      ELSE
         NIP R> 1+
      THEN
   LOOP
   ;

: SQR ( N ... N1 ) ( return square root of single number )
   0 0 (SQR)
   NIP NIP
   ;

: DSQR ( D -- N ) ( return square root of double number )
   0 0 (SQR)
   ROT DROP (SQR)
   NIP NIP
   ;