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