%%%牛頓法尋優
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;