blob: 1cbd30bd0c363807cfb4ca4f229ee90b422d8ea1 (
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
47
48
49
50
51
52
53
54
55
|
\documentclass{article}
\usepackage{axiom}
\begin{document}
\title{\$SPAD/src/input expexpan.input}
\author{The Axiom Team}
\maketitle
\begin{abstract}
\end{abstract}
\eject
\tableofcontents
\eject
\section{License}
<<license>>=
--Copyright The Numerical Algorithms Group Limited MCMXCII.
@
<<*>>=
<<license>>
--% ExponentialExpansion
-- We compute exponential expansions and right hand limits.
-- Author: Clifton J. Williamson
-- Date Created: 27 August 1992
-- Date Last Updated: 27 August 1992
-- Keywords: Exponential expansion, limit, right-hand limit
-- References:
xxp f == exprToXXP(f,true)$FS2EXPXP(INT,EXPR INT,x,0)
f1 := (a**2 + 1) * exp(1/x**3 + 2/x**2) - exp(b) * exp(1/x**3 + 3/x**2)
x1 := xxp f1
limitPlus x1 -- %minusInfinity
f2 := (a**2 + 1) * exp(1/x**3 + 2/x**2) - exp(b) * exp(-1/x**3 + 3/x**2)
x2 := xxp f2
limitPlus x2 -- %plusInfinity
f3 := (a**2 + 1) * exp(1/x**3) - exp(b) * exp(c/x**2)
x3 := xxp f3
limitPlus x3 -- %plusInfinity
f4 := (a**2 + 1) * exp(-1/x**3) - exp(b) * exp(c/x**2)
x4 := xxp f4
limitPlus x4 -- "failed"
p5 := tan(x) * exp(1/x**2) - tan(x) * exp(1/x**2 - 1/x) + sin(x) * exp(1/x)
q5 := -4 * exp(-1/x**2 - 1/x) + sin(x) * exp(-1/x**2 + 1/x)
f5 := p5 / q5
x5 := xxp f5
limitPlus x5 -- %plusInfinity
@
\eject
\begin{thebibliography}{99}
\bibitem{1} nothing
\end{thebibliography}
\end{document}
|
