python求积分不用外部包_高等数学被Python轻松化解？数学：我这么不要面子吗？...

说起高数大家一定不陌生，那真的是一个字，难，难，难，重要的事情说三遍！不过有人居然用Python轻松解决了,相信大家一定很好奇吧!

问题背景

高等数学应用非常广，基本上涉及到函数的地方都要用到微积分，还有在几何方面也是如此，计算机的应用让我们能简单快速处理各种高等数学中的计算，比如极限、导数、积分、微分方程等的计算。

实验目的

使用 Python 通过计算与作图，加强对极限、导数、积分等概念的理解，并掌握它们计算方法，以及求微分方程和方程组解析解的方法。

实验原理与数学模型函数极限的求解讨论以及两个重要极限的验证。

导数概念和导数的几何意义，以及计算多元函数偏导数和全微分的方法。

一元函数积分学和多元函数积分学。

微分方程和方程组在有无初始条件的分析。

实验所用软件Python 3.7

NumPy 1.16.4

SymPy 1.4

Matplotlib 3.1.1

主要内容函数极限的求解和两个重要极限的探究；

导数、高阶导数以及隐函数、参数方程定义函数导数的求解，多元函数偏导数和全微分的求解；

计算定积分和不定积分以及重积分的方法；

求解微分方程以及方程组解析解的方法。

实验过程

1. 函数极限的求解和两个重要极限

在这个实验中我们通过对简单的函数进行单侧极限的求解，并且分析两个重要极限。

例 1：考虑函数

解：编写Python代码如下：

import matplotlib.pyplot as plt

import numpy as np

import sympy as sp

# 求函数 y=arctan(1/x) 的左右极限

x = sp.Symbol('x')

fr = sp.atan(1 / x)

xl = sp.limit(fr, x, 0, dir='-')

xr = sp.limit(fr, x, 0, dir='+')

print('%s 左极限是：%s' % (fr, xl))

print('%s 右极限是：%s' % (fr, xr))

# 绘制函数 y=arctan(1/x) 的图像

x = np.arange(-6, 6, 0.01)

y = np.arctan(1 / x)

plt.title('y=arctan(1/x)')

plt.plot(x, y)

plt.show()

运行代码输出结果和绘制图像：

atan(1/x) 左极限是：-pi/2

atan(1/x) 右极限是：pi/2

根据计算结果和绘制的图像分析求得出题中函数的左右极限分别为 -pi/2 和 pi/2 。

例 2：两个重要极限的验证。

解：编写Python代码如下：

import matplotlib.pyplot as plt

import numpy as np

import sympy as sp

# 分析两个重要极限

x = sp.Symbol('x')

f1 = sp.sin(x) / x

f2 = (1 + 1 / x) **