13 parameter(half=0.5d0,one=1.d0,two=2.d0,three=3.d0)
17 Subroutine normmax(vektor, n, normaxi)
18 REAL*8 normaxi, vektor(0:n-1)
30 betrag=dabs(vektor(i))
31 if (betrag > norm) norm = betrag
35 end subroutine normmax
37 Subroutine copymat(ziel, quelle, n, m)
38 REAL*8 ziel(0:n-1,0:m-1),quelle(0:n-1,0:m-1)
48 ziel(i,j) = quelle(i,j)
52 end subroutine copymat
57 Subroutine readvec (u, n, x)
60 read(u,*) (x(i),i=0,n-1)
62 end subroutine readvec
67 Subroutine readvec1 (u, n, x)
70 read(u,*) (x(i),i=1,n)
72 end subroutine readvec1
78 Subroutine setvec(n,x,val)
91 Subroutine setmat(n,m,a,val)
100 end subroutine setmat
102 Subroutine writevec(u, n, x)
121 write(u,10) (x(i),i=0,n-1)
124 end subroutine writevec
126 Subroutine writevec1 (u, n, x)
145 write(u,10) (x(i),i=1,n)
148 end subroutine writevec1
150 Subroutine readmat (u,n, m, a)
174 REAL*8 a(0:n-1,0:m-1)
175 Read(u,*) ((a(i,j),j=0,m-1),i=0,n-1)
177 end subroutine readmat
179 Subroutine readmat1(u,n, m, a)
204 Read(u,*) ((a(i,j),j=1,m),i=1,n)
206 end subroutine readmat1
208 Subroutine writemat(u,n, m, a)
226 REAL*8 a(0:n-1,0:m-1)
228 write(u,10) (a(i,j),j=0,m-1)
232 end subroutine writemat
234 Subroutine writemat1(u,n, m, a)
255 write(u,10) (a(i,j),j=1,m)
259 end subroutine writemat1
261 Subroutine writehead (u, nom)
281 Write(u,*)
'----------------------------------------------------------' 283 Write(u,*)
'----------------------------------------------------------' 285 end subroutine writehead
287 Subroutine writeend(u)
301 Write(u,*)
'----------------------------------------------------------' 303 end subroutine writeend