Составители:
Рубрика:
31
Определим функцию, характеризующую вектор значений
производных в каждой точке решения (t, Y) (рис.1.2):
()
−−
=
10
1
2
5
2
1
2
1
:,
YY
Y
YtD
Зададим дополнительно:
t0:=0 – начальное значение аргумента;
t1:=10 – его конечное значение;
Y0:= – вектор начальных (например нулевых) значений
функций;
0
0
N = 100 – число решений на интервале [t0, t1].
Для решения можно использовать уже знакомую нам из
предыдущего примера функцию rkfixed, но она не единствен-
ная в арсенале встроенных функций Mathcad. Испытаем из них
фунцию
Rkadapt – также по методу Рунге - Кутта, но с адапти-
руемым шагом интегрирования, для чего напишем:
S:= Rkadapt(Y0, t0, t1, N, D)
t:=S
<0>
y0:=S
<1>
y1:=S
<2>
После этого можно ввести:
S =
и на экране высветится матрица значений
t, y0, y1 (всего по
умолчанию 15 строк). Затем строим график (рис. 1.3) действия-
ми, знакомыми из предыдущего примера.
Средства редактирования полученного графика аналогичны
упомянутым в примере 1.14.
Страницы
- « первая
- ‹ предыдущая
- …
- 29
- 30
- 31
- 32
- 33
- …
- следующая ›
- последняя »
