Google
 

Trailing-Edge - PDP-10 Archives - decus_20tap4_198111 - decus/20-0101/model2.pgm
There are 2 other files named model2.pgm in the archive. Click here to see a list.
      SUBROUTINE MODEL2 (NPROB, B, FY, NOB, NC, X, NVARX, NOBMAX, NCMAX,05M20010
     1 KTOU)                                                            05M10011
C                                                                       05M20020
C --- 43 COEFFICIENTS, 14 INDICATOR VARIABLES FOR D, A AND C            05M20030
C --- B IS COMMON.                                                      05M20040
C                                                                       05M20050
      IMPLICIT REAL*8(A-H,O-Z)                                          05M20060
      DIMENSION  B(NCMAX), FY(NOBMAX), X(NVARX,NOBMAX)                  05M20070
C                                                                       05M20080
      RB  = 1.0/B(1)                                                    05M20090
      DO 10  J=1,NOB                                                    05M20100
      D = B( 2)*X( 2,J) + B( 3)*X( 3,J) + B( 4)*X( 4,J) + B( 5)*X( 5,J) 05M20110
     1  + B( 6)*X( 6,J) + B( 7)*X( 7,J) + B( 8)*X( 8,J) + B( 9)*X( 9,J) 05M20120
     2  + B(10)*X(10,J) + B(11)*X(11,J) + B(12)*X(12,J) + B(13)*X(13,J) 05M20130
     3  + B(14)*X(14,J) + B(15)*X(15,J)                                 05M20140
      A = B(16)*X( 2,J) + B(17)*X( 3,J) + B(18)*X( 4,J) + B(19)*X( 5,J) 05M20150
     1  + B(20)*X( 6,J) + B(21)*X( 7,J) + B(22)*X( 8,J) + B(23)*X( 9,J) 05M20160
     2  + B(24)*X(10,J) + B(25)*X(11,J) + B(26)*X(12,J) + B(27)*X(13,J) 05M20170
     3  + B(28)*X(14,J) + B(29)*X(15,J)                                 05M20180
      C = B(30)*X( 2,J) + B(31)*X( 3,J) + B(32)*X( 4,J) + B(33)*X( 5,J) 05M20190
     1  + B(34)*X( 6,J) + B(35)*X( 7,J) + B(36)*X( 8,J) + B(37)*X( 9,J) 05M20200
     2  + B(38)*X(10,J) + B(39)*X(11,J) + B(40)*X(12,J) + B(41)*X(13,J) 05M20210
     3  + B(42)*X(14,J) + B(43)*X(15,J)                                 05M20220
      F = (A/C)**RB                                                     05M20230
      FY(J) = D - C + ( A / ( (X(1,J) + F )**B(1) )    )                05M20240
   10 CONTINUE                                                          05M20250
      RETURN                                                            05M20260
      END                                                               05M20270