- #1
Physics_rocks
- 12
- 0
Hi guys ,
I need to write a code in order to calculate Simpson's Rule.
Here is my code , what's wrong with it ?
Here is the link to the actual formula
"[URL
Thanks a lot :)
I need to write a code in order to calculate Simpson's Rule.
Here is my code , what's wrong with it ?
PHP:
function [integral]=simp(f,a,b,n)
if mod(n,2)~=0
warning('n must be even!')
else
h=(b-a)/n;
integral=0;
for i=1:1:n-1
x=a+i*h;
if (mod(i,2)==0)
integral=integral+2*eval(f);
else
integral=integral+4*eval(f);
end
end
x=a;
integral=integral+eval(f);
x=a+n*h;
integral=integral+eval(f);
integral=(h/3)*integral;
end
"[URL
Thanks a lot :)
Last edited by a moderator: