function  [out,iseed,iy,iv]=randf_1(nmd,alow,ahigh,idum,iy,iv)

for j = 1:nmd
    [ran1,iseed,iy,iv]=randf_2(idum,iy,iv);
    idum=iseed;
   out(j)=(ran1*(ahigh-alow))+alow;
end