VPF::[Pascal] Решение СЛУ методом простых итераций

snook

Дата 21.12.2006, 19:57 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 21.12.2006

Репутация: нет
Всего: нет

Условие такое: написать программу решения систем линейных уравнений методом простых итераций и вывести результаты в файл results.txt

Нашел здесь как примерно делать: http://www.tspu.ru/ebooks/razina789/7.5.htm :

Код


program iter;
const n=4;
a: array[1..n,1..n] of real=((1.87,5.38,1.03,1.17),
(7.03,8.04,9.05,6.08),
(1.11,2.02,2.03,-.04),
(3.41,-4.52,7.28,5.18));
b: array [1..n] of real=(1.75,3.62,2.89,0.96);
label 1;
var i, j: integer;
m, x, y: array [1..n] of real;
p, e: real;
begin
for i:=1 to n do
begin
m[i]:=abs(b[i]);
for j:=1 to n do
m[i]:=m[i] + abs(a[i,j]);
end;
for i:=1 to n do
begin
if a[i,i]<0 then p:=1 else p:=-1;
for j:=1 to n do 
a[i,j]:=p*a[i,j]/m[i];
b[i]:=-p$*$b[i]/m[i];
a[i,i]:=a[i,i]+1;
writeln('b',i,'=',b[i]:4:4);
end;
for i:=1 to n do
begin
for j:=1 to n do
write(' ',a[i,j]:4:4);
writeln;
end;
writeln('e='); read(e);
for i:=1 to n do x[i]:=b[i];
1: for i:=1 to n do
begin
y[i]:=b[i];
for j:=1 to n do
y[i]:=y[i]+a[i,j]*x[j];
end; 
for i:=1 to n do
if abs(y[i]-x[i])>e then begin
for j:=1 to n do
x[j]:=y[j];
goto 1;
end;
for i:=1 to n do
writeln('x', i, '=', x[i]:3:3);
end.

В общем, этот код надо подстроить под моё условие. Данные должны вводиться с клавиатуры. Нужно этот метод, который используется в проге, сделать универсальным и объединить в функцию.

Спасибо за любые ответы. Завтра зачет. Надеюсь на вашу помощь, ибо больше надеяться не на кого.

M
alexeis1

Модератор: не забывайте устанавливать подсветку кода http://forum.vingrad.ru/index.php?showtopic=126445

Это сообщение отредактировал(а) alexeis1 - 21.12.2006, 21:57

snook

Дата 22.12.2006, 22:52 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 21.12.2006

Репутация: нет
Всего: нет

2Модератор

Удаляйте тему, если считаете нужным. Мне больше она не нужна и отписываться тут не буду.

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »