dgamma(50, 5, 5)

## [1] 2.172e-100

x = 0:10
curve(dgamma(x, 2, 1), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 3, lwd = 3, 
    main = "Gráfico distribución Gamma")
curve(dgamma(x, 1, 2), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 5, lwd = 3, 
    add = TRUE)
curve(dgamma(x, 11, 3), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 7, lwd = 3, 
    add = TRUE)
legend("topright", c("a=2,ß=1", "a=1,ß=2", "a=11,ß=3"), col = c(3, 5, 7), lwd = 3, 
    inset = 0.05)

curve(dgamma(x, 2, 1), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 3, lwd = 3, 
    main = "Gráfico distribución Gamma (a=2,ß=1)")

curve(dgamma(x, 1, 2), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 5, lwd = 3, 
    main = "Gráfico distribución Gamma (a=1,ß=2)")

curve(dgamma(x, 11, 3), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 7, lwd = 3, 
    main = "Gráfico distribución Gamma (a=11,ß=3)")

curve(pgamma(x, 2, 1), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 3, lwd = 3, 
    main = "Gráfico de la función de distribución Gamma")
curve(pgamma(x, 1, 2), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 5, lwd = 3, 
    add = TRUE)
curve(pgamma(x, 11, 3), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 7, lwd = 3, 
    add = TRUE)
legend("bottomright", c("a=2,ß=1", "a=1,ß=2", "a=11,ß=3"), col = c(3, 5, 7), 
    lwd = 3, inset = 0.05)

curve(pgamma(x, 2, 1), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 3, lwd = 3, 
    main = "Gráfico de la función de distribución Gamma (a=2,ß=1)")

curve(pgamma(x, 1, 2), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 5, lwd = 3, 
    main = "Gráfico de la función de distribución Gamma (a=1,ß=2)")

curve(pgamma(x, 11, 3), xlab = "x", ylab = "f(x;a,ß)", 0, 10, col = 7, lwd = 3, 
    main = "Gráfico de la función de distribución Gamma (a=11,ß=3)")