clear
all
;close
;clc;
syms
x1
x2
t
;
f=2*x1^2+x2^2-x1*x2;
f_grad=[diff(f,x1);diff(f,x2)];
X0=[1;1];
n=10;
epsonal=0.01;
fx=inline(f);
fx_grad=inline(f_grad);
X=X0;
fx0=fx(X(1),X(2));
fx0_grad=fx_grad(X(1),X(2));
while
1