t = var('t')
x = function('x')(t)
y = function('y')(t)
de1 = diff(x,t) == -3/20*x + 1/12*y
de2 = diff(y,t) == 3/20*x-3/20*y
sol = desolve_system([de1,de2],[x,y],ics=[0,4,0])

f(t) = sol[0].rhs() g(t) = sol[1].rhs()

1 | initial version |

t = var('t')
x = function('x')(t)
y = function('y')(t)
de1 = diff(x,t) == -3/20*x + 1/12*y
de2 = diff(y,t) == 3/20*x-3/20*y
sol = desolve_system([de1,de2],[x,y],ics=[0,4,0])

f(t) = sol[0].rhs() g(t) = sol[1].rhs()

Listed below is my code:

```
t = var('t')
x = function('x')(t)
y = function('y')(t)
de1 = diff(x,t) ==
```~~-3/20~~*x **-3/20*x + *~~1/12~~y
1/12*y
de2 = diff(y,t) == ~~3/20~~*x-3/20*y
3/20*x-3/20*y
sol = ~~desolve_system([de1,de2],[x,y],ics=[0,4,0])~~desolve_system([de1,de2],[x,y],ics=[0,4,0])
~~ ~~f(t) = sol[0].rhs()
g(t) = ~~sol[1].rhs()~~

sol[1].rhs()

Listed below is my code:

```
t = var('t')
x = function('x')(t)
y = function('y')(t)
de1 = diff(x,t) == -3/20*x + 1/12*y
de2 = diff(y,t) == 3/20*x-3/20*y
sol = desolve_system([de1,de2],[x,y],ics=[0,4,0])
f(t) = sol[0].rhs()
g(t) = sol[1].rhs()
```

I need to plot two differential equations on the same set of ~~axes ~~axes
(with a legend showing which is ~~which). ~~which).

My domain has to be ~~0≤t≤100. ~~0 ≤ t ≤ 100.

Listed below is my code, any help is greatly appreciated.~~ ~~

Listed below is my code:

```
t = var('t')
x = function('x')(t)
y = function('y')(t)
de1 = diff(x,t) == -3/20*x + 1/12*y
de2 = diff(y,t) == 3/20*x-3/20*y
sol = desolve_system([de1,de2],[x,y],ics=[0,4,0])
f(t) = sol[0].rhs()
g(t) = sol[1].rhs()
```

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.