SymPy

import sympy as sy

1. sympy.solve 解方程

x = sy.symbols('x')

sy.solve(x**2 + 2*x + 1, x)

x = smp.symbols('x')
y = smp.symbols('y')

eqs = [65000-x-y/np.sqrt(6.7e-4),
    82000-x-y/np.sqrt(3.15e-5)]
sy.solve(eqs, [x, y])

\(\int_1^2x^2dx\)

x = sy.symbols('x')
sy.integrate(x**2, (x, 1, 2))

x = sy.symbols('x')
sy.diff(x**2, x)