Latihan 2 (pertemuan 3)

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 (Ekspektasi^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 (Ekspektasi^2) kontinu = 1607.143 
## Nilai Variansi kontinu = 224.8286

CONTOH 2.3.1.b

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