X = -1:.1:3;
Y = -2:.1:2;
[x, y] = meshgrid(X, Y);
z = exp(2*x - 3*y);
T2 = exp(2)*(1 + 2*(x - 1) - 3*y + ...
    2*(x - 1).^2 - 6*(x - 1).*y + 9/2*y.^2);
figure
hold on
xlabel('x')
ylabel('y')
zlabel('z')
surf(x, y, z)
surf(x, y, T2)