% generate a 100x1 matrix filled with normal distributed values
X = randn(100, 1);

% plot the histogram of generated values
subplot(2, 1, 1); hist(X);

% plot sorted X to generate an approximation of the
% integral function of the gaussian distribution
subplot(2, 1, 2); plot(sort(X), 0:0.01:0.99);