Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Turbo Pascal / Stud_1_1 / Metlab / LABOR3 / LABOR3A

.PAS
Скачиваний:
92
Добавлен:
03.03.2016
Размер:
2.82 Кб
Скачать
Program Labor3a;
Uses Crt;
Const Nmax = 100;
eps = 0.001;
Type Ar = array[1..Nmax] of real;
Var
Sy,Sz : shortint; { §­ ЄЁ н«Ґ¬Ґ­в®ў y[i] Ё z[i] }
i, { Ї а ¬Ґва жЁЄ«  }
n, { Є®«ЁзҐбвў® н«Ґ¬Ґ­в®ў ¬ ббЁў  }
m, { Є®«-ў® б®ўЇ ¤ ойЁе Ї® §­ Єг н«Ґ¬Ґ­в®ў }
k : byte; { ўбЇ®¬®Ј вҐ«м­ п ЇҐаҐ¬Ґ­­ п }
Ly,Lz : real; { ¤«Ё­л ўҐЄв®а®ў }
ch : char; { бЁ¬ў®« ­ ¦ в®© Є« ўЁиЁ }
X,Y,Z : Ar; { ¬ ббЁўл }
F : text; { Ёб室­л© д ©« }
Begin

{ ‚ў®¤ Ёб室­ле ¤ ­­ле }
ClrScr;
Assign(F,'X.dat');
Reset(F);
n:=0;
While not SeekEof(F) do
Begin
Inc(n);
Read(F,x[n]);
End;
Close(F);

{ ‚лў®¤ Ёб室­ле ¤ ­­ле ­  нЄа ­ }
Writeln('Њ ббЁў X n=',n);
k:=0;
For i:=1 to n do
Begin
k:=k+1;
If k<5 then
Write(x[i]:8:2,' ')
Else
Begin
k:=0;
Writeln(x[i]:8:2);
End
End;
If k>0 then
Writeln;

{ ‚лзЁб«Ґ­ЁҐ н«Ґ¬Ґ­в®ў ¬ ббЁў®ў Y Ё Z }
For i:=1 to n do
Begin
if x[i]<i then
y[i]:=2*ln(sqr(x[i]+1)+(sqr(sin(sqr(x[i])))+exp(ln(abs(x[i])/3))));
If (y[i]-1)>eps then
z[i]:=((x[i]+2)*x[i]+1)*x[i]-1
Else
y[i]:=(x[i]*sqrt(sqr(sin(x[i]))+sqr(cos(sqr(x[i]))+1)/(sqrt(2*x[i]+1)*(ln(x[i]+1))+1)));
End;

{ ‚лў®¤ ¬ ббЁў  Y ­  нЄа ­ }
Writeln('Њ ббЁў Y');
k:=0;
For i:=1 to n do
Begin
k:=k+1;
If k<5 then
Write(y[i]:8:2,' ')
Else
Begin
k:=0;
Writeln(y[i]:8:2);
End
End;
If k>0 then
Writeln;

{ ‚лў®¤ ¬ ббЁў  Z ­  нЄа ­ }
Writeln('Њ ббЁў Z');
k:=0;
For i:=1 to n do
Begin
k:=k+1;
If k<5 then
Write(z[i]:8:2,' ')
Else
Begin
k:=0;
Writeln(z[i]:8:2);
End
End;
If k>0 then
Writeln;

{ ЋЇаҐ¤Ґ«Ґ­ЁҐ ¤«Ё­ ўҐЄв®а®ў Y Ё Z }
m:=0;
For i:=1 to n do
S:=(sqrt(sqr(x[c]-x[i])+sqr(y[c]-y[i])));
S:=sqrt(S);
Writeln('S= ',S:8:2);
For i:=1 to n do
if (S>r) and (S<rb);
m:=m+1;
S:=m/n*100
Writeln('S= ',S:8:2);

{ ЋЇаҐ¤Ґ«Ґ­ЁҐ Є®«ЁзҐбвў  б®ўЇ ¤ ойЁе Ї® §­ Єг н«Ґ¬Ґ­в®ў Y Ё Z }
m:=0;
For i:=1 to n do
Begin
If abs(y[i])<=eps then
Sy:=0
Else
If y[i]>eps then
Sy:=1
Else
Sy:=-1;
If abs(z[i])<=eps then
Sz:=0
Else
If z[i]>eps then
Sz:=1
Else
Sz:=-1;
If Sy=Sz then
Inc(m);
End;
Writeln('Љ®«-ў® б®ўЇ ¤ ойЁе Ї® §­ Єг н«-в®ў m = ',m);

{ ЏаЁ®бв ­®ўЄ  а Ў®вл Їа®Ја ¬¬л }
Repeat
ch:=ReadKey;
Until ord(ch)=27;

End.

Соседние файлы в папке LABOR3