两种思路：
第一钟是拟合b1，b2两条曲线，然后利用monte carlo投点法计算积分值。
第二种直接利用给定的数据，结合数值积分里面的积分转化为级数求和的方法，两种均可行。

我现在查到的结果是先用三次样条插值，再用数值积分，现在正在写程序，不造对不？

我现在查到的结果是先用三次样条插值，再用数值积分，现在正在写程序，不造对不？

嗯，可行

library(RODBC)

z<-odbcConnectExcel2007("F:\\***\\ESUN2.xlsx")

foo<-sqlFetch(z,'band1')

attach(foo)

close(z)

s<-spline(waveLenght,responseFactor,xout=seq(430,530,0.5))

m<-s$x

n=s$y

f<-function(m) return(n)

source('E:\\desktop\\code\\r\\simpson.R')

simpson(f,440,520)

我这个代码写完以后run不起来，老是说内存不足“Error: cannot allocate vector of size 16.0 Gb”
电脑内存还挺大的唉，24G的，这是怎么回事呢