Если я правильно понял задачу, тут простая математика. Условие состоит из двух уравнений, в которых всего две неизвестные: y[i]
и z[i]
.
y[i] = y[i-1]*d1 +x[i]*d2 + z[i]*d3
z[i] = z[i-1]*d4 + y[i]*d5
А такие вещи решаются очень легко, методом подстановки выражения из первого уравнения во второе.
Следите за руками:
- Вместо
y[i]
во второе уравнение подставляем его значение:
z[i] = z[i-1]*d4 + d5*(y[i-1]*d1 +x[i]*d2 + z[i]*d3)
- Делим обе части уравнения на
d5
:
z[i]/d5 = z[i-1]*d4/d5 + y[i-1]*d1 +x[i]*d2 + z[i]*d3
- Переносим
z[i]*d3
в левую часть:
z[i]/d5 - z[i]*d3 = z[i-1]*d4/d5 + y[i-1]*d1 +x[i]*d2
- Умножаем обратно на
d5
:
z[i] - z[i]*d3*d5 = z[i-1]*d4 + (y[i-1]*d1 +x[i]*d2)*d5
- В левой части выносим
z[i]
за скобки:
z[i]*(1 - d3*d5) = z[i-1]*d4 + (y[i-1]*d1 +x[i]*d2)*d5
- Делим обе части на
(1 - d3*d5)
:
z[i] = (z[i-1]*d4 + (y[i-1]*d1 +x[i]*d2)*d5) / (1 - d3*d5)
И вот теперь мы знаем как вычислить z[i]
, которое не зависит от y[i]
.
Следовательно, решение:
// инициализация
for i := 0 to n do begin
x[i] := ???
end;
y[0] := ???
z[0] := ???
// вычисление
for i := 1 to n do begin
z[i] := (z[i-1]*d4 + (y[i-1]*d1 +x[i]*d2)*d5) / (1 - d3*d5);
y[i] := y[i-1]*d1 +x[i]*d2 + z[i]*d3;
end;