Форум программистов, компьютерный форум, киберфорум
Наши страницы
Turbo Pascal
Войти
Регистрация
Восстановить пароль
 
David_Villa7
0 / 0 / 0
Регистрация: 24.11.2010
Сообщений: 20
1

Машинный ноль

22.12.2010, 12:22. Просмотров 1327. Ответов 0
Метки нет (Все метки)

Что такое машинный ноль, если не сложно, можно как нибудь попроще, а то как в википедии написано ваще понять не могу и какую роль он играет в этом алгоритме? Помогите пожалуйста, заранее спасибо


Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
const MAX = 10;
type Tmatrix = array [1..MAX, 1..MAX] of real;
     Tmassiv = array [1..MAX] of real;
     Tmassiv2 = array[1..MAX] of integer;
 
var a: Tmatrix;
    b, x, a0: Tmassiv;
    int: Tmassiv2;
    n, maxiter, retcode, i,j,k: integer;
    eps, w: real;
 
procedure NormMatr(n:integer; var a:Tmatrix; var A0:Tmassiv);
begin
   for i:=1 to n do
   begin
      a0[i]:=a[1,i];
      for j:=2 to n do
         if abs(a[j,i])>a0[i] then
            a0[i]:=abs(a[j,i]);
      for j:=1 to n do
         a[j,i]:=a[j,i]/a0[i];
   end;
end;
 
procedure KrautDeComp(eps:real; n:integer; var a:Tmatrix; var int:Tmassiv2; var A0:Tmassiv; var retcode:integer);
var imax: integer;
    amax: real;
begin
   NormMatr(n,a,a0);
   for i:=1 to n do
      int[i]:=i;
   for k:=1 to n-1 do
   begin
      amax:=abs(a[int[k], k]);
      for i:=k+1 to n do
         if abs(a[int[k], i])>amax then
         begin
            amax:=abs(a[int[k], i]);
            imax:=i;
         end;
      if imax<>k then
      begin
         j:=int[k];
         int[k]:=int[imax];
         int[imax]:=j;
      end;
      if abs(a[int[k], k])<n*eps then
      begin
         retcode:=1;
         break;
      end;
      for j:=k+1 to n do
         a[int[k], j]:=a[int[k], j]/a[int[k], k];
      for i:=k+1 to n do
         for j:=k+1 to n do
            a[int[i], j]:=a[int[i], j]-a[int[k], j]*a[int[i], k];
   end;
end;
 
procedure KrautSLU(n:integer; var b:Tmassiv; var a:Tmatrix;var int:Tmassiv2;var A0:Tmassiv; var x:Tmassiv);
var sum, z: real;
begin
   for i:=1 to n do
      x[i]:=b[i];
   for i:=1 to n do
      if int[i]<>i then
      begin
         z:=x[i];
         x[i]:=x[int[i]];
         x[int[i]]:=z;
      end;
   for i:=1 to n do
   begin
      sum:=0;
      for j:=1 to i-1 do
         sum:=sum+a[i,j]*x[j];
      x[i]:=(x[i]-sum)/a[i,i];
   end;
   for i:=n downto 1 do
   begin
      sum:=0;
      for j:=i+1 to n do
         sum:=sum+A[i,j]*x[j];
      x[i]:=x[i]-sum;
   end;
   for i:=1 to n do
      x[i]:=x[i]/a0[i];
end;
 
function KrautSolve(eps:real; n:integer; var a:Tmatrix; var b, x:Tmassiv;var int:Tmassiv2;var A0: Tmassiv; var retcode:integer): integer;
begin
   retcode:=0;
   KrautDeComp(eps,n,A,int,a0,retcode);
   if retcode = 0 then
      KrautSLU(n,b,a,int,A0,x);
end;
 
procedure Vvod(var A: Tmatrix; var B: Tmassiv; var n: integer);
begin
   writeln('‚введите количество неизвестных:');
   readln(n);
   writeln('введите матрицу коээффициентов:');
   for i:=1 to n do
      for j:=1 to n do
         read(a[i, j]);
   writeln('введите правую часть:');
   for i:=1 to n do
      read(b[i]);
end;
procedure Print(var x:Tmassiv; const n: integer; const retcode: integer);
begin
   if retcode = 1 then writeln('нет решения') else
   begin
      writeln('решение:');
      for i:=1 to n do
         writeln('x[',i,']= ',x[i]:0:4);
   end;
end;
 
begin
   Vvod(a,b,n);
   eps:=0.001;
   KrautSolve(eps,n,a,b,x,int,a0,retcode);
   Print(x,n,retcode);
   readln;
   readln;
end.

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

0
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
22.12.2010, 12:22
Ответы с готовыми решениями:

Машинный ноль
program Epsilon; var epsilon: Real; begin epsilon := 1; while...

есть ли ноль?
помогите пожалуйста решить задачу: определить есть ли среди первых трех из...

Деление на ноль
У прогу пишу,для проверки. Вроде нормально, а у меня деление на ноль пишет. ...

0 (ноль) в программировании
Бывают задания посчитать столько положительных элементов в массиве и прочее......

Ошибка деление на ноль
в конкретном случае (a&gt;1/a) при а= 0 естественно ошибка деление на ноль, как...

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
22.12.2010, 12:22

Вещественное деление на ноль
Пишет вещественное деление на 0. Помогите отыскать ошибку пожалуйста. program...

Ошибка деления на ноль
{$M $1000,0,0} Program Scan_code; Uses Crt,Dos; Var OldKey: Procedure;...

ошибка деление на ноль
пишет деление на ноль program l4; const eps=0.1 var ...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
1
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru