代码是错误的，请正确表达M和H的关系式，并重新给我答案

对不起，您提到的M和H的关系应该是根据Stoner-Wohlfarth模型来计算的，这里是一个修正后的版本：

import numpy as np
import matplotlib.pyplot as plt

# 定义模型参数
gamma = 45  # 易轴与M的夹角 (假设固定)
Ms = 8e5   # 磁饱和强度
k = Ms / (2 * np.pi)  # Stoner-Wohlfarth常数

# 创建sita角度数组
sita_angles = np.arange(0, 91, 10)  # 变化的角度，步长为10度

# 正确的磁化强度M与外磁场H的关系
def M(H, sita):
    return Ms * (np.cos(sita) + k * np.sin(sita)**2)

# 计算磁滞回线数据点
hysteresis_data = [(np.sin(sita), M(np.cos(sita), sita)) for sita in sita_angles]

# 绘制磁滞回线
plt.figure(figsize=(6, 4))
for H, M_s in hysteresis_data:
    plt.plot(H, M_s, color='blue', linewidth=1)
plt.xlabel('磁感应强度 (H)')
plt.ylabel('磁化强度 (M)')
plt.title(f"磁滞回线 (gamma={gamma}°)")
plt.grid(True)
plt.show()