Berld[x_, a_] := N[a^x Exp[-a]/Gamma[x+1,a]] /; (x >= 0) && (a > 0)