Cody

Problem 1031. Composite Trapezoidal Rule for Numeric Integration

Solution 161743

Submitted on 13 Nov 2012
This solution is locked. To view this solution, you need to provide a solution of the same size or smaller.

Test Suite

Test Status Code Input and Output
1   Pass
%% Simple Trapezoidal Rule p=[2 0 -4 0 -1 1]; a=-2; b=4; f = @(x) polyval(p,x); assert(isequal(trapezoidal_rule(f,a,b,1),5280))

``` ```

2   Fail
%% Composite Trapezoidal Rule for 2 intervals p=[2 0 -4 0 -1 1]; a=-2; b=4; f = @(x) polyval(p,x); assert(isequal(trapezoidal_rule(f,a,b,2),2634))

```Error: Assertion failed. ```

3   Fail
%% Composite Trapezoidal Rule for 4 intervals p=[2 0 -4 0 -1 1]; a=-2; b=4; f = @(x) polyval(p,x); assert(isequal(trapezoidal_rule(f,a,b,4),1516.875))

```Error: Assertion failed. ```

4   Fail
%% Exact analytical comparison p=[2 0 -4 0 -1 1]; a=-2; b=4; f = @(x) polyval(p,x); P=polyint(p); I_correct=polyval(P,b)-polyval(P,a); I=trapezoidal_rule(f,a,b); assert(abs(I-I_correct)<1)

```Error: Not enough input arguments. ```

5   Fail
%% Exact analytical comparison--higher tolerance p=[2 0 -4 0 -1 1]; a=-2; b=4; f = @(x) polyval(p,x); I = trapezoidal_rule(f,a,b,1000); P=polyint(p); I_correct=polyval(P,b)-polyval(P,a); assert(abs(trapezoidal_rule(f,a,b,1000)-I_correct)<1e-1)

```Error: Assertion failed. ```