1 urna 
>> a=25;			%biele
>> b=75;			%cierne
>> pb=25/100;			%pravedepodobnost ze vytiahnem bielu
>> n=99;  			%pocet opakovany
>> for k=0:n			%pripustne k lebo taham n krat
>> x(k+1) = factorial(n)/factorial(k)/factorial(n-k) * pb^k * (1-pb)^(n-k);       %pravdepodobnost ze z 99 opakovany bude k bielych a n-k ciernych 
>> end                          
>> pr=sum(x(1:25))              %pravdepodobnost ze vytiahneme aspon 75 ciernych
>> y=round(rand*n)              %nahodna premenna
>>
>>

2. meskanie autobusu
>> lambda=0.04;				%zadefinovanie lambdy
>> f=inline('lambda.*exp(-lambda.*x)'); %funkcia hustoty
>> s=sym('s');				%symbolicka premenna s
>> F=int(f(s));				%distribucna funkcia
>> minut30aViac=exp(-lambda*30)    	%pravdepodobnost ze bude cakat hodinu a viac lebo f(30<x<oo)=F(oo)-F(30)=1-F(30)  a F(x)=1-exp(-lambda*x)
>> od0po20=1-exp(-lambda*20)       	%pravdepodobnost ze autobus bude meskat 0-20 minut lebo f(0<x<20)=F(20)-F(0)= F(20)