帮我找绘制Gumowski-Mira 公式的分形图的python代码

以下是绘制Gumowski-Mira公式的分形图的Python代码。

import numpy as np
import matplotlib.pyplot as plt

def gumowski_mira(x, a, b, c, d):
    return (b * x + a * x**3) / (1 + c * x**2) + d * x

def generate_gumowski_mira(a, b, c, d, x0, n):
    x = np.zeros(n)
    x[0] = x0
    for i in range(1, n):
        x[i] = gumowski_mira(x[i-1], a, b, c, d)
    return x

def plot_gumowski_mira(a, b, c, d, x0, n):
    x = generate_gumowski_mira(a, b, c, d, x0, n)
    plt.plot(x[:-1], x[1:], ',')
    plt.title('Gumowski-Mira Fractal')
    plt.xlabel('x')
    plt.ylabel('y')
    plt.show()

plot_gumowski_mira(0.99, 0.002, 0.1, 0.1, 0.1, 100000)

在这个例子中，我们使用Gumowski-Mira公式生成一个分形图，其中a、b、c和d是参数，x0是初始值，n是生成的数据点数。我们使用generate_gumowski_mira函数生成数据，并使用plot_gumowski_mira函数绘制分形图。在这个例子中，我们使用了以下参数：a=0.99，b=0.002，c=0.1，d=0.1，x0=0.1，n=100000。您可以尝试不同的参数值来生成不同的分形图