MathCad 2000. Шейкер Т.Д. - 61 стр.

UptoLike

ВУЗ: 

ДВФУ | Владивосток

Составители: 

Рубрика: 

Иногда оказывается целесообразным попытаться решить
уравнение или систему несколькими методами, чтобы выбрать
более подходящий по точности или скорости получения результата.
12.2. Пример использования функции rkfixed
Решим рассмотренное ранее уравнение с помощью функции
rkfixed(y, x1, x2, n, D). Аналогичным образом для решения
уравнения можно было применить функцию Rkadapt или Bulstoer.
Последовательность действий:
задаём количество шагов
формируем вектор ic начальных условий. Так как
уравнение одно, то вектор содержит один элемент
задаём правые части D(t,Y). Вторым аргументом
должен быть вектор искомой функции
вызываем функцию rkfixed для получения матрицы
решения
выделяем первый столбец матрицы, который содержит
значения независимой переменной
выделяем второй столбец матрицы, который содержит
значения функции-результата
строим график
N : = 100 ic
0
:= y0
DtY,( ) sin t()
1
Y
0
+:=
S rkfixed ic t0, t1, N, D,():=
TS
0
〈〉
:=
YS
1
〈
〉
:=
0 5 10 15 20 25
5
10
Y
T
61
      Иногда оказывается целесообразным попытаться решить
уравнение или систему несколькими методами, чтобы выбрать
более подходящий по точности или скорости получения результата.

12.2. Пример использования функции rkfixed
       Решим рассмотренное ранее уравнение с помощью функции
rkfixed(y, x1, x2, n, D). Аналогичным образом для решения
уравнения можно было применить функцию Rkadapt или Bulstoer.
       Последовательность действий:
      ƒ задаём количество шагов
      ƒ формируем вектор ic начальных условий. Так как
          уравнение одно, то вектор содержит один элемент
      ƒ задаём правые части D(t,Y). Вторым аргументом
          должен быть вектор искомой функции
      ƒ вызываем функцию rkfixed для получения матрицы
          решения
      ƒ выделяем первый столбец матрицы, который содержит
          значения независимой переменной
      ƒ выделяем второй столбец матрицы, который содержит
          значения функции-результата
      ƒ строим график
      N : = 100         ic0 := y0
                              1
     D( t , Y) := sin( t) +
                              Y0
     S := rkfixed( ic , t0 , t1 , N , D)
     T := S 〈0〉           Y := S 〈1〉

                10



       Y          5




                      0       5     10         15   20   25
                                           T

                                    61

Страницы