Curve with ggplot2

    ggplot(data.frame(x = c(0, 1)), aes(x)) +
      stat_function(fun=function(x){
        dgunit(y = x, mu = mu, phi = phi)
      }, color = "red2", size = 1.3) +
      guides(fill = FALSE) +
      theme_minimal() +
      ylab(expression("f(y)")) +
      xlab(expression(y)) +
      theme(
  
        axis.title.x = element_text(size = 16),
        axis.title.y = element_text(size = 16)
      )