curve(dnorm(x, mean = 0, sd = 1), from = -4, to = 4,
      main = "Kurva Distribusi Normal",
      xlab = "X", ylab = "Kepadatan (Density)",
      col = "blue", lwd = 2)

# Parameter untuk 3 distribusi Binomial yang berbeda
n1 = 10; p1 = 0.5   # mu = n*p, varians = n*p*(1-p)
n2 = 20; p2 = 0.3
n3 = 30; p3 = 0.2

x1 = 0:n1
x2 = 0:n2
x3 = 0:n3

y1 = dbinom(x1, n1, p1)   # dbinom = fungsi pmf binomial
y2 = dbinom(x2, n2, p2)
y3 = dbinom(x3, n3, p3)

mu1 = n1*p1;      var1 = n1*p1*(1-p1)
mu2 = n2*p2;      var2 = n2*p2*(1-p2)
mu3 = n3*p3;      var3 = n3*p3*(1-p3)

# Plot distribusi pertama
plot(x1, y1, type = "h", lwd = 2, col = "blue",
     xlim = c(0, 30), ylim = c(0, 0.35),
     main = "Distribusi Binomial dengan Mu dan Varians Berbeda",
     xlab = "X", ylab = "P(X = x)")
points(x1, y1, pch = 16, col = "blue")

# Tambahkan distribusi kedua dan ketiga di plot yang sama
lines(x2, y2, type = "h", lwd = 2, col = "red")
points(x2, y2, pch = 16, col = "red")

lines(x3, y3, type = "h", lwd = 2, col = "darkgreen")
points(x3, y3, pch = 16, col = "darkgreen")

legend("topright",
       legend = c(paste0("n=", n1, ", p=", p1, " (μ=", mu1, ", σ²=", round(var1,2), ")"),
                  paste0("n=", n2, ", p=", p2, " (μ=", mu2, ", σ²=", round(var2,2), ")"),
                  paste0("n=", n3, ", p=", p3, " (μ=", mu3, ", σ²=", round(var3,2), ")")),
       col = c("blue", "red", "darkgreen"), lwd = 2)

# Parameter untuk 3 distribusi normal yang berbeda
mu1 = 0;  sigma1 = 1     # varians = sigma1^2 = 1
mu2 = 0;  sigma2 = 2     # varians = sigma2^2 = 4
mu3 = 3;  sigma3 = 1.5   # varians = sigma3^2 = 2,25

x = seq(-10, 10, length = 500)

y1 = dnorm(x, mean = mu1, sd = sigma1)
y2 = dnorm(x, mean = mu2, sd = sigma2)
y3 = dnorm(x, mean = mu3, sd = sigma3)

# Plot kurva pertama
plot(x, y1, type = "l", lwd = 2, col = "blue",
     ylim = c(0, 0.45),
     main = "Distribusi Normal dengan Mu dan Varians Berbeda",
     xlab = "X", ylab = "Kepadatan (Density)")

# Tambahkan kurva kedua dan ketiga di plot yang sama
lines(x, y2, lwd = 2, col = "red")
lines(x, y3, lwd = 2, col = "darkgreen")

legend("topright",
       legend = c(paste0("μ=", mu1, ", σ²=", sigma1^2),
                  paste0("μ=", mu2, ", σ²=", sigma2^2),
                  paste0("μ=", mu3, ", σ²=", sigma3^2)),
       col = c("blue", "red", "darkgreen"), lwd = 2)

# Data sesuai tabel
kurva = c(1, 2, 3, 4, 5)
rata_rata = c(0, 2, 5, 10, 15)
simpangan_baku = c(1, 1, 5, 3, 0.1)

# Membagi jendela grafik menjadi 2 baris x 3 kolom (5 grafik + 1 kosong)
par(mfrow = c(2, 3))

for (i in 1:5) {
  mu = rata_rata[i]
  sigma = simpangan_baku[i]
  
  # Rentang x disesuaikan per kurva (mean +/- 4*sd) agar bentuk kurva terlihat jelas
  x = seq(mu - 4*sigma, mu + 4*sigma, length = 200)
  y = dnorm(x, mean = mu, sd = sigma)
  
  plot(x, y, type = "l", lwd = 2, col = "blue",
       main = paste0("Kurva ", kurva[i], " (μ=", mu, ", σ=", sigma, ")"),
       xlab = "X", ylab = "Kepadatan (Density)")
}

mu <- c(0, 2, 5, 10, 1.5)      # Nilai Rata-rata (mu)
sigma <- c(1, 1, 5, 3, 0.1)    # Nilai Simpangan Baku (sigma)

x <- seq(-5, 20, length.out = 500)

plot(x, dnorm(x, mean = mu[1], sd = sigma[1]), type = "n", 
     main = "Simulasi Kurva Distribusi Normal Tertumpuk",
     xlab = "X", ylab = "Densitas", ylim = c(0, 4), xlim = c(-4, 18))

warna <- c("blue", "red", "green", "purple", "darkorange")

for (i in 1:length(mu)) {
  y <- dnorm(x, mean = mu[i], sd = sigma[i])
  lines(x, y, col = warna[i], lwd = 2)
  polygon(c(x[1], x, x[length(x)]), c(0, y, 0), 
          col = adjustcolor(warna[i], alpha.f = 0.2), 
          border = warna[i])
}

legend("topright", 
       legend = paste("Kurva", 1:5, "(mu =", mu, ", sigma =", sigma, ")"), 
       fill = adjustcolor(warna, alpha.f = 0.2), 
       border = warna, 
       cex = 0.7)