4.Программа:
var
x10,x20,x30,x40,x11,x21,x31,x41,e,s:real;j:integer;
begin
writeln
(' vv nach priblig');
read(x10,x20,x30,x40);
x11:=x10;x21:=x20;x31:=x30;x41:=x40;j:=0;
repeat
x10:=x11;x20:=x21;x30:=x31;x40:=x41;
x11:=(125-27*x40-22*x30-x20*11)/(65);
x21:=(-113-19*x40+30*x30+7*x11)/65;
x31:=(139+32*x40-34*x21-21*x11)/92;
x41:=(110-35*x31-3*x21+22*x11)/65;
s:=abs(x11-x10)+abs(x21-x20)+abs(x31-x30)+abs(x41-x40);
inc(j)
until
s<0.00001;
write(x11:5:2,'
',x21:5:2,' ',x31:5:2,' ',x41:5:2,' nomer itter ',j);
end.
5.Результаты:
Turbo
Pascal Version 7.0 Copyright (c) 1983,92 Borland International
vv
nach priblig
456
654
777
6
1.00
-1.00 2.00 1.00
nomer itter 26
6.Результаты
вычисления функции в maxima:
(%i4)
solve([-113-19*s+30*z-65*y+7*x,125-27*s-22*z-11*y-65*x,139+32*s-92*z-34*y-21*x,110
-65*s-35*z-3*y+22*x], [x,y,z,s]);
(%o4)
[[x=1,y=-1,z=2,s=1]]
(%i5)
7.Вывод:
Вычисление
системы (*) в maxima и в Паскале приблизительно
равны. Система (*) решена с помощью метода
Зейделя
и получены истинные значения системы(*).
| 
