How to rewrite normrnd?

Hello! Matlab beginner here. I'm trying to adapt a code sent by a colleague which uses the 'normrnd' function, but I don't have the Statistics Toolbox (we're running an old version of Matlab for equipment compatibility reasons). How would I rewrite this function for myself?

I have:

function nrmmatrix = nrmrnd(mu, sigma, sz1, sz2)

nrmmatrix = [];

for i = 1:sz1

for j = 1:sz2

nrmmatrix(i,j) = mu + sigma * rnd(n);

end

end

...But I have a feeling this isn't right!

Thanks very much in advance.

Answer 1

[Jeff Miller]

function nrmmatrix = nrmrnd(mu, sigma, sz1, sz2)
nrmmatrix = mu+sigma*randn(sz1,sz2);
end