thetay1=alpha.*rand(a,b);
     thetay2=alpha.*rand(a,b);
     h_low=min(abs(thetay1),abs(thetay2));
     h_high=max(h_low);
            for ii=1:b
                for jj=1:a
                    if h_high(jj)==h_low(ii,jj)
                       y1(jj)=thetay1(ii,jj);
                       y2(jj)=thetay2(ii,jj);
                    end
                 end
            end