# plot the gamma p.d.f 

x = seq(0,10,0.1)

alpha = 10
lambda = 5

curve(dgamma(x, shape = alpha, rate = lambda), 0, 20)


alpha = 20	  # bigger alpha here	
lambda = 5        # bigger lambda here

curve(dgamma(x, shape = alpha, rate = lambda), 0, 20, add = TRUE, lty = 1)

for(i in 1:10){

alpha = 20 + 10*i	  # bigger alpha here	
lambda = 5                # bigger lambda here

curve(dgamma(x, shape = alpha, rate = lambda), 0, 20, add = TRUE, lty = i)

}