home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / sigm / vol266 / otput3.for < prev    next >
Encoding:
Text File  |  1986-05-19  |  2.4 KB  |  56 lines
  1.         SUBROUTINE      OTPUT(NSET)
  2. C////////////////////////////////////////////////////////////////
  3. C/                                                              /
  4. C/      Program-id.     OTPUT                                   /
  5. C/      Date-written.   Jan. 23rd 1984                          /
  6. C/      File-name.      OTPUT3.FOR                              /
  7. C/      Remarks.        Subroutine OTPUT.FOR page 130           /
  8. C/                      Written by the programmer to perform    /
  9. C/                      calculations and provide additional     /
  10. C/                      output at the end of a simulation run.  /
  11. C/            This is the version for Example-3.    /
  12. C/                                                              /
  13. C////////////////////////////////////////////////////////////////
  14. C
  15.         INTEGER*4       NSET(6,1)
  16. C
  17.       COMMON /C1/ID,IM,INIT,JEVNT,JMNIT,MFA,MSTOP,MX,MXC,NCLCT,NHIST,
  18.      1           NOQ,NORPT,NOT,NPRMS,NRUN,NRUNS,NSTAT,OUT,SCALE,ISEED,
  19.      2           TNOW,TBEG,TFIN,MXX,NPRNT,NCRDR,NEP,VNQ(4)
  20. C
  21.       COMMON /C2/ATRIB(4),ENQ(4),INN(4),JCELS(5,22),KRANK(4),MAXNQ(4),
  22.      1           MFE(4),MLC(4),MLE(4),NCELS(5),NQ(4),PARAM(20,4),
  23.      2           QTIME(4),SSUMA(10,5),SUMA(10,5),NAME(6),NPROJ,MON,
  24.      3         NDAY,NYR,JCLR
  25. C
  26.       COMMON /C3/ XISYS,BUS,XL,XMU
  27. C
  28.     COMMON /C3/ XISYS,BUS
  29. C
  30. C       --- Compute theoretical and simulation values of performance
  31. C           measures for the queuing system.
  32. C
  33.         ETISS = SUMA(1,1) / SUMA(1,3)
  34.         EIDTS = (SSUMA(2,1) - SSUMA(2,2)) / (SUMA(1,3) - 1.0)
  35.         EWTS = SUMA(2,1) / SUMA(2,3)
  36.         EIDTC = PARAM(1,1) - PARAM(2,1)
  37.     EWTC = (1.0/PARAM(1,1)) / ((1.0/PARAM(2,1)) * (1.0/
  38.      1  PARAM(2,1) - 1.0/PARAM(1,1)))
  39.         ETISC = 1.0/(1.0/PARAM(2,1) - 1.0/PARAM(1,1))
  40.         YA = ETISS / (SSUMA(1,2) / SSUMA(1,1))
  41.         YS = ETISS - EWTS
  42.         WRITE(NPRNT,85)
  43.    85     FORMAT(/36X,'Simulated Value',4X,'Theoretical Value'/)
  44.         WRITE(NPRNT,90) EIDTS,EIDTC
  45.    90     FORMAT(10X,'Expected idle time',11X,F8.3,12X,F8.3)
  46.         WRITE(NPRNT,95) EWTS,EWTC
  47.    95     FORMAT(10X,'Expected waiting time',8X,F8.3,12X,F8.3)
  48.         WRITE(NPRNT,96) ETISS,ETISC
  49.    96     FORMAT(10X,'Expected time in system',6X,F8.3,12X,F8.3)
  50.         WRITE(NPRNT,97) YA,PARAM(1,1)
  51.    97     FORMAT(10X,'Expected arrival time',8X,F8.3,12X,F8.3)
  52.         WRITE(NPRNT,98) YS,PARAM(2,1)
  53.    98     FORMAT(10X,'Expected service time',8X,F8.3,12X,F8.3)
  54.         RETURN
  55.         END
  56.