Моделирование технических систем. Дьячков Ю.А - 189 стр.

     Procedure rk4;
     begin
      fst;
      for i:=nn to nd do
      begin
      y[i]:=yy[i]+h*der[i]/2;
      rk[i]:=der[i];
      end;
      fst;
      for i:=nn to nd do
      begin.
      y[i]:=yy[i]+h*der[i]/2;
      rk[i]:=rk[i]+2*der[i];
      end;
      fst;
      for i:=nn to nd do
      begin
      y[i]:=yy[i]+h*der[i];
      rk[i]:=rk[i]+2*der[i];
      end;
      fst;
      for i:=nn to nd do
      begin
      rk[i]:=(rk[i]+der[i])*h;
      y[i]:=yy[i]+rk[i]/6
      end;
     end.



         14.3 Многошаговые численные методы
     В ходе решения задачи появляется дополнительная информа-
ция о рассчитанных точках траектории движения системы, которая
не используется в одношаговых методах. Применение же многоша-
говых методов основано на использовании информации о предыду-
щих вычислениях. Для этого используют две формулы: прогноза и
коррекции. Поэтому такие методы известны под названием методов
прогноза и коррекции.
     В отличие от одношаговых данные методы не обладают свой-
ством «самостартования». Поэтому при их использовании началь-
ные точки расчета определяются одношаговыми методами.

                                 189