2. 代数多項式近似:ぱらぱらめくる『A Short Course on Approximation Theory』
library(pracma) my.function<-function(x){ sin(x)+cos(x^2) } n<-1:20 precision<-rep(0,length(n)) for(i in 1:length(n)){ out<-polyApprox(my.function, -pi/2, pi/2, n[i]) precision[i]<-out$estim.prec } par(mfcol=c(1,2)) plot(n,precision) f <- polyApprox(my.function, -pi/2, pi/2, 100)$f x <- seq(-pi/2, pi/2, length.out = 100) y <- my.function(x) - f(x) plot(x, y, type = "l", col = "blue") grid() par(mfcol=c(1,1))