Gauss Map

 投稿者:しばっち  投稿日:2015年10月28日(水)22時15分13秒
 
!'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
 

戻る