Summing integrands in a loop
显示 更早的评论
I want to integrate the following expression:
A = x + x^2 + x^3;
How can one do it using a for loop? For instance, I am trying something like as follows:
for k = 1:1:3
fun1{k} = @(x)(x.^k);
end
A_sum = cellfun(@sum,fun1{k});
int = integral(A_sum, 0, 1)
But this results in an error.
I know that the above expression can be evaluated easily without using a for loop. But the actual expression that I want to evaluate contains the summation of about 200 complex terms.
采纳的回答
Star Strider
2021-3-16
The loop is not necessary.
Try this:
fun1 = @(x,k)(x.^k);
k = 1:3;
intv = integral(@(x)sum(fun1(x,k)), 0, 1, 'ArrayValued',1)
Check = integral(@(x) x + x.^2 + x.^3, 0, 1, 'ArrayValued',1) % Check That Vectorisation Worked
Both produce the same result:
intv =
1.083333333333333
Check =
1.083333333333333
.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
