Решение уравнения методом Гаусса

@tephnut10 · Регистрация: 09.03.2010

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

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
program gauss;
uses crt;
var s:real; i,j,k,n:integer;
x: array [1..10] of real; a: array [1..10,1..11] of real;
begin
clrscr;
write ('укажите число уравнений n=');
readln (n);
writeln ('задайте значения коэффициентов системы');
for i:=1 to n do
for j:=1 to n do
begin
write ('A(',i:2,',',j:2,')=');
readln (a[i,j]);
end;   writeln;
for k:=1 to n-1 do
begin
for j:=k+1 to n+1 do
a[i,j]:=a[k,j]/a[k,k];
a[k,k]:=1;
for i:=k+1 to n do
begin
for j:=k+1 to n+1 do
a[i,j]:=a[i,j]-a[k,j]*a[i,k];
a[i,k]:=0;
end;
end;
if abs(a[n,n])> 1e-50 then begin
x[n]:=a[n,n+1]/a[n,n];
for i:=n-1 downto 1 do
begin
s:=0;
for j:=i+1 to n do s:=s+a[i,j]*x[j];
x[i]:=a[i,n+1]-s;
end;
writeln ('решение системы');
for i:=1 to n do
begin
writeln ('x(',i:2,')=',x[i]:2:5);
end;
end
else
writeln ('система решения не имеет ');
end.

корни все равны 0 а математически решив не так
найти ошибку заранее спасибо

))))))))))))))))

@Puporev · 29.10.2010, 22:15

tephnut10, Вообще-то у систем кроме коэффициентов есть еще и свободные члены, которые ты в программу не вводишь...

@tephnut10 · 30.10.2010, 17:21 **[ТС]**

все равно корни 0

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
program gauss;
uses crt;
var s:real; i,j,k,n:integer;
x: array [1..10] of real; a: array [1..10,1..11] of real;
begin
clrscr;
write ('укажите число уравнений n=');
readln (n);
writeln ('задайте значения коэффициентов системы');
for i:=1 to n do
for j:=1 to n +1 do
begin
write ('A(',i:2,',',j:2,')=');
readln (a[i,j]);
end;   writeln;
for k:=1 to n-1 do
begin
for j:=k+1 to n+1 do
a[i,j]:=a[k,j]/a[k,k];
a[k,k]:=1;
for i:=k+1 to n do
begin
for j:=k+1 to n+1 do
a[i,j]:=a[i,j]-a[k,j]*a[i,k];
a[i,k]:=0;
end;
end;
if abs(a[n,n])> 1e-50 then begin
x[n]:=a[n,n+1]/a[n,n];
for i:=n-1 downto 1 do
begin
s:=0;
for j:=i+1 to n do s:=s+a[i,j]*x[j];
x[i]:=a[i,n+1]-s;
end;
writeln ('решение системы');
for i:=1 to n do
begin
writeln ('x(',i:2,')=',x[i]:2:5);
end;
end;
else
writeln ('система решения не имеет ');
end.

@Puporev · 30.10.2010, 17:56

tephnut10, Ты неправильно обрабатываешь свободные члены вместе с коэффициентами в одной матрице, их нужно отдельно.
И не все корни у тебя ноли, а один, хотя все равно неверно считает.
Вот посмотри мой код, как-то писал, добавь только точность.

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
uses crt;
var a:array[1..20,1..20] of real;
    b:array[1..20] of real;
    x:array[1..20] of real;
    n,i,j,k:integer;z,r,g:real;
begin
clrscr;
write('Введите количество уравнений n=');
readln(n);
writeln('Введите коэффициенты системы и свободные члены');
for i:=1 to n do
for j:=1 to n+1 do
read(A[i,j]);
readln;
clrscr;
writeln('Матрица коэффициентов и свободных членов');
for i:=1 to n do
  begin
    for j:=1 to n+1 do
    write(A[i,j]:6:2);
    writeln;
  end;
readln;
//прямой ход
for i:=1 to n do
b[i]:=a[i,n+1];   //свободные члены в отдельный массив
for k:=1 to n do
 begin
  for j:=k+1 to n do
    begin
     r:=a[j,k]/a[k,k]; //коэффициент пересчета
     for i:=k to n do
     a[j,i]:=a[j,i]-r*a[k,i]; //пересчитываем коэффициенты
     b[j]:=b[j]-r*b[k]; //свободные члены
    end;
 end;
 //обратный ход
for k:=n downto 1 do
 begin
  r:=0;
  for j:=k+1 to n do
   begin
    g:=a[k,j]*x[j];
    r:=r+g;
   end;
  x[k]:=(b[k]-r)/a[k,k];
 end;
writeln('Корни системы:');
for i:=1 to n do
write('x[',i,']=',x[i]:0:2,'   ');
readln;
end.

Новые блоги и статьи Все статьи Все блоги /
CQRS (Command Query Responsibility Segregation) на Java Javaican 12.05.2025 CQRS — Command Query Responsibility Segregation, или разделение ответственности команд и запросов. Суть этого архитектурного паттерна проста: операции чтения данных (запросы) отделяются от операций. . .	Шаблоны и приёмы реализации DDD на C# stackOverflow 12.05.2025 Когда я впервые погрузился в мир Domain-Driven Design, мне показалось, что это очередная модная методология, которая скоро канет в лету. Однако годы практики убедили меня в обратном. DDD — не просто. . .	Исследование рантаймов контейнеров Docker, containerd и rkt Mr. Docker 11.05.2025 Когда мы говорим о контейнерных рантаймах, мы обсуждаем программные компоненты, отвечающие за исполнение контейнеризованных приложений. Это тот слой, который берет образ контейнера и превращает его в. . .	Micronaut и GraalVM - будущее микросервисов на Java? Javaican 11.05.2025 Облачные вычисления безжалостно обнажили ахиллесову пяту Java — прожорливость к ресурсам и медлительный старт приложений. Традиционные фреймворки, годами радовавшие корпоративных разработчиков своей. . .	Инфраструктура как код на C# stackOverflow 11.05.2025 IaC — это управление и развертывание инфраструктуры через машиночитаемые файлы определений, а не через физическую настройку оборудования или интерактивные инструменты. Представьте: все ваши серверы,. . .
Инъекция зависимостей в ASP.NET Core - Практический подход UnmanagedCoder 11.05.2025 Инъекция зависимостей (Dependency Injection, DI) — это техника программирования, которая кардинально меняет подход к управлению зависимостями в приложениях. Представьте модульный дом, где каждая. . .	Битва за скорость: может ли Java догнать Rust и C++? Javaican 11.05.2025 Java, с её мантрой "напиши один раз, запускай где угодно", десятилетиями остаётся в тени своих "быстрых" собратьев, когда речь заходит о сырой вычислительной мощи. Rust и C++ традиционно занимают. . .	Упрощение разработки облачной инфраструктуры с Golang golander 11.05.2025 Причины популярности Go в облачной инфраструктуре просты и одновременно глубоки. Прежде всего — поразительная конкурентность, реализованная через горутины, которые дешевле традиционных потоков в. . .	Создание конвейеров данных ETL с помощью Pandas AI_Generated 10.05.2025 Помню свой первый опыт работы с большим датасетом — это была катастрофа из неотформатированных CSV-файлов, странных значений NULL и дубликатов, от которых ехала крыша. Тогда я потратил три дня на. . .	C++ и OpenCV - Гайд по продвинутому компьютерному зрению bytestream 10.05.2025 Компьютерное зрение — одна из тех технологий, которые буквально меняют мир на наших глазах. Если оглянуться на несколько лет назад, то сложно представить, что алгоритмы смогут не просто распознавать. . .

@Puporev Почетный модератор 64312 / 47609 / 32742 Регистрация: 18.05.2008 Сообщений: 115,181
	29.10.2010, 22:15
	tephnut10, Вообще-то у систем кроме коэффициентов есть еще и свободные члены, которые ты в программу не вводишь... 0

Опции темы