Algorithm AS 133: Optimization of One-Dimensional Multimodal Functions