% To create a 1D Gaussian Distribution and show its histogram
figure(1); hold on;
% plot pdf of the Gaussisan Distribution
X_theo=[-4:0.2:5];
Y_theo=pdfgauss(X_theo,model);
Y_theo=Y_theo/sum(Y_theo);
plot(X_theo,Y_theo,'r');

[Y,X] = hist(gsamp(model,500),50); % compute histogram
bar(X,Y/500); % plot histogram

% To create a 2D Gaussian Distribution and show its histogram
model.Mean = [1;1]; % Mean vector
model.Cov = [1 0.6; 0.6 1]; % Covariance matrix
figure; hold on;
data=gsamp(model,1000);
ppatterns(data); % plot sampled data
pgauss(model); % plot shape