How to plot a function containing an integral?
显示 更早的评论
Hello all! I've got very simple code:
function y=sc()
a=2;
x=linspace(0,100);
fun1 =quadgk(@(m) m.*a .*x, 0, 2);
y = plot(x,fun1);
end
In the specific case above the value of the integral =2, and also a=2, so one could expect the linear plot as a result o function sc. Instead there is an error: error using .*, matrix dimensions must agree. The fun1 gives exact values, when user provides a and x values. Is there a simple way to create a plot(x,fun1) without making any loops?
采纳的回答
Star Strider
2016-3-31
I had to delve back in the documentation, but I found the quadv function that will do what you want without a loop:
fun1 = quadv(@(m) m.*a .*x, 0, 2);
That works.
If you have R2012a or later, use the integral function with the 'ArrayValued' option:
a=2;
x=linspace(0,100);
fun1 = integral(@(m) m.*a .*x, 0, 2, 'ArrayValued',1);
y = plot(x,fun1);
It produces the same results as quadv for this problem. Otherwise, I doubt there is any way to do it with any of the numerical integration routines without using a loop.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Annotations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
