脚本：

def f(x, y, d, R, n, N):
    if n < N:
        b = zeros(3)
        if n > 3:
            b = [.5,.2,.6]

        X = R * sin(np.deg2rad(d)) + x
        Y = R * cos(np.deg2rad(d)) + y
        plot(X, Y, color=b, linewidth=np.linalg.norm(R)/15)
        A = array([0,-48,-25,43,31])
        for i in range(1, 5):
            r = (i + 1) % 2
            L = .58 + r * .07
            f(X[i], Y[i], d+A[i], L*R, n+1, 7-r)

ax = axes(outerposition=[0,0,1,1], axis=False)
f(0, 0, array([0,0,-5.76,0,6.58,0]), array([0,55,99,55,105,55]), 0, 1)