Метод наискорейшего спуска

[RocBoy-D]

Здравствуйте! Ни у кого нету кода метода наискорейшего спуска в Mathematica?

Нашел код, но в Maple:
restart;
f:= 1.2*x1^2 + 1.8 * x2^2 - 4 * x1 - 4 * x2:
Определяем производные по двум переменным
df1:=diff(f,x1);
df2:=diff(f,x2);
x1prec:=solve(df1,x1):
x2prec:=solve(df2,x2):
Начальное приближение
X1[1]:=-1:
X2[1]:=-1:
prec:=0.01:
i:=1:
x1rez:=0:
x2rez:=0:
while x1rez=0 do
x1:=X1[i]:
x2:=X2[i]:
x1rez:=`if`(df1<prec,(`if`(df1>-prec,X1[i],0)),0):
x2rez:=`if`(df2<prec,(`if`(df2>-prec,X2[i],0)),0):
X1[i+1]:=`if`(df1>0,X1[i]-prec,X1[i]+prec):
X2[i+1]:=`if`(df2>0,X2[i]-prec,X2[i]+prec):
i:=i+1:
end do:
"Точные значения:";
x1:=x1prec;
x2:=x2prec;
"Точка минимума:";
x1:=x1rez;
x2:=x2rez;
"Значение функции в минимуме:";
f;

Помогите переделать в Mathematica.

Еще нашел реализацию одного шага метода наискорейшего спуска в Mathematica, но никак не соображу, как сделать цикл и счетчик итераций: http://math.fullerton.edu/mathews/n2...Mod_lnk_3.html

Помогите, пожалуйста)