import numpy as np
import matplotlib.pyplot as plt
plt.figure(figsize=(5,5),dpi=80)
def F(x0,y0,r,a): #定義遞迴程式
if r <1 :return
x1=x0+r*np.cos(a) ; y1=y0+r*np.sin(a)
X=[x0,x1]
Y=[y0,y1]
plt.plot(X,Y,color="blue",linewidth=2) #畫線
F(x1,y1,0.5*r,a-pi/3)
F(x1,y1,0.6*r,a+pi/3)
#以下內容為主程式
pi=3.14159 #圓周率
(x0,y0,r,a)=(0,0,200,pi/2) #(基準點之x,y坐標,長度,方向角)
F(x0,y0,r,a)
plt.show()