5 subroutine init (v1
, v2
, N
)
14 subroutine check (p
, N
)
16 real, parameter :: EPS
= 0.00001
19 diff
= p(i
) - (i
+ 2.0) * (i
- 3.0)
20 if (diff
> EPS
.or
. -diff
> EPS
) STOP 1
24 subroutine vec_mult (N
)
26 real :: p(N
), v1(N
), v2(N
)
40 use e_50_1_mod
, only
: vec_mult