Lambda=seq(.5,20,by=.5)
#logL=-Lambda+10*log(Lambda)-lgamma(11)
logL=-Lambda+9*log(Lambda)-lgamma(10)
par(mfcol=c(1,2))
plot(Lambda,logL,type="l",ylab="Log Likelihood")
plot(Lambda,exp(logL),type="l",ylab="Likelihood")