Метод Гаусса прямым ходом

@Gromjk · Регистрация: 22.12.2019

Студворк — интернет-сервис помощи студентам

Проблема с получением нижней треугольной матрицы. С верхней треугольной матрицей вроде нет проблем. Вот что получается в ответе:
U
1.0 5e-05 3.3333333333333335e-05 2.5e-05 2e-05 1.6666666666666667e-05
0.0 1.0 3.333000016665e-05 2.49975001249875e-05 1.9998000099990004e-05 1.6665000083325e-05
0.0 0.0 1.0 2.499625037495417e-05 1.9997000299963337e-05 1.6664166916636114e-05
0.0 0.0 0.0 1.0 1.9996333824926536e-05 1.6663611520772113e-05
0.0 0.0 0.0 0.0 1.0 1.6663194999904353e-05
0.0 0.0 0.0 0.0 0.0 1.0

L
9.99999985501611 0.0 0.0 0.0 2.710302164310409e-20 0.0
0.0019998100132212946 19.99999990500661 0.0 0.0 5.4208753466420714e-20 0.0
0.002999815007499694 0.001499907503749847 29.999999938335833 0.0 -2.032804537729561e-24 0.0
0.003999853335999951 0.0019999266679999756 0.001333284445333317 39.999999963334 1.0842021724855044e-19 0.0
0.004999916666666666 0.002499958333333333 0.001666638888888889 0.0012499791666666666 49.99999998333333 0.0
0.006 0.003 0.002 0.0015 0.0012 60.
код:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
import numpy as np
n=6
r=1
p=2
q=1
A=np.zeros((n,n))
LU=np.zeros((n,n))
#заполнение массива А
for i in range(n): 
    for j in range(n): 
        if i==j:
            A[i,j]=10*(i+1)**(p/2)
        else:
            A[i,j]=10**(-3)*((i+1)/((j+1)**r))**(1/q)   
B=np.zeros((n))
 #заполнение массива В
for i in range(n):
    B[i]=9*i**(p/2)
 
for i in range(len(A)): 
    for j in range(len(A[i])):
        print(A[i,j], end = ' ')
    print()
D=np.copy(A)
# прямой ход
for i in range(5):
    A[i,:]/=A[i,i]
    for j in range (i+1,6):
        A[j,:]-=A[i,:]*A[j,i]/A[i,i]
A[5,:]/=A[5,5]
U=A
for i in range(len(U)): 
    for j in range(len(U[i])):
        print(U[i,j], end = ' ')
    print()
# матрица L
for i in range(5,0,-1):
    for j in range (i-1,-1,-1):
        D[j,:]-=D[i,:]*D[j,i]/D[i,i]
L=D
for i in range(len(L)): 
    for j in range(len(L[i])):
        print(L[i,j], end = ' ')
    print()

@Gdez · 29.11.2020, 11:36

Gromjk, в чем именно проблема?

Добавлено через 1 час 11 минут
Посмотри тут, может поможет
https://orion1401.gitbooks.io/... henie.html

Новые блоги и статьи Все статьи Все блоги /
Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .	Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .

@Gdez 8840 / 4492 / 1864 Регистрация: 27.03.2020 Сообщений: 7,312
	29.11.2020, 11:36
	Gromjk, в чем именно проблема? Добавлено через 1 час 11 минут Посмотри тут, может поможет https://orion1401.gitbooks.io/... henie.html 0