Why do my loops only run once

1 次查看(过去 30 天)
function I=simpson13(f,a,b,n)
h=(b-a)/n; % length of the interval
total=0;
x=a;
fa=f;
x=b;
fb=f;
for i=1:2:n-1;
x=i*h;
fn=f+total;
total=fn;
end
for j=2:2:n-2;
x=j*h;
fm=f+total;
total=fm;
end
I=(h/3)*(fa+fb+4*fn+2*fm);
end
  1 个评论
Raquel Terrer van Gool
编辑:Geoff Hayes 2020-4-14
I have n=6:
x=0;
n=6; % n is the number of intervals and must be > 0 and even
a=0; b=pi;% a is the lower bound and b is the upper bound
f=5+13*sin(x); % Function
I=simpson13(f,a,b,n)

请先登录,再进行评论。

采纳的回答

Geoff Hayes
Geoff Hayes 2020-4-14
Raquel - the MATLAB debugger is helpful in cases like this. If you put a breakpoint in your simpson13 you would notice a problem with the input parameter f. Look closely at how it is being assigned
x=0;
n=6; % n is the number of intervals and must be > 0 and even
a=0; b=pi;% a is the lower bound and b is the upper bound
f=5+13*sin(x); % Function
f is not a function in this case but is a scalar value (5) since x is zero. If you want f to be an anonymous function then you need to assign it as
f = @(x)5+13*sin(x);
where nthe @ operator creates the handle, and the parentheses () immediately after the @ operator include the function input arguments. You will need to change how this function is used in your code as well. For example,
fa=f(a);
fb=f(b);
etc.

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Mathematics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by