MODE 32: OFF

R%=500 
FOR X%=0 TO 1600 STEP 8
 FOR Y%=0 TO 1200 STEP 8
  V% = SQR( (X%-800)^2+(Y%-600)^2 )
  IF V%<500 THEN
   V% = (1-(V%/R%)^2) * FNrnd * 256
   GCOL V%, V%, V%
  ELSE
   GCOL FNrnd * 256, FNrnd * 256, FNrnd * 256
  ENDIF
  RECTANGLE FILL X%, Y%, 8
 NEXT
NEXT

PRINT "Press any key to exit"
A%=GET
QUIT

END

DEF FNrnd: IF X%>800: =RND(1) ELSE =FNSuperMegaNutter_b7

DEF FNf(v)=v-INTv

REM SuperMegaNutter, version "b7".

DEF FNSuperMegaNutter_b7
nn1 = nn1 + 0.499743730631690299468026764666361
IF nn1 > 2147483647.0 THEN nn1 = FNf( nn1 )
nn2 = nn2 - 0.619712099029093592809216927916784 + 0.01 * FNf( nn1 )
IF nn2 < 0.0 THEN nn2 = nn2 + 2147483647.0
LOCAL a,b
a = FNf( nn1 * FNf( nn1 * 0.70635556640556476940272083033647 ) * FNf( nn1 * 0.96241840600902081782173280676402 ) )
b = FNf( nn2 * FNf( nn2 * 0.41059769134948776074925592402104 ) * FNf( nn2 * 0.83598420020319692166434447425947 ) )
nn1 -= FNf( a + b ) * 0.249871861
= FNf( a * 2.79218049951371579168124556907984 - a + b + FNf( nn1 + nn2 ) * 1.1235804235832785923589346329 )