function values=normrand(mu,sigma,n,m);
% Usage:   V = normrand( mu, sigma)        % for 1 x 1
%    or    V = normrand( mu, sigma, n)     % for n x n
%    or    V = normrand( mu, sigma, n, m)  % for n x m
%
% computes a matrix of size n x m, normally distributed 
%             with mean mu and standard deviation sigma.
if nargin < 2, help normrand, return, end
if nargin == 2, 
    n=1; m=1;
elseif nargin == 3, 
    m=n;
end
 z=randn(n,m);
 z=mu+z*sigma;
 values=z; 
