Sympy 解方程

2019-11-10 16:48:42

字體：大中小

來源：轉載

供稿：網友

解普通方程

from sympy import *x = Symbol('x')y = Symbol('y')PRint solve([2 * x - y - 3, 3 * x + y - 7],[x, y])

test

解微積分

from sympy import *n = Symbol('n')s = ((n+3)/(n+2))**n#無窮為兩個小寫oprint limit(s, x, oo)

求定積分

from sympy import *t = Symbol('t')x = Symbol('x')m = integrate(sin(t)/(pi-t),(t,0,x))n = integrate(m,(x,0,pi))print n

這里寫圖片描述

解微分方程

#y' = 2xy 的通解from sympy import *f = Function('f')x = Symbol('x')print dsolve(diff(f(x),x) - 2*f(x)*x,f(x))#說明：f = Function('f')x = Symbol('x')#表示f(x)的導：diff(f(x), x, index) >>> diff(sin(x), x, 1)cos(x)dsolve(eq, f(x)) #第一個參數為微分方程（要先將等式移項為右端為0的形式) #第二個參數為要解的函數(在微分方程中)

矩陣化簡

這里寫圖片描述

from sympy import *x1,x2,x3 = symbols('x1 x2 x3')a11,a12,a13,a22,a23,a33 = symbols('a11 a12 a13 a22 a23 a33')m = Matrix([[x1,x2,x3]])n = Matrix([[a11,a12,a13],[a12,a22,a23],[a13,a23,a33]])v = Matrix([[x1],[x2],[x3]])f = m * n * vf[0] 化簡， subs代入計算print f[0].subs({x1:1, x2:1, x3:1})

上一篇：Openstack 實現技術分解 (3) 開發工具 — VIM & dotfiles

下一篇：華為機試在線訓練-牛客網（16）字符串排序