- •Теоретическая часть
- •Методы решения обыкновенных дифференциальных уравнений
- •Метод Эйлера
- •Метод Рунге-Кутта
- •Численное интегрирование систем обыкновенных дифференциальных уравнений
- •Хід виконання роботи
- •Задание
- •Отчет по лабораторной работе должен содержать:
- •Контрольные вопросы
- •Варианты заданий
- •Результат завдання
- •Висновки
Результат завдання
Метод Ейлера
Xe =
0 |
0.0100 |
0.0200 |
0.0300 |
0.0400 |
0.0500 |
0.0600 |
0.0700 |
0.0800 |
0.0900 |
0.1000 |
0.1100 |
0.1200 |
0.1300 |
0.1400 |
0.1500 |
0.1600 |
0.1700 |
0.1800 |
0.1900 |
0.2000 |
0.2100 |
0.2200 |
0.2300 |
0.2400 |
0.2500 |
0.2600 |
0.2700 |
0.2800 |
0.2900 |
0.3000 |
0.3100 |
0.3200 |
0.3300 |
0.3400 |
0.3500 |
0.3600 |
0.3700 |
0.3800 |
0.3900 |
0.4000 |
0.4100 |
0.4200 |
0.4300 |
0.4400 |
0.4500 |
0.4600 |
0.4700 |
0.4800 |
0.4900 |
0.5000 |
0.5100 |
0.5200 |
0.5300 |
0.5400 |
0.5500 |
0.5600 |
0.5700 |
0.5800 |
0.5900 |
0.6000 |
0.6100 |
0.6200 |
0.6300 |
0.6400 |
0.6500 |
0.6600 |
0.6700 |
0.6800 |
0.6900 |
0.7000 |
0.7100 |
0.7200 |
0.7300 |
0.7400 |
0.7500 |
0.7600 |
0.7700 |
0.7800 |
0.7900 |
0.8000 |
0.8100 |
0.8200 |
0.8300 |
0.8400 |
0.8500 |
0.8600 |
0.8700 |
0.8800 |
0.8900 |
0.9000 |
0.9100 |
0.9200 |
0.9300 |
0.9400 |
0.9500 |
0.9600 |
0.9700 |
0.9800 |
0.9900 |
1.0000 |
|
|
|
Ye =
1.0000 |
1.0000 |
1.0002 |
1.0006 |
1.0012 |
1.0020 |
1.0030 |
1.0042 |
1.0056 |
1.0072 |
1.0090 |
1.0110 |
1.0132 |
1.0155 |
1.0181 |
1.0209 |
1.0238 |
1.0270 |
1.0303 |
1.0339 |
1.0376 |
1.0415 |
1.0456 |
1.0499 |
1.0543 |
1.0590 |
1.0638 |
1.0688 |
1.0740 |
1.0793 |
1.0849 |
1.0906 |
1.0964 |
1.1025 |
1.1087 |
1.1151 |
1.1216 |
1.1283 |
1.1352 |
1.1422 |
1.1494 |
1.1568 |
1.1643 |
1.1720 |
1.1798 |
1.1878 |
1.1959 |
1.2042 |
1.2126 |
1.2212 |
1.2300 |
1.2388 |
1.2479 |
1.2571 |
1.2664 |
1.2758 |
1.2855 |
1.2952 |
1.3051 |
1.3151 |
1.3253 |
1.3357 |
1.3461 |
1.3567 |
1.3675 |
1.3784 |
1.3894 |
1.4006 |
1.4119 |
1.4233 |
1.4349 |
1.4467 |
1.4586 |
1.4706 |
1.4828 |
1.4951 |
1.5075 |
1.5201 |
1.5329 |
1.5458 |
1.5589 |
1.5721 |
1.5854 |
1.5990 |
1.6126 |
1.6265 |
1.6405 |
1.6546 |
1.6690 |
1.6835 |
1.6981 |
1.7130 |
1.7280 |
1.7432 |
1.7586 |
1.7741 |
1.7899 |
1.8059 |
1.8220 |
1.8384 |
1.8550 |
|
|
|
Метод Рунге-Кутта
a=0 %Початкове значення інтервалу диференціювання
b=1 %Кінцеве значення інтервалу диференціювання
ya=1 %Початкове значення аргумента
f=inline('tan(x)/y*(2)')
f(x,y) = tan(x)/y*(2)
[Xr,Yr]=ode45(f,a,b,ya)
plot(Xr,Yr) %Візуалізація функції
Xr = 0, 0.0250, 0.0500, 0.0750, 0.1000, 0.1250, 0.1500, 0.1750, 0.2000, 0.2250, 0.2500, 0.2750, 0.3000, 0.3250, 0.3500, 0.3750, 0.4000, 0.4250, 0.4500, 0.4750, 0.5000, 0.5250, 0.5500, 0.5750, 0.6000, 0.6250, 0.6500, 0.6750, 0.7000, 0.7250, 0.7500, 0.7750, 0.8000, 0.8250, 0.8500, 0.8750, 0.9000, 0.9250, 0.9500, 0.9750, 1.0000
