Решение дифференциальных уравнений в системе компьютерной математики Maxima. Губина Т.Н - 92 стр.

UptoLike

Т.Н. Губина, Е.В. Андропова
Теперь заполним матрицу а по формулам:
a
k , k
=h
2
q
k
h p
k
1, a
k , k1
=h p
k
2, a
k ,k 2
=1, k =1, n1
,
a
n ,1
=
1
h
2
, a
n,2
=
2
, a
n1, n
=
2
, a
n 1, n1
=
1
h
2
, k =n , n1
Для заполнения коэффициентами первых двух уравнений системы
воспользуемся циклом с параметром:
В двух последних уравнениях поменяем значения некоторых элементов
с помощью оператора присваивания:
Теперь заполним столбец свободных членов:
Выведем полученные матрицы на экран:
Получили систему линейных уравнений, записанную в матричном виде
ay=b
, где
y
— искомое решение. Найдем его матричным способом.
92
                                                                    Т.Н. Губина, Е.В. Андропова




     Теперь заполним матрицу а по формулам:
                         2
                   a k , k =h q k −h p k 1, a k , k1=h p k −2, a k ,k 2=1, k =1, n−1 ,
        a n ,1 =1 h− 2 , a n ,2=2 , a n1, n=−2 , a n 1, n1=1 h2 , k =n , n1
     Для заполнения коэффициентами первых двух уравнений системы
воспользуемся циклом с параметром:




     В двух последних уравнениях поменяем значения некоторых элементов
с помощью оператора присваивания:


     Теперь заполним столбец свободных членов:



     Выведем полученные матрицы на экран:




     Получили систему линейных уравнений, записанную в матричном виде
 ay=b , где y — искомое решение. Найдем его матричным способом.




                                              92