Hi Chris! I was playing with the adaptive Simpson function and found some problems evaluating numerically a (well-behaved) integral. I've a couple of questions: 1) Why the decision to *not* couple the precision to the value of Digits? Coupling it to Digits would make it available in Ginsh. Wouldn't it do the right thing in almost all cases? 2) Are variations of Simpson's rule really attractive for numerical evaluation? Have you considered/tried something funkier (maybe Romberg's method) or is that a particularly smart implementation? 3) Probably an oversight: the last argument to the adaptivesimpson function is unused. Instead, the code always uses the default. Best wishes -richy. -- Richard B. Kreckel <http://www.ginac.de/~kreckel/>