連續內接正方形

import numpy as np
import matplotlib.pyplot as plt
plt.figure(figsize=(5,5),dpi=90)
R=1  
N=4  
n=8
a=(2*np.pi)/(12.0)  #傾斜程度

for i in range(1,n+1,1):
    t = np.linspace(0.25*np.pi+a*(i-1),2.25*np.pi+a*(i-1),N+1)
    x=R*np.cos(t)
    y=R*np.sin(t)
    plt.plot(x,y,color="green",linewidth=3)
    R=R*0.7320508075688773
plt.xlim(-1, 1)
plt.ylim(-1, 1)
plt.show()