Python绘制Gumowski-Mira分形图代码

以下是绘制Gumowski-Mira分形图的Python代码：

import numpy as np
import matplotlib.pyplot as plt

def Gumowski_Mira(a, b, x, N):
    y = np.zeros(N)
    y[0] = x
    for i in range(1, N):
        y[i] = a * y[i-1] + 2 * b * y[i-1] / (1 + y[i-1]**2) + (1 - a - b) * y[i-1]**2 / (1 + y[i-1]**2)
    return y

a = 0.999
b = 0.001
N = 100000
x = np.zeros(N)
x[0] = 0.1
y = Gumowski_Mira(a, b, x[0], N)

plt.plot(range(N), y, 'b', linewidth=0.5)
plt.axis('off')
plt.show()

在这个代码中，Gumowski_Mira函数实现了Gumowski-Mira公式，并返回一个数组y，表示生成的分形图的每个点的y坐标。a和b是Gumowski-Mira公式中的参数，x是分形图的起始点，N是生成的点的数量。

在主程序中，我们给定了a、b、N和起始点x，然后调用Gumowski_Mira函数得到分形图的y坐标数组y，最后使用Matplotlib绘制出分形图。 plt.axis('off')命令关闭了坐标轴。