>
mtaylor(cos(x*y),[x=a,y=b],3);
(1)
>
subs([x=a+h,y=b+k],%);
(2)
>