共轭梯度法求解二元函数极小值

x0=[1;1];%初始点xk=x0;g0=zeros(2,1);g1=zeros(2,1);g2=zeros(2,1);d0=zeros(2,1);d1=zeros(2,1);symsx1symsx2xk1=xk;%计算xk点的梯度及梯度值fun=fun(x1,x2);fx1=diff(fun,'x1');fx2=diff(fun,'x2');fun=inline(fun);fx1=inline(fx1);fx2=inline(fx2);funval=feval(fun,xk1(1),xk1(2));gradx1=feval(fx