python求解微分方程代码
时间: 2023-10-06 18:13:54 浏览: 50
以下是使用Python求解微分方程的代码示例:
```python
from typing import List
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
# 定义微分方程
def diff_equation(y_list: List\[float\], x: float):
y1, y2, y3 = y_list
return np.array(\[y2 * y3, -y1 * y3, -0.51 * y1 * y2\])
# 设置求解范围
x = np.linspace(0, 12, 100)
# 求解微分方程
result = odeint(diff_equation, \[0, 1, 1\], x)
# 绘制图像
plt.plot(x, result\[:, 0\], label='y1')
plt.plot(x, result\[:, 1\], label='y2')
plt.plot(x, result\[:, 2\], label='y3')
plt.legend()
plt.grid()
plt.show()
```
这段代码使用了`odeint`函数来求解微分方程。首先定义了微分方程的函数`diff_equation`,然后使用`odeint`函数传入微分方程函数、初始值和求解范围,得到微分方程的数值解。最后使用`matplotlib`库绘制了结果的图像。
#### 引用[.reference_title]
- *1* *2* [Python解微分方程(验证数学建模第五版火箭发射模型)](https://blog.csdn.net/qq_48744162/article/details/121234882)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insert_down1,239^v3^insert_chatgpt"}} ] [.reference_item]
- *3* [python——求解微分方程](https://blog.csdn.net/Thattear/article/details/124308546)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insert_down1,239^v3^insert_chatgpt"}} ] [.reference_item]
[ .reference_list ]
相关推荐
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)