For comparison with other CAS choose from: Axiom Derive Macsyma Maple Mathematica
taylor (e**x, x, 0, 4);
load_package taylor;
1 2 1 3 1 4 5 1 + x + ---*x + ---*x + ----*x + O(x ) 2 6 24
taylor (e**(x+y), x, 0, 2, y, 0, 2);
1 2 3 3 1 + y + ---*y + x + y*x + (4 terms) + O(x ,y ) 2
taylorcombine (ws**2);
2 3 3 1 + 2*y + 2*y + 2*x + 4*y*x + (4 terms) + O(x ,y )