Учебная работа № 341475. Тема: Вычислительная математика

Учебная работа № 341475. Тема: Вычислительная математика

Выдержка из похожей работы

…

Лабораторные работы по вычислительной математике

…….xn,
xn-1,
… x1.
Он
начинается с равенства: xn=bn+1;
const
max=10;
type
matrix=array[1..max]
of
real;
matrix_2=array[0..max]
of
real;
procedure
input_matr(var
a:matrix;n:byte;c:char);
var
i:byte;
begin
for
i:=1 to
n do
begin
writeln(‘Введите
‘,i ,’-ый элемент массива ‘,c);
readln(a[i])
end
end;
procedure
process(a,b,c,d:matrix;var
x:matrix;n:byte);
var
alfa,betta:matrix_2;
gamma,fi:matrix;
i:byte;
begin
betta[0]:=0;
alfa[0]:=0;
for
i:=1 to
n do
begin
gamma[i]:=b[i]+a[i]*alfa[i-1];
fi[i]:=d[i]-a[i]*betta[i-1];
alfa[i]:=-c[i]/gamma[i];
betta[i]:=fi[i]/gamma[i]
end;
x[n]:=betta[n];
for
i:=n-1 downto
1 do
x[i]:=alfa[i]*x[i+1]+betta[i]
end;
procedure
out_matr_x(a:matrix;n:byte);
var
i:byte;
begin
for
i:=1 to
n do
writeln(i ,’ корень уравнения равен ‘,a[i]:5:3)
end;
var
i:byte;
a,b,c,d,x,gamma,fi:matrix;
alfa,betta:matrix_2;
n:byte;
begin
writeln(‘Введите
размерность системы ‘);
readln(n);
if
(n>=2) and
(nthen
begin
input_matr(a,n,’a’);
input_matr(b,n,’b’);
input_matr(c,n,’c’);
input_matr(d,n,’d’);
process(a,b,c,d,x,n);
out_matr_x(x,n)
end
else
writeln(‘1
end.
Результат
работы программы:
1
корень уравнения равен 0.526
2
корень уравнения равен 0.628
3
корень уравнения равен 0.640
4
корень уравнения равен 1.200

…