Tag: 积分

快速计算R中的双积分

我正在寻找一种比双重积分更快的解决方案 integrate(function(y) { sapply(y, function(y) { integrate(function(x) myfun(x,y), llim, ulim)$value }) }, llim, ulim) 用例如 myfun <- function(x,y) cos(x+y) llim <- -0.5 ulim <- 0.5 我发现了一篇旧文件 ,提到了一个名为quad2d的FORTRAN程序,但除了matlab的其他一些帮助页面之外我找不到其他任何东西。 所以我正在寻找一个可以快速进行双积分的C或FORTRAN库(即没有sapply循环),并且可以从R调用。所有的想法都非常受欢迎,只要它们都是GPL兼容的。 如果解决方案涉及从R已经附带的库中调用其他function,我也很乐意听取他们的意见。