Simpsons Rule: With for loops

12 views (last 30 days)
Eric
Eric on 10 May 2015
Answered: Malek Arnous on 20 Dec 2017
Hi, So I have a question where I have to use Simpsons rule to integrate (1-x^3)*sin(x) + exp(x^2/20) between -1 and 4 with 20 intervals. The function has 4 inputs, f(x), a,b (start and end points) and n intervals
I know that I can make this code simpler with the sum function but unfortunately I have to use loops for this exercise.
My code looks like this:
function integral = simpsonsrule(f,a,b,n)
h = (b-a)/n;
x = linspace(a,b,n);
x4=0;
x2=0;
for j=2:2:b
x4 = x4 + f(x4);
end
for k=3:2:b
x2= x2 + f(x2);
end
integral = (h/3)*(f(a)+ f(b) + 4*(x4)+ 2*(x2));
end
And I'm calling it like this:
clear;
integral = simpsonsrule((1-x.^3)*sin(x) + exp(x.^2/20),-1,4,20)
But I'm getting the error: Undefined function or variable 'x'. but haven't I defined it with x=linspace(a,b,n)?

Accepted Answer

Walter Roberson
Walter Roberson on 10 May 2015
integral = simpsonsrule(@(x) (1-x.^3)*sin(x) + exp(x.^2/20),-1,4,20)
You need the @(x) to make an anonymous function
  8 Comments
Walter Roberson
Walter Roberson on 10 May 2015
What kind of n are you using?
Eric
Eric on 10 May 2015
I've been using 20,50,100,500,1000.

Sign in to comment.

More Answers (1)

Malek Arnous
Malek Arnous on 20 Dec 2017
yes

Categories

Find more on Numerical Integration and Differential Equations 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!