Решение системы трех линейных уравнений с тремя неизвестными методом Гаусса

@ТаНюШкЕ · Регистрация: 14.09.2014

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

Помогите пожалуйста составить вот такую программу
решения системы трех линейных уравнений с тремя неизвестными методом Гаусса.

@CodeR · 27.10.2014, 13:44

Ввод матрицы вручную или каким либо другим способом сами сделаете)

Pascal
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
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
Const nmax=10;
Type mas=array[1..nmax] of real;
     matr=array[1..nmax] of mas;
Var A:matr;
    b,x:mas;
    i,j,n:integer;
Procedure Gaus(A:matr;B:mas;Var X:mas);
Var i,j,k,imax:integer;
    temp,pr,sum:real;
Begin
    For k:=1 to n-1 do
    Begin
     imax:=k;
     For i:=k+1 to n do
      if abs(A[i][k])>abs(A[imax][k]) then imax:=i;
     if imax<>k then
     Begin
       For j:=1 to n do
       Begin
         temp:=A[k][j];
         A[k][j]:=A[imax][j];
         A[imax][j]:=temp;
       End;
       temp:=B[k];
       B[k]:=B[imax];
       B[imax]:=temp;
     End;
     For i:=k+1 to n do
     Begin
        pr:=a[i][k]/a[k][k];
        For j:=1 to n do
         a[i][j]:=a[i][j]-a[k][j]*pr;
        b[i]:=b[i]-b[k]*pr;
     End;
     End;
     x[n]:=b[n]/a[n][n];
     For k:=n-1 downto 1 do
     Begin
        sum:=0;
        For i:=k+1 to n do
         sum:=sum+x[i]*a[k][i];
        x[k]:=(b[k]-sum)/a[k][k];
     End;
End;
Begin
    Repeat
     write('Размерность матрицы 1<=N<=',nmax,': ');readln(n);
    Until n in [1..nmax];
    writeln('Матрица A');
    For i:=1 to n do
    Begin
      For j:=1 to n do
      Begin
        A[i][j]:=random(10)+1;
        write(A[i][j]:6:2);
        B[i]:=B[i]+A[i][j];
      End;
      writeln(' | ',B[i]:6:2);
    End;
    Gaus(A,B,X);
    write('X = (',X[1]:6:2);
    For i:=2 to n do
     write(';',X[i]:6:2);
    writeln('  )');
    readln;
End.

Добавлено через 11 минут
Вот с ручным вводом и проверкой на вырожденность

Pascal
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
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
Const nmax=10;
Type mas=array[1..nmax] of real;
     matr=array[1..nmax] of mas;
Var A:matr;
    b,x:mas;
    i,j,n:integer;
    flag:boolean;
Procedure Gaus(A:matr;B:mas;Var X:mas;Var flag:boolean);
Var i,j,k,imax:integer;
    temp,pr,sum:real;
Begin
    flag:=true;
    For k:=1 to n-1 do
    Begin
     imax:=k;
     For i:=k+1 to n do
      if abs(A[i][k])>abs(A[imax][k]) then imax:=i;
     if imax<>k then
     Begin
       For j:=1 to n do
       Begin
         temp:=A[k][j];
         A[k][j]:=A[imax][j];
         A[imax][j]:=temp;
       End;
       temp:=B[k];
       B[k]:=B[imax];
       B[imax]:=temp;
     End;
     For i:=k+1 to n do
     Begin
        pr:=a[i][k]/a[k][k];
        For j:=1 to n do
         a[i][j]:=a[i][j]-a[k][j]*pr;
        b[i]:=b[i]-b[k]*pr;
        if a[i][i]=0 then
        Begin
          flag:=false;
          Exit;
        End;
     End;
     End;
     x[n]:=b[n]/a[n][n];
     For k:=n-1 downto 1 do
     Begin
        sum:=0;
        For i:=k+1 to n do
         sum:=sum+x[i]*a[k][i];
        x[k]:=(b[k]-sum)/a[k][k];
     End;
End;
Begin
    Repeat
     write('Введите размерность 1<=N<=',nmax,': ');readln(n);
    Until n in [1..nmax];
    For i:=1 to n do
     For j:=1 to n do
     Begin
       write('A[',i,',',j,'] = ');
       readln(A[i][j]);
     End;
    For i:=1 to n do
    Begin
       write('B[',i,'] = ');
       readln(B[i]);
    End;
    writeln('Ìàòðèöà A');
    For i:=1 to n do
    Begin
      For j:=1 to n do
       write(A[i][j]:6:2);
      writeln(' | ',B[i]:6:2);
    End;
    Gaus(A,B,X,flag);
    if flag then
    Begin
      write('X = (',X[1]:6:2);
      For i:=2 to n do
       write(';',X[i]:6:2);
      writeln('  )');
    End
    else writeln('Система вырождена');
    readln;
End.

Новые блоги и статьи Все статьи Все блоги /
Трассировка корутин Kotlin с OpenTelemetry mobDevWorks 14.07.2025 Асинхронное программирование меняет правила игры, особенно когда речь заходит о трассировке операций. В Kotlin с его корутинами эта проблема приобретает особый оттенок, который я хотел бы детально. . .	Облачные приложения на Rust: руководство по архитектуре микросервисов golander 13.07.2025 Когда я впервые взялся за проектирование облачной платформы для одного из наших клиентов, выбор стоял между привычными Go и Java. Но после нескольких месяцев разработки микросервисной системы,. . .	Как Node.js выполняет асинхронные операции Reangularity 13.07.2025 Каждый раз, когда я рассказываю про Node. js, возникает один и тот же вопрос: "Как эта штука может быть быстрой, если JavaScript — однопоточный язык?" И это действительно кажется парадоксом. Ведь в. . .	Как писать чистый, тестируемый и качественный код на Python py-thonny 12.07.2025 Помню свой первый проект на Python. Работал тогда быстро, грязно, лишь бы работало. Код был похож на запутанный клубок - переменные по одной букве, функции на 200 строк, комментарии отсутствовали как. . .	Blazor и контроллер сервопривода IoT Meadow Maple Wired 11.07.2025 Я решил разобраться, как можно соединить современные веб-технологии с миром "железа". Интересная комбинация получилась из Blazor в качестве веб-интерфейса и микроконтроллера Meadow с его веб-сервером. . .
Генерация OpenQASM из кода Q# EggHead 10.07.2025 Летом 2024-го я начал эксперименты с библиотекой Q# Bridge, и знаете что? Она оказалась просто находкой для тех, кто работает на стыке разных квантовых экосистем. Основная фишка этой библиотеки -. . .	Изучаем новый шаблон ИИ-чата .NET AI Chat Web App stackOverflow 10.07.2025 В . NET появилось интересное обновление - новый шаблон ИИ-чата под названием . NET AI Chat Web App. Когда я впервые наткнулся на анонс этого шаблона, то сразу понял, что Microsoft наконец-то. . .	Результаты исследования от команды ARP (июль 2025 г.) Programma_Boinc 10.07.2025 Результаты исследования от команды ARP (июль 2025 г. ) Африканский проект по дождям (ARP) World Community Grid снова запущен! Мы рады поделиться обновленной информацией о нашем прогрессе с осени. . .	Angular vs Svelte - что лучше? Reangularity 09.07.2025 Сегодня рынок разделился на несколько четких категорий: тяжеловесы корпоративного уровня (Angular), гибкие универсалы (React), прогрессивные решения (Vue) и новая волна компилируемых фреймворков. . .	Code First и Database First в Entity Framework UnmanagedCoder 09.07.2025 Entity Framework дает нам свободу выбора, предлагая как Code First, так и Database First подходы. Но эта свобода порождает вечный вопрос — какой подход выбрать? Entity Framework — это. . .

@ТаНюШкЕ 1 / 1 / 0 Регистрация: 14.09.2014 Сообщений: 72

	Решение системы трех линейных уравнений с тремя неизвестными методом Гаусса 26.10.2014, 16:28. Показов 5002. Ответов 1 Метки нет (Все метки) Помогите пожалуйста составить вот такую программу решения системы трех линейных уравнений с тремя неизвестными методом Гаусса. 0

Опции темы