00001 ! 00002 implicit none 00003 ! 00004 real(8) :: x(4), f(4), y(1:10,1:20) 00005 real(8) :: rrff 00006 integer :: ii, ir0, jj 00007 ! 00008 do ii = 1, 10 00009 x(ii) = dble(ii) 00010 end do 00011 do jj = 1, 10 00012 do ii = 1, 10 00013 y(ii,jj) = 10. + dble(ii) 00014 end do 00015 end do 00016 ! 00017 ir0 = 3 00018 f(1:4) = 2.0*y(ir0:ir0+3,9) 00019 write(6,'(1p,5e11.3)') f(1:4) 00020 end