Python-koodia, esim. iPythonille

Laminaarinen nopeus

import sympy as sym
from sympy import pprint

v = 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

C1 = sym.symbols('C1')
sym.solve( sol1.subs(t,0), C1)

$C1 = m \ln ( -F /\alpha)$

kuten tekstissä.