> mtaylor(cos(x*y),[x=a,y=b],3);
 

`+`(cos(`*`(a, `*`(b))), `-`(`*`(sin(`*`(a, `*`(b))), `*`(b, `*`(`+`(x, `-`(a)))))), `-`(`*`(sin(`*`(a, `*`(b))), `*`(a, `*`(`+`(y, `-`(b)))))), `-`(`*`(`/`(1, 2), `*`(cos(`*`(a, `*`(b))), `*`(`^`(b, ... (1)
 

> subs([x=a+h,y=b+k],%);
 

`+`(cos(`*`(a, `*`(b))), `-`(`*`(sin(`*`(a, `*`(b))), `*`(b, `*`(h)))), `-`(`*`(sin(`*`(a, `*`(b))), `*`(a, `*`(k)))), `-`(`*`(`/`(1, 2), `*`(cos(`*`(a, `*`(b))), `*`(`^`(b, 2), `*`(`^`(h, 2)))))), `*... (2)
 

>