#==============================================================================
# Distribuicao a priori gama com parametros a=4 e b=20000
# Funcao de densidade de probabilidade
#==============================================================================
a    <- 4      # parametro de forma
b    <- 20000  # parametro de escala
xseq <- seq(0.000001,0.001,length=1000) 
yseq <- dgamma(xseq,a,b)

par(mfrow=c(1,1),lwd=2.0,cex.lab=1.5,cex.axis=1.5,lab=c(10,5,5),
    mar=c(5,5,2,2.5),xpd=T,cex.main=2.0)
plot(xseq,yseq,type="l",lwd=2,xlab=expression(lambda),
     ylab=expression(f(lambda)),col="red")