Я хочу вычислить двукратный повторяющийся интеграл (не двойной интеграл) в R, например,
куда
На практике и f(x), и g(x) довольно сложны, но для проведения эксперимента упростим g(x)=1 и f(x)=cos(x), в R я использую integrate
для вычисления:
> phi = function(x){integrate(function(x){cos(x)},lower=x,upper=3)[["value"]]^2}
> foldintegral = integrate(phi,lower=0,upper=3)
И я получил это сообщение об ошибке:
Error in integrate(phi, lower = 0, upper = 3) :
evaluation of function gave a result of wrong length
Любая идея, как это сделать?