計算不能に備える

  • Rには関数を積分する関数integrate()がある
    • ときに計算不能になる
  • そんなときは、別な方法が逃げ道になるかもしれない(ならないかもしれない)ので、確認しておく
install.packages("sfsmisc")
library(sfsmisc)
help(integrate.xy)
x <- 1:4
 integrate.xy(x, exp(x))
 print(exp(4) - exp(1), digits = 10) # the true integral

 for(n in c(10, 20,50,100, 200)) {
   x <- seq(1,4, len = n)
   cat(formatC(n,wid=4), formatC(integrate.xy(x, exp(x)), dig = 9),"\n")
 }