Info
This question is closed. Reopen it to edit or answer.
can someone explain me this code please?
1 view (last 30 days)
Show older comments
clc;
clear all;
close all;
f=@(x)x.^3;
a=1;b=2;
n=1.25;
h=(b-a)/n;
p=0;
for i=a:b
p=p+1;
x(p)=i;
y(p)=i.^3;
end
l=length(x);
x
y
answer=(h/2)*((y(1)+y(l))+2*(sum(y)-y(1)-y(l)))
0 Comments
Answers (1)
Roger Stafford
on 27 May 2015
This looks like a misguided attempt to use trapezoidal approximation to find the integral of y = x^3 over an interval in x. The values assigned to 'a', 'b', and 'h' are inappropriate for such an interpretation however.
Anyhow, there is a much easier method for such an integral using calculus:
answer = (b^4-a^4)/4;
Why do numerical approximation when such an easy formula exists with the exact answer?
0 Comments
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!