clc;
syms?x1?x2?bc?a?b
f=2*x1^2+6*x2^2+2*x1*x2+2*x1+3*x2+3;
e=0.01;
a=10;
b=1;
fx1=diff（f‘x1‘）;
fx2=diff（f‘x2‘）;
g=[fx1fx2]‘;
x1=a;
x2=b;
g1=subs（g）
x0=[x1?x2]‘;
k=0;
norm（g1）
while?（abs（norm（g1））>=e）
????