The fitdist (link) 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 (link) or Plot the pdf of a Standard Normal Distribution (link) 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
