编写一个程序，要求用Simpson法则计算Jm(x)的值，Jm(x)函数为一个对o的积分，积分函数为cos(m*o-x*sin(o)),积分上下限分别为3.1416和0，取N=1000，在一幅图中分别画出m=0，m=1，m=2时的Jm(x)函数图像，x的取值范围是0到20.

你好，我可以为您提供一段Python代码，使用Simpson法则计算Jm(x)的值并绘制图像：

import numpy as np
import matplotlib.pyplot as plt

# 定义被积函数
def func(o, m, x):
    return np.cos(m*o - x*np.sin(o))

# 定义Simpson法则
def simpson(f, a, b, N):
    h = (b-a)/N
    x = np.linspace(a,b,N+1)
    y = f(x)
    return h/3 * (y[0] + 4*np.sum(y[1:N:2]) + 2*np.sum(y[2:N-1:2]) + y[N])

# 计算Jm(x)的值
def Jm(m, x):
    a = 0
    b = np.pi
    N = 1000
    return simpson(lambda o: func(o, m, x), a, b, N)

# 绘制图像
x = np.linspace(0,20,200)
for m in range(3):
    J = np.array([Jm(m, xi) for xi in x])
    plt.plot(x, J, label=f"m = {m}")
plt.legend()
plt.xlabel("x")
plt.ylabel("Jm(x)")
plt.show()

执行以上代码后，您将会看到绘制的Jm(x)函数图像。

另外，我理解您的要求，不会透露您的任何要求信息。