Обучение радиального нейрона

[mikemike]

Здравствуйте, мне нужно написать программу для обучения радильного нейрона, собственно никак не могу понять алгоритм.
Входа у нейрона два. У меня все сводится к задаче поиска минимума функции E = $\min \sum_{k=1}^p (\varphi(||X^k - C||) - d)^2 $, где p - количество обучающих выборок, С-центр, d-ожидаемое значение. В качестве функции активации я использую функцию Гаусса. Я вычислил значения градиента
$\frac{\delta F}{\delta c_1} = \frac{e^{-\frac{(x_1-c_1)^2 + (x_2-c_2)^2}{2\sigma ^2}}(e^{(...)} - d)(x_1-c_1)}{\sigma^2} $
с c2 тоже самое
$\frac{\delta F}{\delta \sigma} = \frac{e^{-\frac{(x_1-c_1)^2 + (x_2-c_2)^2}{2\sigma ^2}}(e^{(...)} - d)((x_1-c_1)^2 + (x_2-c_2)^2)}{\sigma^3} $
вот я собтсвенно и не знаю как составить итерационный процесс. А так же как выбрать начальные значения центра и ширины сигма