65 integer,
intent(in) :: n
66 real(kind=kind(1.d0)),
intent(inout) :: a(n,n)
68 real(kind=kind(1.d0)) :: r
87 integer,
intent(in) :: n
88 real(kind=kind(1.d0)),
intent(inout) :: a(n,n)
100 integer,
intent(in) :: n
101 real(kind=kind(1.d0)),
intent(inout) :: b(n,n)
103 real(kind=kind(1.d0)) :: rj,dj,ri,d
107 dj = (0.9 + (0.001 + 0.05*rj)*rj)*rj
110 b(i,j) = 0.2 + (0.3 + (0.4 + 0.01*ri)*ri)*ri - dj
128 subroutine initx(n,x)
130 integer,
intent(in) :: n
131 real(kind=kind(1.d0)),
intent(out) :: x(n)
133 real(kind=kind(1.d0)) :: d
146 integer,
intent(in) :: n
147 real(kind=kind(1.d0)),
intent(in) :: a(n)
148 real(kind=kind(1.d0)),
intent(out) ::b
153 end module initialise
subroutine to_test(n, a, b)
subroutine initx(n, x)
X initialization.