python用slove解微分方程

Question

请问这个该怎么解？

Answer 1

先贴一下前两题的代码，后面的正在写：

导入库

from scipy.integrate import odeint, solve_bvp, solve_ivp
import numpy as np
import matplotlib.pyplot as plt

第一题 

#定义求解的函数
def fun(t,y):
  return t+y*y

t_span = (0,0.5) #定义自变量t的取值范围，（t0, tf)
y0 = [1]  #初始值 y0, y(t0) = y0
t = np.linspace(0,0.5,100) #方程数值解的点

ans = solve_ivp(fun, t_span, y0, t_eval=t) #求解
plt.plot(ans.t, ans.y[0]) #ans.y是二维的

第二题 

#定义求解的函数
def fun(t,y):
  return np.e**t/y

t_span = (0, 5) #定义自变量t的取值范围，（t0, tf)
y0 = [2]  #初始值 y0, y(t0) = y0
t = np.linspace(0, 5, 1000) #方程数值解的点

ans = solve_ivp(fun, t_span, y0, t_eval=t) #求解
plt.plot(ans.t, ans.y[0]) #ans.y是二维的