Generate a random variable between 0 and 1 uniformly and use inverse cdf of a gaussian to compute the samples from a normal distribution.

Use Central Limit Theorem. Generate samples of size 'n' with replacement using rand() and compute mean. The distribution of such means will be gaussian irrespective of distribution being followed in rand().

rand(m,n) produces a m*n matrix, where each element is a random observation from N(0,1).