Решение ДУ в таблице Excel.
Решение ДУ в программе Visual Basic. Метод Эйлера. Sub Metod_E1() x0 = Val(InputBox("Ввод х0", "Начальное значение ")) y0 = Val(InputBox("Ввод y0", "Начальное значение")) h = Val(InputBox("Ввод h", "Шаг интегрирования")) hp = Val(InputBox("Ввод hp", "Шаг печати")) b = Val(InputBox("Ввод b", "Конец интервала")) x = x0: y = y0 i = 12 While x < b f = x + Cos(y / 3.14) y = y + h x = x + h Do While x > hp i = i + 1 Worksheets("лаба6").Cells(i, 1).Value = x Worksheets("лаба6").Cells(i, 2).Value = y hp = hp + 1 Loop Wend End Sub Метод Рунге-Кутта. Sub Runge_Kytta() h = 0.01 hp = 0.1 x0 = Val(InputBox(“Введите начальное значение x0", "Ввод х0")) xk = Val(InputBox(“Введите начальное значение xk", "Ввод xk")) y0 = Val(InputBox("Введите начальное значение y0", "Ввод y0")) x = x0: y = y0 i = 20: j = 1 While x < xk k0 = z(x, y) * h k1 = z(x + h / 2, y + k0 / 2) * h k2 = z(x + h / 2, y + k1 / 2) * h k3 = z(x + h / 2, y + k2) * h y = y + (k0 + 2 * k1 + 2 * k2 + k3) / 6 x = x + hp Do While x > hp i = i + 1 Worksheets("лаба6").Cells(i, j).Value = x Worksheets("лаба6").Cells(i, j + 1).Value = y hp = hp + 1 Loop Wend End Sub
Самостоятельная работа. Задание: Решить в программе Visual Basic методом Эйлера систему дифференциальных уравнений:
|
при начальных значениях х0=0; у0(х0)= -83; z0(x0)=17; на интервале [0; 10], шаг интегрирования h = 0,5.
