%%坐標輪換法
clear
clc
e=0.04;??%要求精度
xm1=[0?3]‘;???%初始點
e1=[1?0]‘;
e2=[0?1]‘;
syms?x1
syms?x2
syms?s
syms?s2
%%方向1
xm01=xm1+s*e1;
f=（x1-2）^4+（x1-2*x2）^2;
fs1=subs（f{x1x2}{xm01（11）xm01（21）}）;
dfs1=diff（fs1）;
s1=solve（dfs1）;
xs1=imag（s1）;
ns01=size（xs1）;
ns1=ns01（11）;
n01=0;
for?i=1:ns1
????if?xs1（i）==0
????????n01=n01+1;
????????ss1（n01）=s1（i）;
????end
end
s1=max（double（ss1））;
xm01=xm1+s1*e1;
%%%方向2
xm02=xm01+s2*e2;
fs2=subs（f{x1x2}{xm02（11）xm02（21）}）;
dfs2=diff（fs2）;
s02=solve（dfs2）;
xs2=imag（s02）;
ns02=size（xs2）;
ns2=ns02（11）;
n02=0;
for?j=1:ns2
????if?xs2（j）==0
????????n02=n02+1;
????????ss2（n02）=s02（j）;
????end
end
s02=max（double（ss2））;
xm02=xm01+s02*e2;
xm2=xm02;
%迭代條件判斷
w=xm2-xm1;
t=1;??%迭代次數
while?norm（w）>=e
????xm1=xm2;
?