How to plot graph of Simpson's rule and Trapezoidal rule between error and n?
11 次查看(过去 30 天)
显示 更早的评论
Hi there,
I have written 2 programs of Simpson's rule and Trapezoidal rule to compute the integral of a function f. Now I would like to plot graphs of both rules between the errors and n but I don't know how to do it. Could anyone help me with this problem please?
Here is the code of my program:
a=0.1;
b=0.2;
n=30;
x=linspace(a,b,n);
h=(b-a)/n;
f=inline('x.^4-x.*5+1');
y=x.^4-x.*5+1;
H=h*ones(1,n); %trapezoid
H([1,n])=h/2;
tr=sum (H.*f(x));
tr
H(1:2:n)=2*h/3; %simpson
H(2:2:n)=4*h/3;
H([1,n])=h/3;
simp=sum(H.*f(x));
simp
Thank you in advance!
0 个评论
回答(2 个)
Alan Stevens
2020-12-13
If you use cumsum instead of just sum then you can plot tr and simp against x.
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!