2 курс / Численные методы / числ_методы / Занятие7

Занятие 7

Задание:

1. Решить дифференциальное уравнение: y′=sin( y t) +0.1y2 ,

y(0)=0.5 методом Рунге-Кутты Фельберга 4-5 порядка (метод rkf), используя процедуру решения дифференциальных уравнений ode в системе SciLab (она описана в главе 8 учебника по SciLab). Вывести график решения на экран.

2. Решить систему дифференциальных уравнений методом Рунге-Кутты Фельберга 4-5 порядка (метод rkf):

где t — время, r=r(t), f=f(t) и α - положительная константа, R =5000, α=0.1, r0=300 и f0=150. Вывести график решения на экран.

3. Реализовать собственную компьютерную программу решения уравнения данного в пункте 1 методом РунгеКутты Фельберга 4-5 порядка не в Scilab, а на языке высокого уровня (Pascal, Basic или др.). Алгоритм РунгеКутты Фельберга был описан на лекции.