{Описание констант} Const phi0=30; // Начальный угол наклона маятника w0=0; // Начальная угловая скорость маятника l=400; // Длина нити r=8; // Радиус груза light=2; // Радиус блика k=0.02; // Коэффициент сопротивления m=1; // Масса груза xp=320; // Координата x точки подвеса yp=10; // Координата y точки подвеса g=9.8; // Ускорение свободного падения tau=0.4; // Шаг по времени
{Описание переменных} Var x, y, xx, yy, w, eps, phi: Real;
{Процедура инициализации графического режима 640x480} Procedure OpenGraph; Var Driver, Mode, ErrorCode: Integer; Begin Driver:=Detect; InitGraph(Driver, Mode, 'c:\bp\bgi'); ErrorCode:=GraphResult; If ErrorCode <> grOK Then Halt(1); End;
{Процедура исключения мерцаний изображений} Procedure Glimmer; Begin Repeat Until Port[$3da] And 8 <> 0; End;
Const CSize1=7; CSize2=6; CSize3=5; type TMas=array[,] of integer; var A1:TMas; A2:TMas; A3:TMas; x, y:byte;
procedure zapoln(var T:TMas;Size:integer); begin var i,j:integer; setlength(T,Size+1,Size+1); for i:=1 to Size do begin for j:=1 to Size do begin T[i,j]:=random(50); write(T[i,j]:3); end; writeln; end; writeln; end;
function zero(T:TMas):integer; var i,j,n:integer; begin n:=length(T,0)-1; for i:=1 to n do for j:=1 to n do if T[i,j]=0 then result:=result+1; end;
function interval(T:TMas):integer; var i,j,n:integer; begin n:=length(T,0)-1; for i:=1 to n do for j:=1 to n do if (1<=T[i,j]) and (T[i,j]<=12) then result:=result+1; end;
begin zapoln(A1,CSize1); zapoln(A2,CSize2); zapoln(A3,CSize3); x:=zero(A1)+zero(A2)+zero(A3); y:=interval(A1)+interval(A2)+interval(A3); writeln('Kol-vo 0 =', x); writeln('Kol-vo v intervale=', y); end.
{Программа
СВОБОДНЫЕ КОЛЕБАНИЯ КРУГОВОГО МАТЕМАТИЧЕСКОГО МАЯТНИКА}
Program Math_01;
{Подключение модулей подпрограмм}
Uses Crt, Graph;
{Описание констант}
Const phi0=30; // Начальный угол наклона маятника
w0=0; // Начальная угловая скорость маятника
l=400; // Длина нити
r=8; // Радиус груза
light=2; // Радиус блика
k=0.02; // Коэффициент сопротивления
m=1; // Масса груза
xp=320; // Координата x точки подвеса
yp=10; // Координата y точки подвеса
g=9.8; // Ускорение свободного падения
tau=0.4; // Шаг по времени
{Описание переменных}
Var x, y, xx, yy, w, eps, phi: Real;
{Процедура инициализации графического режима 640x480}
Procedure OpenGraph;
Var Driver, Mode, ErrorCode: Integer;
Begin
Driver:=Detect;
InitGraph(Driver, Mode, 'c:\bp\bgi');
ErrorCode:=GraphResult;
If ErrorCode <> grOK Then Halt(1);
End;
{Процедура исключения мерцаний изображений}
Procedure Glimmer;
Begin
Repeat
Until Port[$3da] And 8 <> 0;
End;
{Процедура вывода изображения маятника на экран}
Procedure Display;
Begin
Glimmer;
SetColor(Black);
SetFillStyle(SolidFill, Black);
Line(xp, yp, Round(xx), Round(yy));
FillEllipse(Round(xx), Round(yy), r, r);
SetLineStyle(SolidLn, 0, NormWidth);
SetColor(DarkGray);
Line(xp, yp, Round(x), Round(y));
SetLineStyle(SolidLn, 0, NormWidth);
SetColor(DarkGray);
SetFillStyle(SolidFill, DarkGray);
FillEllipse(Round(x), Round(y), r, r);
SetColor(LightGray);
SetFillStyle(SolidFill, White);
FillEllipse(Round(x)-3, Round(y)-3, light, light);
End;
{Основная часть программы}
Begin
OpenGraph;
SetFillStyle(SolidFill, DarkGray);
Bar(xp-20, yp-5, xp+20, yp-1);
phi:=phi0*Pi/180;
w:=w0;
x:=xp+l*sin(phi);
y:=yp+l*cos(phi);
Repeat
Begin
xx:=x;
yy:=y;
eps:=-k*w/m-g*sin(phi)/l;
w:=w+eps*tau;
phi:=phi+w*tau;
x:=xp+l*sin(phi);
y:=yp+l*cos(phi);
Display;
End;
Until KeyPressed;
CloseGraph;
End.