46740 (607927), страница 3
Текст из файла (страница 3)
Додаток А
Текст програми:
Program MGVS;
{$N+}
Uses crt;
Const size = 60;
type
TMatrix=array[1..size, 1.. size] of double;
VAR lm,Tsr,d,Vs,Vn:double;
l,n,m,i,j,MaxN,MaxM: integer;
A:TMatrix;
FRez:text;
fname:string;
function deg(x:real;y:integer):double;
begin
deg:=Exp(Ln(x)*y);
end;
function fac(x:byte):real;
var i:integer;
rez:real;
begin
rez:=1;
for i:=1 to x do rez:= rez * i;
fac:=rez;
end;
function Pr(p_lm,p_Tsr,p_d,p_Vs,p_Vn:real; p_l,p_n,p_m:byte):double;
var P1,P2,p3,p4,p5:double;
i:Integer;
Begin
P1:=deg(p_Tsr*p_lm, p_n+p_l*p_m) / ( deg(p_n,p_l*p_m) * fac(p_n) );
p3:=0;
for i:=1 to n do p3:=p3+( deg(p_Tsr*p_lm,i)/(fac(i) ) );
p4:=deg(p_Tsr*p_lm,p_n+1)/( p_n*fac(p_n) ) ;
p5:= ( 1 - deg((p_Tsr*p_lm)/n ,p_l*p_m) ) / (1- ((p_Tsr*p_lm)/n) );
P2:=1/(1+p3+p4*p5);
Pr:=(p_d*p_lm)*(1-P1*P2)-(p_Vs*p_n)-(p_Vn*p_m);
end;
BEGIN
ClrScr;
Write('Vvedite fail rezultatov: '); Readln(Fname);
assign(frez,Fname);
Rewrite(frez);
Write('Enter lm = ');readln(lm);
Write('Enter Tsr = ');readln(Tsr);
Write('Enter d = ');readln(d);
Write('Enter Vs = ');readln(Vs);
Write('Enter Vn = ');readln(Vn);
Write('Enter l = ');readln(l);
Write('Enter max n = ');readln(n);
Write('Enter max m = ');readln(m);
for i:=1 to n do
for j:=1 to m do A[i,j]:=Pr(lm,Tsr,d,Vs,Vn,l,i,j);
for i:=1 to n do















