left riemann sum integration

38 views (last 30 days)
Harjot Purewal
Harjot Purewal on 4 Dec 2016
Commented: Sam Zavala on 5 Apr 2021
so I have to make a numerical integrator using a few methods, just made a left Riemann sum calculator but I seem to be missing something
my code seems to be more accurate when I multiply my y array by dx in the for loop, but should I be multiplying dx at each iteration or at the very end?
also is it just my code or does Riemann sum require a lot of iteration to be accurate?
% code
clear;
clc;
x1 = 1; x2 = 3;
f=@(x) 2*x.^5 - 3*x.^2 - 5;
c=1;
n=500;
dx=(x2-x1)/n;
for i=0:dx:n-1
y(c)=dx*f(1+i*dx);
c=c+1;
s=sum(y)*dx
z=integral(f,x1,x2)
end

Accepted Answer

Torsten
Torsten on 5 Dec 2016
Edited: Torsten on 6 Dec 2016
Try
x1 = 1;
x2 = 3;
f=@(x) 2*x.^5 - 3*x.^2 - 5;
n=500;
dx=(x2-x1)/n;
summe=0.0;
for i=1:500
summe=summe+f(x1+dx*(i-1));
end
summe=summe*dx;
z=integral(f,x1,x2);
Best wishes
Torsten.
  2 Comments
David Araujo
David Araujo on 20 Oct 2020
How can i plot this?
Sam Zavala
Sam Zavala on 5 Apr 2021
I hope this is somewhat of a helpful answer but what we have been doing is build an array using a for loop and our time as telling us how many times to repeat it. Then you simply plot the array :)

Sign in to comment.

More Answers (0)

Categories

Find more on Mathematics in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!