function pravd = zajimava_uloha(n)
% pravd = zajimava_uloha(n)
%
% vypocte pravdepodobnost toho, ze kdyz nahodne zamichate prirozena
% cisla 1 az n, nezustane ani jedno cislo na sve puvodni pozici

if n <= 0 | n - floor(n) > eps
    error('n neni prirozene cislo')
end

perm = permutace(n); % generuje vsechna mozna usporadani n cisel
pocet_perm=size(perm,1); % pocet vsech moznych usporadani n cisel

for i=1:pocet_perm
    pom(i,:)=[1:n];
end

 P=(pom==perm); 
 pom2=sum(P,2);
 pom2=~pom2;
 suma=sum(pom2); % pocet vsech usporadani, u kterych neni ani jedno cislo na svem puvodnim miste
 
 pravd=suma/pocet_perm;
 
 disp(['Pravdepodobnost: ', num2str(pravd)])

% konec

