Problema dos Desencontros
11/10/2007
Cálculo do número D_n de desencontros de ordem n
| > | Des := proc(n::integer) |
| > | local k; |
| > | RETURN(sum((-1)^k * (n!/k!), k=0..n)); |
| > | end; |
| > | Des(5); |
| > | seq(Des(i),i=2..15); |
Cálculo da probabilidade
| > | Digits := 20; |
| > | for i from 2 to 20 do |
| > | i, evalf(Des(i)/i!); |
| > | od; |
Gráfico da probabilidade
| > | with(plots); |
| > | grafprob := proc(n::integer) |
| > | local w,i,j; |
| > | for i from 2 to n do |
| > | w[i]:=Des(i)/i!; |
| > | od; |
| > | pointplot({seq([j,w[j]],j=2..n)},symbol=circle,symbolsize=8,color=red,axes=boxed,title=`Probabilidade para n em [2,20]`,style=LINE): |
| > | end: |
| > |
Warning, the name changecoords has been redefined
| > | grafprob(20); |
![[Maple Plot]](desencontros/desencontros28.gif)
| > |