library(ggplot2)
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
# cho p = 0.4 , cho số lần thử là 20 lần
# vẽ 2 biểu đồ
n <- 20
p <- 0.4
par(mfrow = c(1,2))
# Phân phối xác suất (PMF)
plot(0:n, dbinom(0:n, n, p),
xlab = "x",
ylab = "P(X = x)",
main = "Binomial Distribution",
col = "darkblue")
# Phân phối tích lũy (CDF)
plot(0:n, pbinom(0:n, n, p),
xlab = "x",
ylab = expression(P(X <= x)),
main = "Binomial Cumulative Distribution",
col = "red")

# Số lượng theo đề bài
lang <- 8
goc_moi_lang <- 8
nha_moi_goc <- 8
vo_moi_ong <- 8
con_moi_vo <- 8
# Số nhà trong 1 làng
so_nha_1_lang <- goc_moi_lang * nha_moi_goc
# Số người trong 1 nhà
nguoi_1_nha <- 1 + vo_moi_ong + vo_moi_ong * con_moi_vo
# 1 ông + 8 vợ + 64 con
# Số người trong 1 làng
nguoi_1_lang <- so_nha_1_lang * nguoi_1_nha
# Số người trong 8 làng
tong_nguoi <- lang * nguoi_1_lang
tong_nguoi
## [1] 37376
lang<-tong_nguoi*8
lang
## [1] 299008