# intsimp¶

## Purpose¶

Integrates a specified function using Simpson’s method with end correction. A single integral is computed in one function call.

## Format¶

y = intsimp(&f, xlims, tol)
Parameters: &f (scalar) – pointer to the procedure containing the function to be integrated. xlims (2x1 vector) – the limits of x. The first element is the upper limit and the second element is the lower limit. tol (scalar) – The tolerance to be used in testing for convergence y (scalar) – The estimated integral of $$f(x)$$ between $$xlims[1]$$ and $$xlims[2]$$.

## Examples¶

// Function to be integrated
proc f(x);
retp(sin(x));
endp;

// Define limits
xlims = { 1, 0 };

// Integrate using Simpson's method
y = intsimp(&f, xl, 1e-8);
print y;


The code above returns the following:

0.45969769


This will integrate the function between 0 and 1.

intsimp.src