mu=10;
N=25;
al=1;
pauze=1;


x=zeros(1,N);
Emu_new=0;

figure(1);
plotGauss1D(0,1);
axis([-3 15 0 0.5]);
pause;

for j=1:2
for i=1:N
   
   x(i)=randn+mu;
   Emu_old=Emu_new;
   Emu_new=sum(x)/(i+(1/al^2));
   figure(1);
   plotGauss1D(Emu_old,1);hold on;
   plotGauss1D(Emu_new,1);
   plot(x(1:i),zeros(1,i),'ro');hold off;
   axis([-3 15 0 0.5]);
   drawnow;
   if pauze == 1
      pause;
   end  
   
end
N=50;
pauze=0;
x=zeros(1,N);
Emu_new=0;

end