IPB
ЛогинПароль:

> Прочтите прежде чем задавать вопрос!

1. Заголовок темы должен быть информативным. В противном случае тема удаляется ...
2. Все тексты программ должны помещаться в теги [code=pas] ... [/code].
3. Прежде чем задавать вопрос, см. "FAQ", если там не нашли ответа, воспользуйтесь ПОИСКОМ, возможно такую задачу уже решали!
4. Не предлагайте свои решения на других языках, кроме Паскаля (исключение - только с согласия модератора).
5. НЕ используйте форум для личного общения, все что не относится к обсуждению темы - на PM!
6. Одна тема - один вопрос (задача)
7. Проверяйте программы перед тем, как разместить их на форуме!!!
8. Спрашивайте и отвечайте четко и по существу!!!

> Ввод и вывод матрицы, Простой вопрос
Shmaniche
сообщение 14.04.2013 13:19
Сообщение #1


Пионер
**

Группа: Пользователи
Сообщений: 79
Пол: Мужской

Репутация: -  0  +


Застопорился на простом – ввод и вывод матрицы – когда составлял одну программу. Хочу сделать так чтобы матрица красиво выводилась, но после ввода значений в вектор B, в нем выводятся одни нули. То есть матрица A выводиться правильно, а в векторе B одни нули. Где допустил ошибку не могу найти. Помогите разобраться. Вот код ввода и вывода:


program vv;

uses crt;
const c=10;
type mat=array [1..c, 1..c] of real;
vec=array [1..c] of real;
var A: mat;
B: vec;
n: integer;

{Процедура ввода матрицы}
procedure Matr_Vvod(n:integer; var A:mat; B:vec);
var i, j: integer;
begin
writeln('Вводите коэфф-ты матpицы A по стpокам нажимая ENTER:');
for i:=1 to n do
for j:=1 to n do
begin
write ('A[',i,',',j,']=');
read(A[i,j]);
end;
writeln;
writeln('Введите вектоp В свободных членов нажимая ENTER:');
for i:=1 to n do
begin
write ('B[',i,']=');
read(B[i]);
end;
end;

{Процедура вывода матрицы}
procedure Matr_Vivod(n: integer; var A:mat; B:vec);
var i, j: integer;
begin
writeln('Исходная матpица:');
for i:=1 to n do
begin
for j:=1 to n do
write (A[i, j]:5:3,' ':2); {Вывод матрицы с отступами}
write (B[i]:5:3,' ':2); {Вывод вектора с отступами}
writeln;
end;
end;

begin
clrscr;

writeln('Вывод матрицы.');
writeln;
writeln('Введите размерность матрицы и нажмите ENTER:');
readln(n);
writeln;

Matr_Vvod(n,A,B); {Ввод матрицы}
writeln;
Matr_Vivod(n,A,B); {Вывод матрицы}
writeln;

readkey;

end.


 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
 
 Ответить  Открыть новую тему 
Ответов
Shmaniche
сообщение 16.04.2013 14:01
Сообщение #2


Пионер
**

Группа: Пользователи
Сообщений: 79
Пол: Мужской

Репутация: -  0  +


Покопался вы коде и изменил всего одну строчку:
 
s:=A[i, n+1];


на:
 
s:=A[i, n] + B[i];


И после этого программа начала правильно находить два корня из трех. Третий корень выводит на единицу больше чем нужно. Например, имеем след. систему уравнений:
x1+4x2+2x3=3;
2x1+7x2+3x3=4;
8x2+x3=-5.
Итераций: 4.
Точность: 1.
Выводиться: x1=1; x2=-1; x3=4.
А должно быть: x1=1; x2=-1; x3=3.
И так практически во всех урав-ях. Третий «икс» выводиться на 1 больше. Помогите определить, что не так?
Полный код:
 
program Zeyd;
{Вычисление СЛАУ методом Зейделя}

uses crt;
const c=10;
type mat=array [1..c, 1..c] of real;
vec=array [1..c] of real;
var A: mat;
B, X: vec;
i,l,m,n: integer;
e: real;

{Процедура ввода матрицы}
procedure Matr_Vvod(n:integer; var A:mat; var B, X:vec);
var i, j: integer;
begin
writeln('Вводите коэфф-ты матpицы A по стpокам нажимая ENTER:');
for i:=1 to n do
for j:=1 to n do
begin
write ('A[',i,',',j,']=');
read(A[i,j]);
end;
writeln;
writeln('Введите вектоp В свободных членов нажимая ENTER:');
for i:=1 to n do
begin
write ('B[',i,']=');
read (B[i]);
end;
end;

{Процедура вывода матрицы}
procedure Matr_Vivod(n: integer; var A:mat; var B, X:vec);
var i, j: integer;
begin
writeln('Исходная матpица:');
for i:=1 to n do
begin
for j:=1 to n do
write (A[i, j]:5:3,' ':2); {Вывод матрицы с отступами}
write (B[i]:5:3,' ':2); {Вывод вектора с отступами}
writeln;
end;
end;

{Процедура вычисления СЛАУ методом Зейделя}
procedure Seid(var n,m,l:integer; var A:mat; var B, X:vec; var e:real);
var i,j,k,k1,n1:integer;
s:real;
begin
for k:=1 to m do
begin
l:=k;
for i:=1 to n do
begin
s:=A[i, n] +B[i];
for j:=1 to n do
s:=s - A[i,j] * X[j];
s:=s / A[i,i];
X[i]:=X[i] + s;
if abs(s) > e then l:=0;
end;
if l<>0 then exit;
end;
end;


begin
clrscr;

writeln('Вычислить СЛАУ методом Зейделя.');
writeln;
writeln('Введите колич-во неизвестных величин "n" и нажмите ENTER:');
readln(n);
writeln;
writeln('Введите максимально допустимое кол-во интераций "m":');
readln(m);
writeln;
writeln('Задайте погрешность решений "e":');
readln(e);
writeln;

Matr_Vvod(n,A,B,X); {Ввод матрицы}
writeln;
Matr_Vivod(n,A,B,X); {Вывод матрицы}
writeln;

Seid(n,m,l,A,B,X,e); {Вычисление СЛАУ методом Зейделя}

if l<>0 then
begin
writeln('ОТВЕТ:');
for i:=1 to n do
writeln('X', i,' = ',X[i]:3:4);
writeln;
writeln(l:2,' интерации(я).');
end
else
writeln('Интерации все.');

readkey;

end.



 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 

Сообщений в этой теме


 Ответить  Открыть новую тему 
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 



- Текстовая версия 23.06.2025 22:27
Хостинг предоставлен компанией "Веб Сервис Центр" при поддержке компании "ДокЛаб"