Реферат Рішення алгебраїчних і диференціальних рівнянь у пакеті MatLab

> | 4.83 | 1.0989 | 3.5437 |

| 4.84 | 1.1000 | 3.5547 |

| 4.85 | 1.1011 | 3.5657 |

| 4.86 | 1.1023 | 3.5768 |

| 4.87 | 1.1034 | 3.5878 |

| 4.88 | 1.1045 | 3.5988 |

| 4.89 | 1.1057 | 3.6099 |

| 4.90 | 1.1068 | 3.6210 |

| 4.91 | 1.1079 | 3.6320 |

| 4.92 | 1.1091 | 3.6431 |

| 4.93 | 1.1102 | 3.6542 |

| 4.94 | 1.1113 | 3.6654 |

| 4.95 | 1.1124 | 3.6765 |

| 4.96 | 1.1136 | 3.6876 |

| 4.97 | 1.1147 | 3.6988 |

| 4.98 | 1.1158 | 3.7099 |

| 4.99 | 1.1169 | 3.7211 |

| 5.00 | 1.1180 | 3.7323 |

В 

Метод Рунге-Кутта 4-го порядку:

Алгоритм рішення диференціального рівняння

методом Рунге-Кутта 4-го порядку.

В 

Файл-функція:

f = faaan (x) = x;

Файл-програма:

a = 5; = 10; = 0.05; = 0; = 5; = 0; = 0; ('Таблиця значень за МРГ n'); ('------- ------------------- n '); (' | t | x | y | n ');

fprintf ('-------------------------- n'); = t0; = (a/b ) * sqrt (t); ('|% 4.2f |% 4.4f |% 4.4f | n', t, x, y); t

f = faaan (x);

k = h * f;

s = k;

y = y + k/2;

t = t + h/2;

x = (a/b) * sqrt (t);

for j = 1:2

f = faaan (x);

k = h * f;

s = s +2 * k;

y = y1 + k/2;

end

t = t + h/2;

x = (a/b) * sqrt (t);

y = y + k/2;

f = faaan (x);

k = h * f;

s = s + k;

y = y1 + s/6;

fprintf ('|% 4.2f |% 4.4f |% 4.4f | n', t, x, y);

plot (t, x, 'r *')

grid on

hold on

plot (t, y, 'k +')

x...