Вывод результата на экран или в другую программу

PROGRAM CHTR_TINH_CU;
{$N+,E+,F+}
    uses crt,dos,printer,vaosl1,sogiann1,TextCrt,
         thtotd1,inkqtk1;
    type
    mat1c=array[1..10] of real;
    mat1c1=array[1..20] of real;
    mat1c2=array[1..146] of real;
    mat1c3=array[1..10] of real;
    mat1c4=array[1..240] of real;
    mat1c5=array[1..180] of real;
    mat1c6=array[1..30] of real;
    mat1c7=array[1..300] of real;
    mat1c8=array[1..150] of real;
    mat2c1=array[1..10,1..7] of real;
    mat2c2=array[1..10,1..3] of real;
    mat2c3=array[1..4,1..2] of real;
    mat2c4=array[1..6,1..60] of real;
   { mat2c5=array[1..100,1..6] of real;}
    const
    ngc=10;mgc=7;           ngc1=10;mgc1=7;
    nqc=10;mqc=3;           nvlc=4;mvlc=2;
    nac=120;nbc=20;         nkc=58;nnc=150;
    var
    slc:mat1c;              sld,sld1:mat2c1;
    sltt:mat2c2;            slvl:mat2c3;
    s:mat1c2;               hsl:mat1c1;
    qs:mat1c3;              rs:mat1c5;
    a:mat1c4;               b:mat1c6;
    sr:mat2c4;              g,r:mat1c7;
    {cm:mat2c5;}              al,m,d1,d2:mat1c8;
    ni:longint;
    ng,ng1,nq,n:integer;
    hc,ha,dh,ksn,epp,hh,ghmm:real;
    mmax,mmin,w,rr,rt,f,sf:real;
    nit1,nit,hskt,i,j,k,iik,jjk,pts,bkij:integer;
    ch,mni:char;
    tenfail,fal:string;
    kmen:string[80];
   PROCEDURE GIAI_CU_NN;
   label
   30;
   var
   a1,a2,a3:real;
   rm,m2,m3:real;
   i,j,k,sx:integer;
   ijk:longint;
   hbeg,mbeg,sbeg,s100beg:word;
   hend,mend,send,s100end:word;
   gio,phut,giay,phgiay:real;
   gio1,phut1,giay1,phgiay1:integer;  
   PROCEDURE GIAI_CU_KHNN;
    label
    10,15;
    var
    i,j,k:integer;
    ka,att:real;
  PROCEDURE TAO_CAC_THSOKHNN;{Khong dung so nn}
    var
      i,j,k,jj,jk:integer;
      ch:char;
    BEGIN
      jk:=0;jj:=1;
      for i:=1 to 3*(ng+ng1) do
           begin 
           jk:=jk+1;if jk>3 then begin jj:=jj+8;jk:=1;end;
           A[jj+jk]:=S[2*i-1];
           end;
      for i:=1+3*(ng+ng1) to 3*(ng+ng1)+nq+4 do
          begin
            jk:=i-3*(ng+ng1);
            if i<=3*(ng+ng1)+nq then B[2*jk-1]:=S[2*i-1];B[2*jk]:=QS[jk];
            if i=3*(ng+ng1)+nq+4 then f:=pi*s[2*i-1]*s[2*i-1]/4;
            if i=3*(ng+ng1)+nq+3 then w:=s[2*i-1];
            if i=3*(ng+ng1)+nq+2 then rr:=s[2*i-1];
            if i=3*(ng+ng1)+nq+1 then rt:=s[2*i-1];
         end;
    END;{Tao xong cac thong so onn}
    {--------------------------------------------------------}
   BEGIN
    TAO_CAC_THSOKHNN;
    XAY_DUNG_SO_LIEU_TT(ng,ng1,hc,a,hsl);
    GIAI_CU3(a,b,g,r,al,m,sr,hskt,
            hc,ha,dh,ksn,hh,w,rr,rt,f,sf,ng,ng1,nq,n);
    if hskt=1 then
    begin
      writeln;write('ѓ«гЎЁ** Ї®Ја㦥**®Ј® иЇг*в* *Ґ ¤®бв®в®з**');
      ch:=readkey;halt;
    end;
    ghmm:=sr[2,56];
    i:=1;
    Repeat
    sr[i,51]:=0.25*sr[i,56];
    sr[i,51]:=int(sr[i,51]);j:=1;
    ka:=0.08*sr[i,56];
    if abs(ka)>=1 then att:=round(ka) else att:=ka;
    k:=1;
    repeat
    if abs(ka)>=1 then if abs(att)<10 then
    begin sr[i,52]:=att*j;goto 15;end else
    begin att:=round(att/10);j:=j*10;end;
    if abs(ka)>=1 then if abs(att)>1 then goto 10 else
    begin sr[i,52]:=att*j;goto 15;end;
    if abs(ka)<1 then begin att:=att*10;j:=j*10;end;
    if abs(att)<1 then goto 10 else
    begin att:=round(att);sr[i,52]:=att/j;goto 15;end;
    10: k:=k+1;
    until k>3;
    15:i:=i+1;
    until i>6;
 {   for k:=1 to 6 do
    write(sr[k,51]:10,sr[k,52]:10);writeln;
    for k:=1 to 6 do
    write(sr[k,56]:10);writeln;
    ch:=readkey;halt;}
   END;
{-----------------------------------------------------------}
   BEGIN
   for i:=1 to 6 do
   for j:=1 to 60 do
     sr[i,j]:=0;
     bkij:=0;hskt:=0;
   GIAI_CU_KHNN;
   nit:=0;nit1:=0;
   gettime(hbeg,mbeg,sbeg,s100beg);
   putwindow(2,2,79,24,white,blue,doubleborder,2);
   gotoxy(15,1);write('CHUONG TRINH TINH CU THEO PHUONG PHAP MONTE-CARLO');
   gotoxy(17,21);write('CHUONG TRINH LAM VIEC CUNG VOI FILE SL1.DAT');
   putwindow(15,5,65,20,lightcyan,lightred,doubleborder,3);cursoroff;
 For ijk:=1 to ni do  {bat dau vong thu}
   Begin
   textcolor(11+128);
   gotoxy(15,2);write('CHuong trinh dang chay');
   textcolor(11);
   gettime(hend,mend,send,s100end);
   gio:=hend-hbeg;gio1:=trunc(gio);
   phut:=mend-mbeg;phut1:=trunc(phut);
   giay:=send;giay1:=trunc(giay);
   phgiay:=s100end;phgiay1:=trunc(phgiay);
   gotoxy(18,4);
   write(gio1:2,':',phut1:2,':',giay1:2,'.',phgiay1:2);
   gotoxy(17,5);writeln('So lan thu   :',ijk:8);
   TAO_CAC_THSONN1(ng,ng1,nq,s,a,rs,b,qs,w,rr,rt,f);
   XAY_DUNG_SO_LIEU_TT(ng,ng1,hc,a,hsl);
 
  GIAI_CU3(a,b,g,r,al,m,sr,hskt,
           hc,ha,dh,ksn,hh,w,rr,rt,f,sf,ng,ng1,nq,n);
   gotoxy(8,7);write('Ђ*ЄҐа.аҐ*ЄжЁп     :',sr[1,56]:20);
   gotoxy(8,8);write('Њ*Є.¬®¬Ґ*в        :',sr[2,56]/0.65:20);
   gotoxy(8,10);write('Њ®¬Ґ*в *Ґб.бЇ®б®Ў.:',sr[3,56]:20);
   gotoxy(8,11);write('Ђ*ЄҐа.бЁ«* *Ґб.бЇ.:',sr[4,56]:20);
   gotoxy(8,13);write('—Ёб«® б®ЎлвЁ© *Ґгбв®©зЁў®бвЁ =',nit:4);
   if hskt=1 then
   begin nit:=nit+1;hskt:=0;goto 30;end;
   for j:=1 to 4 do
   begin
   a1:=sr[j,56];a2:=a1;
   for k:=57 to 60 do
   begin sr[j,k]:=sr[j,k]+a2;a2:=a2*a1; end;
   end;
   for j:=1 to 4 do
   begin
   sr[j,53]:=(sr[j,56]-sr[j,51])/sr[j,52];
   sx:=trunc(sr[j,53]);
   if sr[j,53]<0 then sr[j,54]:=sr[j,54]+1;
   if sr[j,53]>50 then sr[j,55]:=sr[j,55]+1;
   if sx>=0 then if sx<=50 then sr[j,sx]:=sr[j,sx]+1;
   end;
  30:for j:=5 to 6 do
   begin
   a1:=sr[j,56];a2:=a1;
   for k:=57 to 60 do
   begin sr[j,k]:=sr[j,k]+a2;a2:=a2*a1; end;
   end;
   for j:=5 to 6 do
   begin
   sr[j,53]:=(sr[j,56]-sr[j,51])/sr[j,52];
   sx:=trunc(sr[j,53]);
   if sr[j,53]<0 then sr[j,54]:=sr[j,54]+1;
   if sr[j,53]>50 then sr[j,55]:=sr[j,55]+1;
   if sx>=0 then if sx<=50 then sr[j,sx]:=sr[j,sx]+1;
   end;
 end;{ket thuc vong thu ni}
   getwindow(3);getwindow(2);
   for i:=1 to 4 do
   begin
   sr[i,57]:=sr[i,57]/(ni-nit);
   rm:=sr[i,57]*sr[i,57];
   m2:=sr[i,58]/(ni-nit)-rm;m3:=sqrt(abs(m2));
   sr[i,60]:=(sr[i,60]/(ni-nit)-4*sr[i,59]*sr[i,57]/(ni-nit)+6*rm*sr[i,58]/(ni-nit)-3*rm*rm)/(m2*m2)-3;
   sr[i,59]:=(sr[i,59]/(ni-nit)-3*sr[i,57]*sr[i,58]/(ni-nit)+2*sr[i,57]*rm)/(m2*m3);
   sr[i,58]:=m3;
   end;
   for i:=5 to 6 do
   begin
   sr[i,57]:=sr[i,57]/(ni);
   rm:=sr[i,57]*sr[i,57];
   m2:=sr[i,58]/(ni)-rm;m3:=sqrt(abs(m2));
   sr[i,60]:=(sr[i,60]/(ni)-4*sr[i,59]*sr[i,57]/(ni)+6*rm*sr[i,58]/(ni)-3*rm*rm)/(m2*m2)-3;
   sr[i,59]:=(sr[i,59]/(ni)-3*sr[i,57]*sr[i,58]/(ni)+2*sr[i,57]*rm)/(m2*m3);
   sr[i,58]:=m3;
   end;
   for i:=0 to 3*(ng+ng1)+nq+4-1 do
   begin
   k:=1+2*i;rs[k]:=rs[k]/(ni);
   rs[k+1]:=sqrt(abs(rs[k+1]/(ni)-rs[k]*rs[k]));
   end;
  END;    {xong tinh cu nn}
  {----------------------------------------------------------}
  procedure cuasodautien;
  begin
  putwindow(1,10,80,16,white,blue,singleborder,1);
  gotoxy(5,2);write('            Џа®Ја*¬¬* а*бзҐв* Ў®«мўҐаЄ* Ї® ЊЋЌ’… ЉЂђ‹Ћ');
  gotoxy(5,3);write('Ќ®¬Ґа ў*аЁ**в*:');readln(fal);
  gotoxy(5,4);write('Љ®¬¬Ґ*в*аЁЁ ');readln(kmen);cursoroff;
  getwindow(1);
  end;
  procedure cuasomenieu;
  begin
  putwindow(10,5,70,20,white,blue,singleborder,1);
  gotoxy(2,2);write(' [F1]: ‚ў®¤ Ёб室*ле ¤***ле');
  gotoxy(2,3);write(' [F2]: ђ*бзҐв');
  gotoxy(2,4);write(' [F3]: ЏҐз*вм १г«мв*в®ў ');
  gotoxy(2,7);write(' [F9]: „агЈ®© ў*аЁ**в ');
  gotoxy(2,8);write(' [F10]:‚л室 ');cursoron;
  gotoxy(2,13);write(' ‚лЎЁа*©вҐ Ё **¦¬ЁвҐ ўлЎа***лҐ Є«*ўЁиЁ');cursoron;
  end;
 
BEGIN
    clrscr;textbackground(0);pts:=0;
    cuasodautien;cuasomenieu;tenfail:='kq'+fal+'.res';
 
 repeat
 mni:=readkey;if mni=#0 then
 Begin
 mni:=readkey;
 case mni of
 #59:begin
         vaotufai(ngc,mgc,ngc1,mgc1,nqc,mqc,nvlc,mvlc,'sl1.dat',
            slc,sld,sld1,sltt,slvl);
{        vao(ngc,mgc,ngc1,mgc1,nqc,mqc,nvlc,mvlc,
            slc,sld,sld1,sltt,slvl);}
         ni:=trunc(slc[1]);ng:=trunc(slc[2]);
         ng1:=trunc(slc[3]);nq:=trunc(slc[4]);
         sf:=slc[5];hc:=slc[6];ha:=slc[7];dh:=slc[8];
         ksn:=slc[9];hh:=slc[10];
         xemsua(mgc,mgc1,mqc,nvlc,mvlc,slc,sld,sld1,sltt,slvl,
                     ng,ng1,nq,ni,sf,hc,ha,dh,ksn,hh);
         catsl(ngc,mgc,ngc1,mgc1,nqc,mqc,nvlc,mvlc,'sl1.dat',
            slc,sld,sld1,sltt,slvl);
         LAP_CAC_BANG_SO_LIEU(mgc,mgc1,mqc,nvlc,mvlc,
                   slc,sld,sld1,sltt,slvl,hsl,s,qs,ng,ng1,nq);
         cuasomenieu;
     end;
  #60:begin pts:=1;
      for iik:=1 to 2*(3*(ng+ng1)+nq+4) do
      RS[iik]:=0;
      for iik:=1 to 6 do
      for jjk:=1 to 60 do
      sr[iik,jjk]:=0;
      n:=trunc(hh/dh)+ng+nq+ng1;
      putwindow(1,1,80,25,red,cyan,emptyborder,1);
      GIAI_CU_NN;
      LUU_KET_QUA_NN(ng,ng1,nq,nit,nit1,ni,slc,rs,sr,tenfail);
      getwindow(1);cuasomenieu;
      end;
 #61:begin
      DOC_KET_QUA_NN(ng,ng1,nq,ni,nit,nit1,slc,rs,sr,tenfail);
      inkqnn(ng,ng1,nq,nit,nit1,ni,slc,rs,sr,kmen);cuasomenieu;
     end;
 
 #67:begin cuasodautien;cuasomenieu;
           tenfail:='kq'+fal+'.res';
     end;
 end;    {ket thuc lenh case of}
 End;
 until mni=#68;
END.