Plotting integral over an interval

5 次查看（过去 30 天）
Angelica Mercy Calip
回答： infinity 2019-7-9
I want to integrate the following over the interval x0 = 0:2:8 . I have tried using arrayfun:
x0 = 0:2:8;
cai= @(theta) ((4.*pi.*x0.*cos(theta).*sin(theta)+8.*pi.*cos(theta).*cos(theta).*sin(theta))./(1+0.5.*x0+cos(theta)));
xval = arrayfun(@(uplim) integral(cai, 0, uplim, 'ArrayValued',true), x0);
plot(x0, xval)
grid
But I got the following error code:
Error using arrayfun
Non-scalar in Uniform output, at index 1, output 1.
Set 'UniformOutput' to false.
Error in PlottingFx (line 6)
xval = arrayfun(@(uplim) integral(cai, 0, uplim, 'ArrayValued',true), x0);
Can anyone help me?

请先登录，再回答此问题。

采纳的回答

infinity
infinity 2019-7-9
Hello,
The way you try to plot the integral may be not correct, for x0 you can compute the integral by using the integral fuction with repsect to theta. This code below may help you understand this well
x0 = 0:2:8;
cai= @(theta) ((4.*pi.*x0.*cos(theta).*sin(theta)...
+8.*pi.*cos(theta).*cos(theta).*sin(theta))...
./(1+0.5.*x0*cos(theta)));
xval = integral(cai, 0, pi,'ArrayValued',true);
plot(x0, xval)
grid
Also, I fix an mismatch between your code and your provided function
In your code:
./(1+0.5.*x0+cos(theta)));
I change it to
./(1+0.5.*x0*cos(theta)));
which is correct with your given equation.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by