Mathcad : математический практикум. Часть 2. Есипенко Д.Г - 8 стр.

UptoLike

Рубрика: 

8
ORIGIN 1:= y
1
1:= Dxy,()sinxy
1
(
)
:=
Yrkfixedy 0, π, 20, D,
()
:=
01234
1
2
Y
2
〈〉
Y
1
〈〉
Y
1 2
1
2
3
4
5
6
7
8
9
10
01
0.1571.012
0.3141.05
0.4711.115
0.6281.208
0.7851.331
0.9421.477
1.11.633
1.2571.774
1.4141.874
=
Указание. Присвойте переменной ORIGIN значение, равное 1, чтобы
нумерация компонент вектора начиналась с единицы . При решении задачи
Коши для уравнения первого порядка вектор решения имеет единственную
компоненту
1
y
. Вектор правых частей также содержит одну компоненту.
Перед обращением к функции rkfixed присвойте переменной
1
y
начальное
значение, равное единице, а переменной D(х,у) - выражение для правой части
уравнения , равное
)sin(
1
xy
. В приведенном фрагменте результаты вычислений
функции rkfixed(y, 0, π ,20, D) присвоены матрице Y; она содержит в первом
столице координаты 20 узлов равномерной сетки на отрезке
]
,
0
[
π
, а во
втором - приближенные значения решения в этих узлах (в приведенном
документе представлена только часть таблицы ). Для того чтобы вывести в
рабочий документ таблицу значений решения , нажмите на клавиатуре
клавиши <Y><=>, и таблица значений будет выведена в рабочем документе
справа от знака равенства . Щелкните по полю таблицы и просмотрите ее,
используя стрелки прокрутки. Для того чтобы построить график решения,
щелкните по кнопке в панели графиков и введите в качестве
переменной на оси абсцисс
>< 1
Y
(столбец координат узлов сетки), а на оси
ординат -
>< 2
Y
(столбец значений решения в узлах сетки). Чтобы ввести
номер столбца, введите имя матрицы , щелкните по кнопке в панели
и введите в помеченной позиции номер столбца.
                                                     8
      ORIGIN := 1         y1 := 1          D ( x , y) := sin(x ⋅ y1)             1       2
      Y := rkfixed(y , 0 , π , 20 , D)                                      1        0       1
                                                                            2   0.157 1.012
                            2                                               3   0.314    1.05
                                                                            4   0.471 1.115
                                                                       Y=   5   0.628 1.208
                   Y ��
                      2     1                                               6   0.785 1.331
                                                                            7   0.942 1.477
                                                                            8     1.1 1.633
                                0    1     2     3       4                  9   1.257 1.774

                                         Y ��
                                            1                               10 1.414 1.874



     Указание. Присвойте переменной ORIGIN значение, равное 1, чтобы
нумерация компонент вектора начиналась с единицы. При решении задачи
Коши для уравнения первого порядка вектор решения имеет единственную
компоненту y 1 . Вектор правых частей также содержит одну компоненту.
Перед обращением к функции rkfixed присвойте переменной y 1 начальное
значение, равное единице, а переменной D(х,у) - выражение для правой части
уравнения, равное sin(xy1). В приведенном фрагменте результаты вычислений
функции rkfixed(y, 0, π ,20, D) присвоены матрице Y; она содержит в первом
столице координаты 20 узлов равномерной сетки на отрезке [0, π ] , а во
втором - приближенные значения решения в этих узлах (в приведенном
документе представлена только часть таблицы). Для того чтобы вывести в
рабочий документ таблицу значений решения, нажмите на клавиатуре
клавиши <=>, и таблица значений будет выведена в рабочем документе
справа от знака равенства. Щелкните по полю таблицы и просмотрите ее,
используя стрелки прокрутки. Для того чтобы построить график решения,

щелкните по кнопке                  в панели графиков                  и введите в качестве
                                     <1>
переменной на оси абсцисс Y                (столбец координат узлов сетки), а на оси
            <2 >
ординат - Y      (столбец значений решения в узлах сетки). Чтобы ввести

номер столбца, введите имя матрицы, щелкните по кнопке                           в панели
и введите в помеченной позиции номер столбца.