Практикум по методам решения задачи Коши для систем ОДУ. Даутов P.З. - 19 стр.

§ 3. Программирование методов Рунге-Кутта при решении систем ОДУ.    19


Выше за символом процента указан комментарий; показано, как со-
ответствующая команда должна выглядеть в программе MatLab (это
связано со спецификой передачи в функции имен функций).
    Для каждой задачи Коши должна быть написана своя программа
вычисления правой части; макет этой программы выглядит следую-
щим образом:
f u n c t i o n dy=f c t ( t , y )
n=numel ( y ) ;
dy=z e r o s ( 1 , n ) ;

dy (1)= . . . ;
dy (2)= . . . ;
... .......
dy ( n)= . . . ;

Здесь вместо многоточия в dy(i) = ...; надо записать выражение
функции fi (t, y) в терминах переменных t, y(1), . . . , y(n). Конечно,
при программировании на языках типа Си, Паскаль и т.д. часть про-
граммы, соответствующий команде dy = zeros(1, n), является лиш-
ней.