rng(1);% For å kunne reprodusere eksempelet
mu=180;% ukjend
sigma= 10;
 
n=25;
% Lengde konfidensinterval
delKI=1.96*sigma/sqrt(n);
 
% Gjenta B ganger
B=100;
 
% konfidensintervall
figure(1);
clf;
plot(mu*ones(1,B),'r');
hold;
nOutside = 0;
pause
for b=1:B
  x=mu+sigma*randn(n,1);
  xbar=mean(x);
 
  % Øvre og nedre grense
  mu_U=xbar+delKI;
  mu_L=xbar-delKI;
 
  plot(b,mu_L,'.b','MarkerSize',16);
  plot(b,mu_U,'.b','MarkerSize',16);
  plot([b b],[mu_L mu_U],'k');
 
  if (mu_U < mu || mu_L > mu)
    nOutside = nOutside + 1;
  end
 
  if (b < 10)
    pause
  else
    pause(0.2)
  end
end
xlabel('Replikasjoner');
ylabel('Konfidensintervall for mu');
 
disp(['Antall konfidensintervall som bommer på mu:',num2str(nOutside)])