# Membuat data dan kurva distribusi normal
x <- seq(-4, 4, length = 100)
y <- dnorm(x, mean = 0, sd = 1)

plot(x, y, type = "l", lwd = 2, col = "blue",
     main = "Kurva Distribusi Normal",
     xlab = "Nilai x",
     ylab = "Kepadatan Peluang")

# Kurva 1: mu = 0, sigma = 1
x1 <- seq(-4, 4, length = 100)
y1 <- dnorm(x1, mean = 0, sd = 1)

# Kurva 2: mu = 10, sigma = 2
x2 <- seq(2, 18, length = 100)
y2 <- dnorm(x2, mean = 10, sd = 2)

# Plot kurva pertama dulu (sebagai dasar)
plot(x1, y1, type = "l", lwd = 2, col = "blue",
     xlim = c(-6, 20), ylim = c(0, 0.45),
     main = "Perbandingan Kurva Distribusi Normal",
     xlab = "Nilai x",
     ylab = "Kepadatan Peluang")

# Tambahkan kurva kedua di atas plot yang sama
lines(x2, y2, lwd = 2, col = "darkgreen")

# Garis mean masing-masing kurva
abline(v = 0, col = "blue", lty = 2)
abline(v = 10, col = "darkgreen", lty = 2)

# Legenda biar jelas kurva mana punya siapa
legend("topright",
       legend = c("μ=0, σ=1", "μ=10, σ=2"),
       col = c("blue", "darkgreen"),
       lwd = 2, lty = 1)

# Data untuk masing-masing kurva
x1 <- seq(-4, 4, length = 100);   y1 <- dnorm(x1, mean = 0,  sd = 1)
x2 <- seq(2, 18, length = 100);   y2 <- dnorm(x2, mean = 10, sd = 2)
x3 <- seq(-2, 10, length = 100);  y3 <- dnorm(x3, mean = 4,  sd = 1.5)
x4 <- seq(10, 30, length = 100);  y4 <- dnorm(x4, mean = 20, sd = 3)

# Plot kurva pertama sebagai dasar
plot(x1, y1, type = "l", lwd = 2, col = "blue",
     xlim = c(-6, 32), ylim = c(0, 0.45),
     main = "Perbandingan Beberapa Kurva Distribusi Normal",
     xlab = "Nilai x",
     ylab = "Kepadatan Peluang")

# Tambahkan kurva-kurva lainnya ke plot yang sama
lines(x2, y2, lwd = 2, col = "darkgreen")
lines(x3, y3, lwd = 2, col = "orange")
lines(x4, y4, lwd = 2, col = "purple")

# Garis mean masing-masing kurva
abline(v = 0,  col = "blue",      lty = 2)
abline(v = 10, col = "darkgreen", lty = 2)
abline(v = 4,  col = "orange",    lty = 2)
abline(v = 20, col = "purple",    lty = 2)

# Legenda
legend("topright",
       legend = c("μ=0, σ=1", "μ=10, σ=2", "μ=4, σ=1.5", "μ=20, σ=3"),
       col = c("blue", "darkgreen", "orange", "purple"),
       lwd = 2, lty = 1)