# Section 6; Question 3:
#
lambda = 1.5/0.75
r = lambda * 1.5
p = pgamma( 2.5, shape=r, scale=1/lambda ) - pgamma( 1.0, shape=r, scale=1/lambda )
n = 100
print( c( p, n*p ) )