読者です 読者をやめる 読者になる 読者になる

ブッダブロのような何か

言語:mathematica

random[] := 4*RandomComplex[] - 2 - 2 I

f[z_, c_] := z^2 + c

onestep[c_, limit_] := 
 NestWhileList[f[#, c] &, 0, Abs[#] < 2 &, 1, limit]

c2pair[z_] := {Re[z], Im[z]}

step[limit_, max_] := onestep[#, limit] & /@ Table[random[], {max}]

makelist[list_] := c2pair /@ Flatten[list]

budda[limit_, max_] := ListPlot[makelist[step[limit, max]]]


f:id:mikuwaorenoyome:20170329104843p:plain