blob: e5e4ab8c6aa317b40543aed213e8d46ef4da8f24 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
\documentclass{article}
\usepackage{axiom}
\begin{document}
\title{\$SPAD/src/input robidoux.input}
\author{Nicolas Robidoux}
\maketitle
\begin{abstract}
\end{abstract}
\eject
\tableofcontents
\eject
\section{License}
<<license>>=
--Copyright The Numerical Algorithms Group Limited 1994.
@
<<*>>=
<<license>>
)clear all
-- from the paper:
-- Does OpenAxiom Solve Systems of O.D.E.'s Like Mathematica?
-- Nicolas Robidoux University of New Mexico 24 June 1993
X1:=operator 'X1
deq1:=D(X1 t,t)=(1+ cos t /(2+sin t)) * X1 t
solve(deq1,X1,t)
C1*%.basis.1
function(%,'x1,'t)
x1
X2:=operator 'X2
deq2:=D(X2 t,t)=x1 t - X2 t
solve(deq2,X2,t)
%.particular
simplify %
%+C2*%%(-3).basis.1
function(%,'x2,'t)
x2
-- we can stop here because the answer is improved over the
-- one given in the paper.
x1 t
@
\eject
\begin{thebibliography}{99}
\bibitem{1} nothing
\end{thebibliography}
\end{document}
|
