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