最优化-牛顿法求最优解matlab程序

%%%牛顿法寻优
clc
clear
syms x1
syms x2
f=（x1-2）^4+（x1-2*x2）^2 ;%目标函数
v=[x1 x2];
df=jacobian（fv）;  %一阶偏导
df=df‘;
G=jacobian（dfv）;  %二阶偏导
G=G‘;
es=0.04;
x0=[0 3]‘;
g1=subs（df{x1x2}{x0（11）x0（21）}）; %初始值
G1=subs（G{x1x2}{x0（11）x0（21）}）;  %初始判断值
k=0;