Файл: Эксперименты лаба8,9(2курс).pdf

f:=false;
for z from 1 to n do

if abs(x1[z]-x[z])>eps then f:=true end if;

end do;

end do;
for i from 1 to n do

print(x1[i]);

end do;
end proc;

# ввод коэффициентов матрицы

A

, вектора

B

и начального произ-

вольного приближения к решению системы
> a[1,1]:=23.81837: a[1,2]:=2.33377: a[1,3]:=-0.73872: a[1,4]:=0.85797:
a[1,5]:=-0.67457:
> a[2,1]:=-1.59171: a[2,2]:=29.48598: a[2,3]:=1.81520: a[2,4]:=3.90203:
a[2,5]:=-0.57008:
> a[3,1]:=-0.15720: a[3,2]:=1.44805: a[3,3]:=21.33172: a[3,4]:=2.64148:
a[3,5]:=3.36273:
> a[4,1]:=-3.29147: a[4,2]:=2.64912: a[4,3]:=2.16798: a[4,4]:=26.72888:
a[4,5]:=1.21258:
> a[5,1]:=-1.55251: a[5,2]:=2.73236: a[5,3]:=2.00201: a[5,4]:=-1.97732:
a[5,5]:=29.53347:
> b[1]:=1.77836: b[2]:=0.10234: b[3]:=3.82426: b[4]:=-2.71881: b[5]:=0.72844:
> x[1]:=1: x[2]:=1: x[3]:=1: x[4]:=1: x[5]:=1:

Проверим работу процедуры
> Zeidel(a,b,x,5,0.01);

0.08365479557
0.01070750424

0.1914082408

-0.1083504947

0.007842401791

16