Trapezoidal rule in a for loop

Question

Mathias Pettersen 2019-9-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/481777-trapezoidal-rule-in-a-for-loop

评论： James Tursa 2019-9-24

采纳的回答： James Tursa

d

2 个评论
显示无隐藏无

awezmm 2019-9-23

Are you getting an error or the values are just not coming out right

James Tursa 2019-9-24

Mathias, it is considered rude to delete your question after it has received answers. This also can render this thread useless to other readers who might be having similar problems.

请先登录，再进行评论。

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

Answer 1

James Tursa 2019-9-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/481777-trapezoidal-rule-in-a-for-loop#answer_393206

编辑：James Tursa 2019-9-24

在 MATLAB Online 中打开

You made a good start, but there are several problems with your summation. You don't accumulate into your trapez variable, you don't use your index i in your function evaluation, you have your endpoints evaluated multiple times inside your loop, etc.

It would probably be best for you to code the logic up exactly as it appears in the formula. So, first do the summation. Then add the endpoints. Then multiply by (b-a)/n. Also, why not make your life simpler and use the same index variable k that is used in the formula? So, the outline would be this:

trapez = 0;  % Initialize your summation to 0
for k = 1: n - 1 % use k as the index so it looks like the formula
    trapez = trapez + calculateFormula(_____); % you fill in the value here based on formula
end
trapez = _____ + trapez + _____; % you add in the end points here
trapez = _____ * trapez; % you fill in the multiplying factor here

Try to fill in the blanks above to see if you can get it to work. If you have more problems, come back and ask for more help with the specific problems you are having.

2 个评论
显示无隐藏无

Mathias Pettersen 2019-9-24

在 MATLAB Online 中打开

thanks so much. Does this look right?

a = -3;
b = 3;
n = 10;
trapez = 0;
for k = 1: n - 1
    
trapez = trapez +  sum(calculateFormula(a + (k.*(b-a)./n)) + calculateFormula(b)./2);
end
trapez = trapez + calculateFormula(a)/2;
trapez = trapez * ((b-a)/n)

James Tursa 2019-9-24

编辑：James Tursa 2019-9-24

No, it doesn't. But you are getting closer. You don't need the sum( ) function inside your loop since you are in fact doing the summation manually. Also, the f(b)/2 is one of the endpoints that gets added after the loop (look at your formula closely and you will see that it is not intended to be inside the summation).

请先登录，再进行评论。