Python在数学建模中的简单应用

1.方程求导

from __future__ import print_functionfrom __future__ import divisionimport numpy as npimport scipy as spimport scipy.miscdef f(x): return 2*x*x + 3*x + 1print(sp.misc.derivative(f, 2))

2.求不定积分

from __future__ import print_functionfrom __future__ import divisionimport numpy as npimport scipy as spimport scipy.integratef = lambda x : x**2print(sp.integrate.quad(f, 0, 2))print(sp.integrate.fixed_quad(f, 0, 2))

3.求解非线性方程组

from __future__ import print_functionfrom __future__ import divisionimport numpy as npimport scipy as spimport scipy.optimizedef f(x):    return [5*x[1] + 3, 4*x[0]*x[0], x[1]*x[2] - 1.5]ans = sp.optimize.fsolve(f, [0, 0, 0])print(ans)print(f(ans))

4.求解线性方程组

from __future__ import print_functionfrom __future__ import divisionimport numpy as npimport scipy as spimport matplotlib.pylab as pltimport scipy.linalga = np.array([[1, 3, 5], [2, 5, 1], [2, 3, 8]])b = np.array([10, 8, 3])print(sp.linalg.solve(a, b))# print(sp.linalg.inv(a).dot(b))