%KOVANEC ocenjevanje nakljucnosti meta kovanca
%r=KOVANEC(n,q);
%r=1 ce zavrnemo hipotezo, n stevilo poskusov, q verjetnost pravilnosti
%zavrnitve
function r=kovanec(n,q);
  a=rand(1,n); % n nakljucnih metov
  N=[sum(a<=0.40),sum(a>0.40)]
  pricak=n/2;
  c=sum((N-pricak).^2/pricak)%vrednost X^2
  P=gammainc(c/2,1/2) %verjetnost X^2<=c
  if P>=q
      r=1;
      sprintf('Hipotezo zavrnemo z verjetnostjo %2.2f!',...
      q)
  else 
      r=0;
      sprintf('Hipoteze ne moremo zavrniti z verjetnostjo %2.2f!',...
      q)
  end
  