尝试用任意曲线拟合Python中的multiple Moyal函数

import numpy as np
from scipy import optimize
import matplotlib.pyplot as plt

def multiple_moyal(x, a, b, c, d):
    return a * np.exp(-b * (x - c)**2) + d

x = np.array([1, 2, 3, 4, 5])
y = np.array([0.5, 0.8, 1.2, 1.5, 1.9])

params, params_covariance = optimize.curve_fit(multiple_moyal, x, y)

x_fit = np.linspace(0, 6, 100)
y_fit = multiple_moyal(x_fit, params[0], params[1], params[2], params[3])

plt.scatter(x, y, label='Data')
plt.plot(x_fit, y_fit, label='Fit')
plt.legend()
plt.show()