Python+NumbaでWikipediaのList of chaotic mapsを１つずつ描く（その2） Tinkerbell Map。あのティンカーベルに似ている？

ネタがないのでWikipediaのList of Chaotic Mapsをいろいろ見ていくシリーズ。

https://en.wikipedia.org/wiki/List_of_chaotic_maps

今日はTinkerbell Mapだ。

https://en.wikipedia.org/wiki/Tinkerbell_map

x[i+1]=x[i]**2-y[i]**2+a*x[i]+b*y[i]
y[i+1]=2*x[i]*y[i]+c*x[i]+d*y[i]

こんな写像。これはとても簡単に

綺麗に描ける。うーん、ティンカーベルかー。

前にも書いたけど阪神高速のロゴに見える、、、

Pythonのコードはこちら：

import numpy as np
import matplotlib.pyplot as plt
from numba import jit

a=0.9
b=-0.6013
c=2.0
d=0.50

@jit
def tinker_calc(x,y,N):

    for i in range(N-1):
        x[i+1]=x[i]**2-y[i]**2+a*x[i]+b*y[i]
        y[i+1]=2*x[i]*y[i]+c*x[i]+d*y[i]

N=20000

fig=plt.figure(figsize=(10,10))
plt.xlim(-1.6,0.6)
plt.ylim(-1.6,0.6)

x=np.zeros(N)
y=np.zeros(N)
x[0],y[0]=-0.72,-0.64
tinker_calc(x,y,N)
plt.scatter(x,y,s=1)

plt.show()