y := operator 'y
x := operator 'x
eq1 := D(y t, t) = x(t)
eq2 := D(x t, t) = -y(t)
solve([eq1,eq2], [x,y], t=0,[x(0)=1,y(0)=0])
For this particular problem, you could have used
(1) -> x := operator 'x
(1) x
Type: BasicOperator
(2) -> eq := D(x t,t, 2) = - x t
,,
(2) x (t)= - x(t)
Type: Equation(Expression(Integer))
(3) -> solve(eq,x,t=0, [1,0])
(3) cos(t)
Type: Union(Expression(Integer),...)
Reference:
http://axiom-wiki.newsynthesis.org/uploads/chapter-8.xhtml#subsec-8.10.1
I also must agree with Martin. I haven't easily found a solve command
that deals with initial values for systems of ODEs.
Ralf
--
You received this message because you are subscribed to the Google Groups "FriCAS -
computer algebra system" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/fricas-devel?hl=en.