(If in scope) add extrema

[eprovst]

It would be useful to provide implementations for minimum, maximum and extrema. I guess (going by Chebfun) this would effectively boil down to:

extrema(f) = extrema(map(x->f[x], [-1; findall(iszero, diff(f)); 1]))

or something similar.