figure
m=-1;  sd=1;
x=-5:.05:5;
z=x(x<0);
y=normpdf(x,m,1);
plot(x,y,'k')
hold on
plot(x,min(y,normpdf(-x,m,sd)),'k')
%fill(x,min(y,normpdf(-x,m,sd)),'r')
mu=cmean(z,normpdf(z,m,sd)-normpdf(-z,m,sd));
    m=m-.0125;
    sd=sd*1.01;
    y=normpdf(x,m,sd);
hold on
plot(x,min(y,normpdf(-x,m,sd)),'k')
%fill(x,min(y,normpdf(-x,m,sd)),'y')