& gt& gtsyms x y；f=(x^2-2*x)*exp(-x^2-y^2-x*y)；
& gt& gtF=maple('mtaylor '，f，'[x，y]'，8)
& gt& gtF=collect(F，x)
F =
1/3*x^7+(y+1/2)*x^6+(y-1+2*y^2)*x^5+(-1-2*y+3/2*y^2+7/3*y^3)*x^4+(2*y^4-y+y^3-3*y^2+2)*x^3+(2*y+1-y^2+y^5-2*y^3+1/2*y^4)*x^2+(-2+1/3*y^6-y^4+2*y^2)*x
多元函數的泰勒展開
f(x1，x2，…xn)=f(a1，a2，…an)+
[(x 1-a 1)]f(a 1，…，an)+
+…+
+…
其中a1，…an是泰勒級數展開的中心點。為了避免分歧，這裏的公式應該理解為先取f函數的導數，再取點a1，a2，…an的導函數值。Matlab沒有直接提供計算函數，但是可以通過調用Maple語言中的mtaylor()函數直接得到。呼叫格式是:
F = maple ('mtaylor '，f，'[x1，…，xn]'，k)按原點展開。
F = maple ('mtaylor '，f，'[x1 = a1，…，xn = an]'，k)根據(a1，…，an)展開。
其中k-1為最高階展開式，f為原多元函數。註意函數調用時自變量的引號不能省略，調用格式會原封不動的傳遞給Maple。