TUGAS 3 DI RPUBS

#1a

mclaurin_exp <- function(x,n){
  sum <- 0
  fact <- 1
  for (i in 0:n){
    pangkat <- x^i
    if(i!=0){
      fact = fact * i
    }
    suku <- pangkat/fact
    sum = sum + suku
  }
  return(sum)
}

mclaurin_exp(-1,5)

#1b

sin_mcl <- function(x){
  hasil = 0
  suku = 1
  k = 0
  while(suku > 10^(-5)){
  p <- 2*k +1
  fact <- 1
  for (i in 1:p){
    fact <- fact * i
  }
  suku <- x^p / fact
  
  if (suku < 10^(-5)){
    break
  }
  hasil <- hasil + ((-1)^k * suku)
  
  k = k + 1
  }
  return(hasil)
}

print(sin_mcl(pi/6))

#2

D <- c(4.1, 4.9, 6.2, 6.9, 6.8, 4.4, 5.7, 5.8, 6.9, 4.7, 6, 4.9)
k = 3

result <- sapply(3:length(D), function(i){
  summ <- D[i]+D[i-1]+D[i-2]
  mean <- summ/k
})
result