/p>
В
Файл-функція:
function [f3, f4, f5] = mad (x) = 2 * sin (3 * x) +3 * cos (2 * x); = (sqrt (x +1) +2 * x)/((x ^ 2) +1);
f5 = -4 * (x ^ (-2));
Файл-програма:
x = 7;
[f3, f4, f5] = mad (x)
Результат роботи програми:
f3 = 2.0835 = 0.3366 = -0.0816
Рішення диференціальних рівнянь
,
а) з використанням класичних алгоритмів чисельних методів Ейлера і Рунге-Кутта4-го порядку; Метод Ейлера:
Алгоритм рішення диференціального рівняння методом Ейлера
В
Файл-функція:
function f = fun (x) = x;
Файл-програма:
t0 = 0; = 5; = 0.05; = 5; = 10; = 0; ('Таблиця значень за МЕ n'); ('.......... ................ n '); (' | t | x | y | n ');
fprintf ('.......................... n'); t = t0: h: tk
x = (a/b) * sqrt (t);
f = faaan (x);
y = y + h * f;
fprintf ('|% 4.2f |% 4.4f |% 4.4f | n', t, x, y);
plot (t, x, 'r +')
grid on
hold on
plot (t, y, 'b *')
xlabel ('t');
ylabel ('y');
Результат роботи програми:
Таблиця значень за МЕ
..........................
| t | x | y _ |
.................................... ..
| 0.00 | 0.0000 | 0.0000 |
| 0.01 | 0.0500 | 0.0005 |
| 0.02 | 0.0707 | 0.0012 |
| 0.03 | 0.0866 | 0.0021 |
| 0.04 | 0.1000 | 0.0031 |
| 0.05 | 0.1118 | 0.0042 |
| 0.06 | 0.1225 | 0.0054 |
| 0.07 | 0.1323 | 0.0067 |
| 0.08 | 0.1414 | 0.0082 |
| 0.09 | 0.1500 | 0.0097 |
| 0.10 | 0.1581 | 0.0112 |
| 0.11 | 0.1658 | 0.0129 |
