352 - Topic 17 - Solving 1st order systems -- Sage

# adapted directly from the sage documentation t = var('t') x = function('x')(t) y = function('y')(t) de1 = diff(x,t) - y == 0 de2 = diff(y,t) + x == 0 soln = desolve_system([de1, de2], [x,y]) for sol in soln: sol.show()

$\newcommand{\Bold}[1]{\mathbf{#1}}x\left(t\right) = \cos\left(t\right) x\left(0\right) + \sin\left(t\right) y\left(0\right)$

$\newcommand{\Bold}[1]{\mathbf{#1}}y\left(t\right) = -\sin\left(t\right) x\left(0\right) + \cos\left(t\right) y\left(0\right)$

# now how about a linear predator-prey model with (x_0,y_0) given t = var('t') x = function('x')(t) y = function('y')(t) de1 = diff(x,t) + y - x == 0 de2 = diff(y,t) - x + y == 0 soln = desolve_system([de1, de2], [x,y], ics=[0,2,3]) for sol in soln: sol.show()

$\newcommand{\Bold}[1]{\mathbf{#1}}x\left(t\right) = -t + 2$

$\newcommand{\Bold}[1]{\mathbf{#1}}y\left(t\right) = -t + 3$

# now how about a non-linear predator-prey model with (x_0,y_0) given t = var('t') x = function('x')(t) y = function('y')(t) de1 = diff(x,t) + y*x - x == 0 de2 = diff(y,t) - x*y + y == 0 soln = desolve_system([de1, de2], [x,y], ics=[0,2,3],ivar=t) for sol in soln: sol.show()

$\newcommand{\Bold}[1]{\mathbf{#1}}x\left(t\right) = \mathcal{L}^{-1}\left(-\frac{\mathcal{L}\left(x\left(t\right) y\left(t\right), t, g_{2821}\right) - 2}{g_{2821} - 1}, g_{2821}, t\right)$

$\newcommand{\Bold}[1]{\mathbf{#1}}y\left(t\right) = \mathcal{L}^{-1}\left(\frac{\mathcal{L}\left(x\left(t\right) y\left(t\right), t, g_{2821}\right) + 3}{g_{2821} + 1}, g_{2821}, t\right)$

# now how about a non-linear predator-prey model # with (x_0,y_0) given. # Notice, you can rotate the graph to see # the phase plane x vs y as well as t vs x and t vs y reset() t,x,y = var('t,x,y') de1 = -y*x + 3*x de2 = x*y - y soln = desolve_system_rk4([de1, de2], [x,y], ics=[0,0.6,0.7],ivar=t,end_points=20) points(soln,size=20)

# Here is a problem that is partially decoupled and # which you can solve by hand analytically. But, Sage # (through the package Maxima) can't do it. var('t') x = function('x',t) y = function('y',t) de1 = diff(y,t)+3*y == 0 de2 = diff(x,t)-x/y == 0 soln =desolve_system([de1,de2],[y,x],ivar=t)

Traceback (click to the left of this block for traceback)
...
TypeError: unable to make sense of Maxima expression
'x(t)=ilt((laplace(x(t)/y(t),t,?g3057)+2)/?g3057,?g3057,t)' in Sage

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "_sage_input_17.py", line 10, in <module>
    exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("dmFyKCd0JykKeCA9IGZ1bmN0aW9uKCd4Jyx0KQp5ID0gZnVuY3Rpb24oJ3knLHQpCmRlMSA9IGRpZmYoeSx0KSszKnkgPT0gMApkZTIgPSBkaWZmKHgsdCkteC95ID09IDAKc29sbiA9ZGVzb2x2ZV9zeXN0ZW0oW2RlMSxkZTJdLFt5LHhdLGl2YXI9dCk="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))
  File "", line 1, in <module>
    
  File "/tmp/tmpqRgnEa/___code___.py", line 8, in <module>
    exec compile(u'soln =desolve_system([de1,de2],[y,x],ivar=t)
  File "", line 1, in <module>
    
  File "/home/sageserver/sage/local/lib/python2.6/site-packages/sage/calculus/desolvers.py", line 743, in desolve_system
    soln[i] = sol.sage()
  File "/home/sageserver/sage/local/lib/python2.6/site-packages/sage/interfaces/expect.py", line 1788, in sage
    return self._sage_()
  File "/home/sageserver/sage/local/lib/python2.6/site-packages/sage/interfaces/maxima.py", line 1824, in _sage_
    maxima=self.parent())
  File "/home/sageserver/sage/local/lib/python2.6/site-packages/sage/calculus/calculus.py", line 1654, in symbolic_expression_from_maxima_string
    raise TypeError, "unable to make sense of Maxima expression '%s' in Sage"%s
TypeError: unable to make sense of Maxima expression 'x(t)=ilt((laplace(x(t)/y(t),t,?g3057)+2)/?g3057,?g3057,t)' in Sage

352 - Topic 17 - Solving 1st order systems

23 days ago by Professor352