In [ ]:
In [ ]:
Zbadaj ruch punktu materialnego w poniższym potencjale:
In [3]:
U(x) = x^4+1/4*x-x^2
plot(U(x),(x,-1.2,1.2),figsize=5)
Out[3]:
In [ ]:
In [ ]:
Rozwiązać numerycznie równanie Fishera-Kolomogorowa. $$\frac{ \partial u(x,t)}{\partial t} = u (1-u) + \frac{ \partial^2 u(x,t)}{\partial x^2} $$
Model ten wykazuje rozwiązanie będące frontem falowym poruszającym się ze stałą prędkością. Obliczyć z otrzymanego rozwiązania numerycznego prędkość frontu falowego w modelu i porównać z wartością przewidywaną przez teorię.
In [7]:
var('x,t')
b=5/sqrt(6)
cs = 1/(1+ exp((x-b*t)/sqrt(6)))^2 ## rozwiązanie dokładne
show(cs)
bool(cs.diff(t)==cs*(1-cs)+cs.diff(x,2))
Out[7]:
In [ ]: