simp {broman} | R Documentation |

## Numerical integration

### Description

Perform numerical integration by Simpson's rule or the trapezoidal rule.

### Usage

```
simp(f, a, b, tol = 0.00000001, max.step = 1000, ...)
```

### Arguments

`f` |
The integrand; must be a vectorized function. |

`a` |
Lower limit of integration. |

`b` |
Upper limit of integration. |

`tol` |
Tolerance for choosing the number of grid points. |

`max.step` |
Log base 2 of the total number of grid points. |

`...` |
Other arguments passed to the integrand, |

### Details

Iterately doubles the number of grid points for the numerical
integral, stopping when the integral decreases by less than
`tol`

.

### Value

The integral of `f`

from `a`

to `b`

.

### See Also

### Examples

```
f <- function(x) x*x*(1-x)*sin(x*x)
I1 <- trap(f,0,2)
I2 <- simp(f,0,2)
```

[Package

*broman*version 0.84 Index]