Can I use taylor(...) with extra massaging to generate a formal expansion like; f(x + e) --> f(x) + e* 'diff( f(x), x, 1) + e^2 * 'diff( f(x), x, 2) / 2 + ... ?? or can this be done with pdiff ? What about 2D? Ted Woollett