CONTOH 2.3.1a dan 2.3.1b

CONTOH 2.3.1a

contoh.2.3.1a <- function(age, n) {
  exp.kontinu <- function(t) {
  (1-(0.01*(age+t))^2)/(1-(0.01*age)^2) }
  exp2.kontinu <- function(t) {
  t*(1-(0.01*(age+t))^2)/(1-(0.01*age)^2) }
  output.exp.kontinu <- integrate(exp.kontinu, 0, n)$value
  output.exp2.kontinu <- 2*integrate(exp2.kontinu, 0, n)$value
  output.var.kontinu <- output.exp2.kontinu-(output.exp.kontinu^2)
  cat("Nilai Ekspetasi Kontinu = ", output.exp.kontinu, "\n")
  cat("Nilai (Ekspetasi^2) Kontinu = ", output.exp2.kontinu, "\n")
  cat("Nilai Variansi Kontinu = ", output.var.kontinu, "\n")}

contoh.2.3.1a(age = 30, n=50)

## Nilai Ekspetasi Kontinu =  37.17949 
## Nilai (Ekspetasi^2) Kontinu =  1607.143 
## Nilai Variansi Kontinu =  224.8286

CONTOH 2.3.1b

contoh.2.3.1b <- function(age, n) {
  k <- 1:n
  
  exp.diskret <- (1 - (0.01 * (age + k))^2) / (1 - (0.01 * age)^2)
  exp2.diskret <- (2 * k - 1) * (1 - (0.01 * (age + k))^2) / 
                  (1 - (0.01 * age)^2)
  output.exp.diskret <- sum(exp.diskret)
  output.exp2.diskret <- sum(exp2.diskret)
  output.var.diskret <- output.exp2.diskret-
                        (output.exp.diskret^2)
  cat("Nilai Ekspektasi diskret =", 
      output.exp.diskret, "\n")
  cat("Nilai (Ekspektasi^2) diskret =", 
      output.exp2.diskret, "\n")
  cat("Nilai Variansi diskret =", 
      output.var.diskret, "\n")}

contoh.2.3.1b(age=30, n=50)

## Nilai Ekspektasi diskret = 36.87637 
## Nilai (Ekspektasi^2) diskret = 1589.799 
## Nilai Variansi diskret = 229.9325