|
!'Gauss Map
OPTION ARITHMETIC COMPLEX
SET WINDOW -1,1,1,-1
SET POINT STYLE 1
FOR Y=-1 TO 1 STEP 1/512
FOR X=-1 TO 1 STEP 1/512
LET Z=COMPLEX(X,Y)
WHEN EXCEPTION IN
LET C=ABS(F(Z))*8
USE
END WHEN
SET POINT COLOR INT(C)
PLOT POINTS : X,Y
NEXT X
NEXT Y
END
EXTERNAL FUNCTION F(Z)
OPTION ARITHMETIC COMPLEX
LET F=COMPLEX(RE(1/Z)-INT(RE(1/Z)),IM(1/Z)-INT(IM(1/Z)))
END FUNCTION
|
|