- 284 名前：Exp <- function(x,n) {
k=0:n
sum(x^k/factorial(k))
}
Exp(1,N)
exp(1)
> N=20
> Exp(1,N)
[1] 2.7182818284590451
> exp(1)
[1] 2.7182818284590451
N=10
Sin <- function(x,n){
k=0:n
sum((-1)^k*x^(2*k+1)/factorial(2*k+1))
}
Sin(1,N)
sin(1)
Cos <- function(x,n) {
k=0:n
sum((-1)^k*x^(2*k)/factorial(2*k))
}
Cos(1,N)
cos(1)
> Sin(1,N)
[1] 0.8414709848078965
> sin(1)
[1] 0.8414709848078965
> Cos(1,N)
[1] 0.54030230586813977
> cos(1)
[1] 0.54030230586813977
e で２０，sin cos で１０も計算すれば２０桁くらいの合致は得られるようだ。
マウント猿ってこういうのも手計算するんだろうかなぁ。 [] - [ここ壊れてます]
|
 |
