Evaluate a distribution programmatically

The fitdist  function will fit a distribution to data, producing a probability distribution object:

pd = fitdist(x,distname)

You can then calculate the probabilities with it using the Plot Standard Normal Distribution cdf or Pl

ot the pdf of a Standard Normal Distribution functions:

data = 5 + 2*randn(1, 100);                         % Create Data
pd = fitdist(data(:), 'Normal');

x = linspace(-2, 12);
Npdf = pdf(pd, x);
Ncdf = cdf(pd, x);

figure
plot(x, Npdf,   x, Ncdf)
grid