(*podla staqtistik  74percent domacich firiem vlastni chladnicku ,ktora sluyi pre yamestnancov nahodne vzberieme 30 firiem a mam zistii
a - prave 22 vlastni chladnicku
b - najviac 19 vlastni chladnicku 
c aspon 15 vlastni chladikuc

PDF - nahodne rozdelenie funkcie
PDF[rozdelenie,x] = P(X=x)   za x dame hladanu
probability density function

CDF[rozdelenie,x] = P(x≤x)   cumulative

nacitanie binomickeho rozdelenia 
*)

bi=BinomialDistribution[n,p]

(*Mean[bi]  zrata sa stredna hodnota*)

(*Variance[bi]    = np(1-p)   =npq  *)
Mean[bi]

n=30;

p=0.74;

PDF[bi,22]

(*najviac 19      p0+p1+.....+p19 = CDF[bi,19]*)

bi=BinomialDistribution[n,p]

CDF[bi,19]

(*aspon15 vlastni*)
(*p15 + ...+p30 *)
(* 1-(p0+....+p14)  *)
(*inak ako takzmto sposobom to nejde*)

1-CDF[bi,14]

0.9987051566809445`

(*aspon 3 firmy ale najviac 11 firiem     p03+..........+p11*)

(*------------vysledok ma byt 0.00014------------*)

n=30
p=0.74
(∑_(i=3)^11)PDF[bi,i]

(*nakreslite graf rozdelenia pravdepodobnosti   treba nacitat baliky <<"BarCharts = `";   *)

<<"BarCharts`";

(*demonstracia jej sposobu mz mozeme mat vlastnz dame si 31prvkove pole..nachystame..tam budu hodnoty p0 az p30*)
tab= Array[w,{31}];

(*vo w[1]=p0 preto to treba s tymi posunmi lebo vo w1 je p0 atd*)
Do[w[i+1]=PDF[bi,i],{i,0,30}]
tab

tab

(∑_(i=1)^30)w[i+1]   (*zosunujeme hodnoty*)

BarChart[tab]

(*----------------------samostatna praca -----------------*)



(*v krajine X ma 90percent domacnosti pripojenie na internet   nahodne vybereme 20 domacnosti ..aka je pravdepodobnost 
a---ze prave 15 ma net
  b-- najviac 16 ma net
c--- aspon 18 ma net
d-- graficky rozdelenie pravdepodobnost*)


n=20;

p=0.9;
bi=BinomialDistribution[n,p]

PDF[bi,15]

CDF[bi,16]

1-CDF[bi,17]

tab= Array[w,{21}];

Do[w[i+1]=PDF[bi,i],{i,0,20}]

tab

BarChart[tab]

clear n

(*hypergonomicke rozdelenie *)
hy=HypergeometricDistribution[n,M,N]

(*   (3)n  - pocet vybranych - kolko vyberame   
	M- pocet s danou vlastnostou (v celom vybere  -- napr teray 5zlych)
		N- celkovy rozsah   --teraz 15

 *)

vyrobok=HypergeometricDistribution[3,5,15]

(*a - vsetky 3 su dobreeee teda ani jeden zly *)
PDF[vyrobok,0]

(*pocet desatinnych miest cez percento sa odvolame na posledny vystup*)
%//N

(*najviac 2 dobre ale nas zaunjimaju zle----------takze aspon 1 zly  teda p1+p2+p3----este chyba p0...takze 1-p0 to jest 1-PDF[....*)
1-PDF[vyrobok,0]//N

(*najviac 2 zle vyrobky *)
CDF[vyrobok,2]//N

(*firma vyrabas 500ks kosiel rocne..v dosledku starych strojov je 70 z nich roznej akosti...nahodne vyberiem 7 na kontrolu..aka je pravdepodobnost,ze....
a- dve budu inej akosti
b-najviac jedna bude inej akosti*)

vyrobok=HypergeometricDistribution[7,70,500]

PDF[vyrobok,2]//N