# Spike integral

Nick Hale, October 2010

We demonstrate the adaptive capabilities of Chebfun by integrating the 'spike function'

```f = @(x) sech(10*(x-0.2)).^2 + sech(100*(x-0.4)).^4 + ...
sech(1000*(x-0.6)).^6 + sech(1000*(x-0.8)).^8;
```

(which appears as F21F in ) over [0 1].

We create a Chebfun representation and plot the function, increasing 'minsamples' so that the spikes are not missed by an overly coarse initial sample:

```ff = chebfun(f,[0 1], 'splitting','on','minsamples',129);
plot(ff,'b','linewidth',1.6,'numpts',1e4)
title('Spike function','FontSize',16)
``` Now we compute the integral. In order to estimate the time for this computation, we create the chebfun again without plotting it.

```tic
ff = chebfun(f,[0 1], 'splitting','on','minsamples',129);
sum(ff)
```
```ans =
0.211717021214835
```

Time for creating this chebfun and integrating it:

```toc
```
```Elapsed time is 0.357464 seconds.
```

References:

 D. K. Kahaner, "Comparison of numerical quadrature formulas", in J. R. Rice, ed., Mathematical Software, Academic Press, 1971, 229-259.