( 8 July 88 -- MFB )
( this is identical to the algorithm implemented in CIR.ASM
   but it uses double precision FORTH words)

off printload
forget new
: new ; 

4 block x
4 block y
2 block x1
2 block y1 

,20 ,const aa
,100 ,const rad
140h const xc 
afh const yc 

: fred egam  ,0 x ,! rad y ,!  150 0 do 
       x ,@ y ,@ aa ,/ ,- ,dup x ,! single xc + x1 !
       x ,@ aa ,/ y ,@ ,+ ,dup y ,! single yc + y1 !

       x1 @ y1 @ 3 egadot loop ;

( traceon fred
 3 tracel )
fred