Сделаем первую итерацию:Начальные значения:
\(a=\left(...,0,0,x_o,0,0,...\right)\)
\(S=\left(...,0,0,x_o^2,0,0,...\right)\)
\(Z=\left(...,0,0,x_o^3,0,0,...\right)\)
\(T=\left(...,0,0,0,0,0,...\right)\)
\(Q=\left(...,0,0,0,0,0,...\right)\)
Получим:
\(\bar a=\left(...,0,0,x_o^2+\frac{1}{3}x_o^3,0,0,...\right)\)
При итерациях гармоники так и останутся равными нулю.
Поэтому правильно немного "пошуметь" над начальными значениями.
Нужно алгоритм писать.