用python写程序，利用scipy.integrate.solve_bvp求解边值问题 y" + sin(x)y'+y·exp(x)=x^2, y(0)= 0,y(5)= 3

好的，我可以帮你写一个用 `scipy.integrate.solve_bvp` 求解该边值问题的 Python 程序。

import numpy as np
from scipy.integrate import solve_bvp

def fun(x, y):
    return np.vstack((y[1], -np.sin(x)*y[1]-y[0]*np.exp(x)+x**2))

def bc(ya, yb):
    return np.array([ya[0], yb[0]-3])

x = np.linspace(0, 5, 101)
y = np.zeros((2, x.size))

sol = solve_bvp(fun, bc, x, y)

if sol.success:
    import matplotlib.pyplot as plt
    xx = np.linspace(0, 5, 1001)
    yy = sol.sol(xx)[0]
    plt.plot(xx, yy)
    plt.xlabel('x')
    plt.ylabel('y')
    plt.show()
else:
    print('求解失败：', sol.message)

该程序使用 `scipy.integrate.solve_bvp` 求解边值问题：

y" + sin(x) y' + y exp(x) = x^2, y(0) = 0, y(5) = 3

其中 `fun` 函数是微分方程组的右端函数，`bc` 函数是边界条件函数。`x` 是自变量的取值范围，`y` 是解的初值矩阵。程序运行后，会输出求解结果的图像。