n <- 3
k <- 0.2
ks <- 0.1
kl <- 0.3
x <- seq(0.1,3,by=0.05)
y <- 1/(k*gamma(n))*(x/k)^(n-1)*exp(-x/k)
plot(x,y,col=2, xlab= " ", ylab="Q", xlim=range(0:3), ylim=range(0:3))
abline( h = seq( 5, 25, 5 ),lty = 2,col = colors()[440])
abline( v = seq( 1, 4, 1 ),lty = 2,col = colors()[440])
curve(1/(ks*gamma(n))*(x/ks)^(n-1)*exp(-x/ks), add = TRUE, col="blue")
curve(1/(k *gamma(n))*(x/k )^(n-1)*exp(-x/k ), add = TRUE, col="red")
curve(1/(kl*gamma(n))*(x/kl)^(n-1)*exp(-x/kl), add = TRUE, col="green")