As per my understanding you want to plot the probability density function with “mu” which must be average of your data. You can do this with the help of “mean” function. You can create gaussian pdf for your data and see how they overlap by setting the distribution name as “Normal”. Here is a sample code:
x1 = [-2 -1 0 1 2];
mu=mean(x1);
pd = makedist('Normal','mu',mu,'sigma',sigma);
y = pdf(pd,x1);
plot(x1,y);
hold on
x2 =(-9:5);
mu=mean(x2);
pd = makedist('Normal','mu',mu,'sigma',sigma);
y = pdf(pd,x2);
plot(x2,y);
