#1a
macluarin<-function(x, n){
faktori<-function(x){
faktor<-1
for(i in x:1){
faktor<-faktor*i
}
return (faktor)
}
hasil<-0
for(i in 1:n){
hasil<-hasil+(x^i/faktori(i))
}
return(1+hasil)
}
macluarin(-1, 5)
## [1] 0.3666667
#1b
macluarin_sin<-function(x){
faktori<-function(x){
faktor<-1
for(j in x:1){
faktor<-faktor*j
}
return (faktor)
}
hasil<-0
i<-0
suku<-x
while(abs(suku)>10^-5){
hasil<-hasil+suku
i <- i+1
suku<-(-1)^i*x^(2*i+1)/faktori(2*i+1)
}
return(hasil)
}
macluarin_sin(pi/6)
## [1] 0.5000021
#2
data<-c(4.1, 4.9, 6.2, 6.9, 6.8, 4.4, 5.7, 5.8, 6.9, 4.7, 6, 4.9)
rata_bergerak<-function(data, k){
n<-length(data)
hasil <- rep(NA, n)
for(t in k:n){
jumlah<-0
for(i in t:(t-k+1)){
jumlah<-jumlah+data[i]
}
hasil[t]<-jumlah/k
}
return(hasil)
}
mt<-rata_bergerak(data, 3)
tabel<-data.frame(
t=1:length(data),
data=data,
Mt=mt
)
tabel