def randf_1(nmd,alow,ahigh,idum,iy,iv):
   from randf_2 import randf_2
   import numpy as np

   out=np.zeros(nmd+1)
   for j in range(1,nmd+1):

      ran1,iseed,iy,iv =randf_2(idum,iy,iv)
      idum=iseed
      out[j]=(ran1*(ahigh-alow))+alow

   return  out,iseed,iy,iv