用MATLAB的fminunc命令编写脚本求解f(x)=x1^4/4+x2^2/2-x1x2+x1-x2

首先，定义目标函数：

function [f,grad] = myfun(x)
f = x(1)^4/4 + x(2)^2/2 - x(1)*x(2) - x(1) + x(2);
if nargout > 1
    grad = [x(1)^3 - x(2) - 1; x(2) - x(1)];
end

然后，使用fminunc命令求解：

options = optimoptions('fminunc','Display','iter','Algorithm','quasi-newton');
x0 = [0;0];
[x,fval,exitflag,output] = fminunc(@myfun,x0,options);

其中，options是fminunc的选项结构体，可以指定一些算法参数，例如显示迭代信息等；x0是初始点；@myfun是目标函数的句柄，即函数名前加上@符号；x是最优解；fval是最优解对应的函数值；exitflag是退出标志，1表示收敛到最优解，其它值表示未收敛或出现错误；output是输出结构体，包含算法的统计信息。