>>714
意外にも答が二個あった。

> (lambda1=uniroot(function(lambda) dpois(mu,lambda) - dbinom(mu,n,p) ,c(0,9))$root)
[1] 8.661711
> (lambda2=uniroot(function(lambda) dpois(mu,lambda) - dbinom(mu,n,p) ,c(9,10))$root)
[1] 9.34699
> dpois(mu,lambda1)
[1] 0.1308992
> dpois(mu,lambda2)
[1] 0.1308991
> dbinom(mu,n,p)
[1] 0.1308991