函数x=FX steep (f，e，a，b ) )。
x2=b；
q=Fxhesson(f，x1，x2 )；
FX1=diff(f，‘x1’)；
FX2=diff(f，‘x2’)；
d=- g1；
while(ABS(norm ) G1 ) (e ) ) ) ) ) ) ) ) ) )。)
lamda=(-d ) ) d/(-d ) ) q*d )；
a=[ 1，0 ] * x0；
b=[ 0，1 ] * x0；
x2=b；
d=- g1；
结束
x=temp；
结束
Functionq=Fxhesson(f，a，b ) ) ) )。
x2=b；
FX1=diff(f，1，‘x1’)；
FX2=diff(f，1，‘x2‘)；
FX 1x1=diff (fx1，1，‘x1’)；
fx1x2=diff (fx1，1，‘x2’)；
fx2x2=diff (fx2，1，‘x2‘)；
fx2x1=diff (fx2，1，‘x1’)；
FX1x1=Subs(FX1x1)；
FX1x2=Subs(FX1x2)；
FX2x1=Subs(FX2x1)；
FX2x2=Subs(FX2x2)；
Q=[fx1x1，fx1x2； fx2x1、FX2x2)；
结束
syms x1 x2；
X=[x1，x2]；
FX=x(1) ) 24*x ) )2) ) 2；
z=FXsteep(FX，0.0001，1，1 ) ) ) ) ) ) )。
ans
1
1
48/65
-3/65
36/325
36/325
1728/21125
-108/21125
1296/105625
1296/105625
62208/6865625
-3888/6865625
46656/34328125
46656/34328125
2239488/2231328125
-139968/2231328125
1679616/11156640625
1679616/11156640625
80621568/725181640625
-5038848/725181640625
60466176/3625908203125
60466176/3625908203125
2902376448/235684033203125
-181398528/235684033203125