Python-koodia, esim. iPythonille
Laminaarinen nopeus
1 2 3 4 5 6 7 8 9 10 | import sympy as symfrom sympy import pprintv = sym.Symbol('v')F,alpha = sym.symbols('F alpha')a,m,t = sym.symbols('a m t')eq1L = m*sym.diff( v(t), t)eq1R = F - alpha*v(t)sol1 = sym.dsolve( eq1L - eq1R, v(t) )pprint( sol1 ) |
Antaa vastauksen
$v(t) = \frac F\alpha + \frac1\alpha e^{\frac{\alpha(C_1 -t)}m}$
Integroimisvakio $C_1$ ratkaistaan alkuehdosta $v(0)=0$, josta saadaan
1 2 | C1 = sym.symbols('C1')sym.solve( sol1.subs(t,0), C1) |
$C1 = m \ln ( -F /\alpha)$
kuten tekstissä.
