本題屬於無約束優化問題，用fminunc求解。
需要註意的是，函數顯然最小值為-inf，所以只能是求極大值：
f=@(x)4*(x(1)-x(2))-x(1)^2-x(2)^2;
x=fminunc(@(x)-f(x),[0 0])
ezmesh('4*(x-y)-x^2-y^2')
hold on
plot3(x(1),x(2),f(x),'b.','markersize',20)