program bh
implicit none
real :: w(100)
call random_seed ()     ! 系统根据日期和时间随机地提供种子
call random_number (w)  ! 每次的随机数就都不一样了
print*,w
end
这个程序可以随机产生100个0-1之间的数，如果下面还有程序，要随机读取这100个数值中的某几个数，该怎么办？比如：
program main
  use IMSL
  implicit none
  real::ra 
  complex:: c(4,4) 
  real, parameter::pi=3.14159265
  complex, parameter:: fi=(0.0, 1.0)
  integer ::n0
  real w,t,f

  f=0.125
  n0=1 
  c=0.0     
  w=0.5  
  t=1.0

  call random_seed ()     ! 系统根据日期和时间随机地提供种子
  call random_number (ra)  ! 每次的随机数就都不一样了
  w=w-ra            
      c(1,1)=w
      c(2, 1) = t*exp(fi*(-pi+2.0*pi/3.0*f*(1.5*n0+0.25))) 
      c(1, 2) = t*exp(-fi*(-pi+2.0*pi/3.0*f*(1.5*n0+0.25)))     
      c(2, 2) = w
      c(2, 3) = t
      c(3, 2) = t
      c(3, 3) = w
      c(3, 4) = t*exp(fi*(-pi+2.0*pi/3.0*f*(1.5*(n0+1)+0.25)))
      c(4, 3) = t*exp(-fi*(-pi+2.0*pi/3.0*f*(1.5*(n0+1)+0.25)))
      c(4, 4) = w 
 end program main

这个程序的思路只能使得w取某一个-0.5---0.5之间的值，即c(1,1)=c(2,2)=c(3,3)=c(4,4)=w,但是如果我想让每个w都取一个-0.5---0.5之间的随机值，随机选取，即c(1,1),c(2,2),c(3,3),c(4,4)分别随机选取一个-0.5---0.5之间的随机值，这样办行不？

program main
  use IMSL
  implicit none
  real::ra(100) ！种子，用于产生随机数用
  complex:: c(4,4) 
  real, parameter::pi=3.14159265
  complex, parameter:: fi=(0.0, 1.0)
  integer ::n0 
  real w,t,f

  f=0.125
  n0=1 
  c=0.0     
  w=0.5  
  t=1.0

  call random_seed ()     ! 系统根据日期和时间随机地提供种子
  call random_number (ra)  ! 每次的随机数就都不一样了
  w=w-ra        ！产生100个-0.5---0.5之间的随机数    
      c(1,1)=w
      c(2, 1) = t*exp(fi*(-pi+2.0*pi/3.0*f*(1.5*n0+0.25))) 
      c(1, 2) = t*exp(-fi*(-pi+2.0*pi/3.0*f*(1.5*n0+0.25)))     
      c(2, 2) = w
      c(2, 3) = t
      c(3, 2) = t
      c(3, 3) = w
      c(3, 4) = t*exp(fi*(-pi+2.0*pi/3.0*f*(1.5*(n0+1)+0.25)))
      c(4, 3) = t*exp(-fi*(-pi+2.0*pi/3.0*f*(1.5*(n0+1)+0.25)))
      c(4, 4) = w 
 end program main
这样能让每个w都取到100个随机数中的某个数不？