Специалист
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 27,807
Сказал(а) спасибо: 340
Поблагодарили 583 раз(а) в 208 сообщениях
Репутация: 113184
|
Цитата:
Сообщение от brokilone
Нужно написать две слишком сложные для меня программы в Pascal.
|
Ну вот, первая (черт, весь день с этой ерундой провозился):
Код:
Uses CRT;
Const
Xc=32.0;
Yc=-40.0;
Xe=-21.0;
Ye=40.0;
R=33.0;
Ae=16.0;
Be=38.0;
Var
X0c,Y0c,X0e,Y0e,X1,Y1,X2,Y2,X3,Y3,Step:Real;
i,Nt,k:Byte;
N:Word;
Rect,Trap,Simp:Array[1..100] of Real;
Np:Array[1..100] of Word;
Function Ysl(Xsl:Real):Real;
begin
Ysl:=Y1+(Xsl-X1)/(X2-X1)*(Y2-Y1);
end;
Function Ycr1(Xcr:Real):Real;
begin
Ycr1:=Y0c+Sqrt(Sqr(R)-Sqr(Xcr-X0c));
end;
Function Ycr2(Xcr:Real):Real;
begin
Ycr2:=Y0c-Sqrt(Sqr(R)-Sqr(Xcr-X0c));
end;
Function Yel(Xel:Real):Real;
begin
Yel:=Y0e-Be*Sqrt(1-Sqr((Xel-X0e)/Ae));
end;
Function F1(Xf:Real):Real;
begin
F1:=Ycr1(Xf)-Ysl(Xf);
end;
Function F2(Xf:Real):Real;
begin
F2:=Yel(Xf)-Ycr2(Xf);
end;
Function F3(Xf:Real):Real;
begin
F3:=Ysl(Xf)-Ycr2(Xf);
end;
Function Rectangle1:Real;
var
i:Integer;
Sum,h,a,b:Real;
begin
a:=X1;
b:=0;
h:=(b-a)/N;
Sum:=0;
for i:=0 to N-1 do
Sum:=Sum+h*F1(a+h*i);
Rectangle1:=Sum;
end;
Function Rectangle2:Real;
var
i:Integer;
Sum,h,a,b:Real;
begin
a:=0;
b:=X2;
h:=(b-a)/N;
Sum:=0;
for i:=0 to N-1 do
Sum:=Sum+h*F2(a+h*i);
Rectangle2:=Sum;
end;
Function Rectangle3:Real;
var
i:Integer;
Sum,h,a,b:Real;
begin
a:=X2;
b:=X3;
h:=(b-a)/N;
Sum:=0;
for i:=0 to N-1 do
Sum:=Sum+h*F3(a+h*i);
Rectangle3:=Sum;
end;
Function Trapezium1:Real;
var
i:Integer;
Sum,h,a,b:Real;
begin
a:=X1;
b:=0;
h:=(b-a)/N;
Sum:=(F1(a)+F1(b))/2*h;
for i:=1 to N-1 do
Sum:=Sum+h*F1(a+h*i);
Trapezium1:=Sum;
end;
Function Trapezium2:Real;
var
i:Integer;
Sum,h,a,b:Real;
begin
a:=0;
b:=X2;
h:=(b-a)/N;
Sum:=(F2(a)+F2(b))/2*h;
for i:=1 to N-1 do
Sum:=Sum+h*F2(a+h*i);
Trapezium2:=Sum;
end;
Function Trapezium3:Real;
var
i:Integer;
Sum,h,a,b:Real;
begin
a:=X2;
b:=X3;
h:=(b-a)/N;
Sum:=(F3(a)+F3(b))/2*h;
for i:=1 to N-1 do
Sum:=Sum+h*F3(a+h*i);
Trapezium3:=Sum;
end;
Function Simpson1:Real;
var
i,Ns:Integer;
Sum,h,a,b:Real;
begin
a:=X1;
b:=0;
Ns:=N div 2;
h:=(b-a)/(2*Ns);
Sum:=(F1(a)+F1(b))/3*h;
for i:=1 to Ns do
Sum:=Sum+4.0*h/3.0*F1(a+h*(2*i-1));
for i:=2 to Ns do
Sum:=Sum+2.0*h/3.0*F1(a+h*(2*i-2));
Simpson1:=Sum;
end;
Function Simpson2:Real;
var
i,Ns:Integer;
Sum,h,a,b:Real;
begin
a:=0;
b:=X2;
Ns:=N div 2;
h:=(b-a)/(2*Ns);
Sum:=(F2(a)+F2(b))/3*h;
for i:=1 to Ns do
Sum:=Sum+4.0*h/3.0*F2(a+h*(2*i-1));
for i:=2 to Ns do
Sum:=Sum+2.0*h/3.0*F2(a+h*(2*i-2));
Simpson2:=Sum;
end;
Function Simpson3:Real;
var
i,Ns:Integer;
Sum,h,a,b:Real;
begin
a:=X2;
b:=X3;
Ns:=N div 2;
h:=(b-a)/(2*Ns);
Sum:=(F3(a)+F3(b))/3*h;
for i:=1 to Ns do
Sum:=Sum+4.0*h/3.0*F3(a+h*(2*i-1));
for i:=2 to Ns do
Sum:=Sum+2.0*h/3.0*F3(a+h*(2*i-2));
Simpson3:=Sum;
end;
BEGIN
ClrScr;
X0c:=-1.0;
Y0c:=Yc+Sqrt(sqr(R)-sqr(X0c));
Y0e:=2.0;
X0e:=Xe+Ae*Sqrt(1.0-sqr(Y0e/Be));
Writeln('X0c=',X0c:9:5,' - X координата центра окружности');
Writeln('Y0c=',Y0c:9:5,' - Y координата центра окружности');
Writeln('X0e=',X0e:9:5,' - X координата центра эллипса');
Writeln('Y0e=',Y0e:9:5,' - Y координата центра эллипса');
Writeln;
Step:=1;
X1:=X0e-Ae;
For i:=1 to 5 do
begin
Step:=Step/10;
Repeat
X1:=X1+Step;
Y1:=Y0e+Be*Sqrt(1-sqr((X1-X0e)/Ae));
Until Sqr(X1-X0c)+Sqr(Y1-Y0c)>Sqr(R);
X1:=X1-Step;
end;
Y1:=Y0e+Be*Sqrt(1-sqr((X1-X0e)/Ae));
Writeln('X1= ',X1:9:5,' - X координата первой точки пересечения');
Writeln('Y1= ',Y1:9:5,' - Y координата первой точки пересечения');
Y2:=0;
X2:=X0e+Ae*Sqrt(1-sqr(Y0e/Be));
Writeln('X2= ',X2:9:5,' - X координата второй точки пересечения');
Writeln('Y2= ',Y2:9:5,' - Y координата второй точки пересечения');
Step:=1;
X3:=X2;
For i:=1 to 5 do
begin
Step:=Step/10;
Repeat
X3:=X3+Step;
Y3:=Ysl(X3);
Until Sqr(X3-X0c)+Sqr(Y3-Y0c)>Sqr(R);
X3:=X3-Step;
end;
Y3:=Ysl(X3);
Writeln('X3= ',X3:9:5,' - X координата третьей точки пересечения');
Writeln('Y3= ',Y3:9:5,' - Y координата третьей точки пересечения');
Writeln;
Writeln('Нажмите "Enter" для продолжения...');
Readln;
ClrScr;
k:=0;
Repeat
Write('Введите количество точек (0 для выхода из программы) ');
Readln(N);
If N>0 then
Begin
ClrScr;
Inc(k);
Np[k]:=N;
For i:=1 to 58 do Write('-');
Writeln;
Writeln('| Число разбиений | Результат |');
Writeln('| |--------------------------------------|');
Writeln('| n | A | B | C |');
For i:=1 to 58 do Write('-');
Writeln;
Rect[k]:=Rectangle1+Rectangle2+Rectangle3;
Trap[k]:=Trapezium1+Trapezium2+Trapezium3;
Simp[k]:=Simpson1+Simpson2+Simpson3;
For i:=1 to k do
Writeln('| ',Np[i]:4,' |',Rect[i]:12:6,'|',Trap[i]:12:6,'|',Simp[i]:12:6,'|');
For i:=1 to 58 do Write('-');
Writeln;
End;
Until N=0;
END.
|