function [A] = uniform(m,n,a,b)
%  generate a uniformly distributed mxn matrix with entries in [a,b]

A = a*ones(m,n) + (b-a)*rand(m,n);

return;
