f =( 1/(A*sqrt(2*pi)) ) * exp(-0.5*((x-B)/A).^2 );
plot(f)
hold on
[a,b] = max(f);
plot(b,a ,'*r')