1subroutine pr40421 (j, q, r)
2 double precision :: q(1,1), r(1,1,3)
3 save
4 integer :: i, j, m, n
5 double precision :: s, t, u
6 do i=1,2
7 do m=1,j
8 do n=1,1
9 s=q(n,m)*r(n,m,1)
10 t=q(n,m)*r(n,m,2)
11 u=q(n,m)*r(n,m,3)
12 end do
13 end do
14 end do
15end
16