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
56
57
|
\documentclass{article}
\usepackage{axiom}
\begin{document}
\title{\$SPAD/src/input f04mbf.input}
\author{The Axiom Team}
\maketitle
\begin{abstract}
\end{abstract}
\eject
\tableofcontents
\eject
\section{License}
<<license>>=
--Copyright The Numerical Algorithms Group Limited 1994.
@
<<*>>=
<<license>>
)cl all
showArrayValues true
showScalarValues true
n := 10
b := new(1,n,0.0)$Matrix DoubleFloat;
b(1,1) := 6
for i in 2..n-1 repeat
b(1,i) := 4
b(1,n) := 6
precon := true
shft := 0
-- The following two parameters are not used
-- by the ASPs. They may be useful if a user
-- supplies FORTRAN code directly.
lrwork:=1
liwork:=1
rtol := 0.00001
itnlim := 100
msglvl := 1
ifail := 1
z:=new(10,10,0.0)$Matrix MFLOAT;
z(1,1):=2;z(1,2):=1;z(1,10):=3
z(10,1):=3;z(10,9):=1;z(10,10):=2
for i in 2..(n-1) repeat
z(i,i-1):=1
z(i,i):=2
z(i,i+1):=1
z
m :Matrix MFLOAT:= [[2,1,0,0,0,0,0,0,0,0],[1,2,1,0,0,0,0,0,0,0],_
[0,1,2,1,0,0,0,0,0,0],[0,0,1,2,1,0,0,0,0,0],[0,0,0,1,2,1,0,0,0,0],_
[0,0,0,0,1,2,1,0,0,0],[0,0,0,0,0,1,2,1,0,0],[0,0,0,0,0,0,1,2,1,0],_
[0,0,0,0,0,0,0,1,2,1],[0,0,0,0,0,0,0,0,1,2]]
result:=f04mbf(n,b,precon,shft,itnlim,msglvl,lrwork,liwork,rtol,ifail,z::ASP28(APROD),(m::ASP34(MSOLVE)))
@
\eject
\begin{thebibliography}{99}
\bibitem{1} nothing
\end{thebibliography}
\end{document}
|
