Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-06 - 43,50402/model3.pgm
There are 2 other files named model3.pgm in the archive. Click here to see a list.
      SUBROUTINE MODEL3 (NPROB, B, FY, NOB, NC, X, NVARX, NOBMAX, NCMAX,06M30010
     1 KTOU)                                                            06M10011
C                                                                       06M30020
C --- 31 COEFFICIENTS, 14 INDICATOR VARIABLES FOR D AND A               06M30030
C --- B, G AND H ARE COMMON.                                            06M30040
C                                                                       06M30050
      IMPLICIT REAL*8(A-H,O-Z)                                          06M30060
      DIMENSION  B(NCMAX), FY(NOBMAX), X(NVARX,NOBMAX)                  06M30070
C                                                                       06M30080
      RB  = 1.0/B(1)                                                    06M30090
      DO 10  J=1,NOB                                                    06M30100
      D = B( 2)*X( 2,J) + B( 3)*X( 3,J) + B( 4)*X( 4,J) + B( 5)*X( 5,J) 06M30110
     1  + B( 6)*X( 6,J) + B( 7)*X( 7,J) + B( 8)*X( 8,J) + B( 9)*X( 9,J) 06M30120
     2  + B(10)*X(10,J) + B(11)*X(11,J) + B(12)*X(12,J) + B(13)*X(13,J) 06M30130
     3  + B(14)*X(14,J) + B(15)*X(15,J)                                 06M30140
      A = B(16)*X( 2,J) + B(17)*X( 3,J) + B(18)*X( 4,J) + B(19)*X( 5,J) 06M30150
     1  + B(20)*X( 6,J) + B(21)*X( 7,J) + B(22)*X( 8,J) + B(23)*X( 9,J) 06M30160
     2  + B(24)*X(10,J) + B(25)*X(11,J) + B(26)*X(12,J) + B(27)*X(13,J) 06M30170
     3  + B(28)*X(14,J) + B(29)*X(15,J)                                 06M30180
      C = B(30)*D - B(31)                                               06M30190
      F = (A/C)**RB                                                     06M30200
      FY(J) = D - C + ( A / ( (X(1,J) + F )**B(1) )    )                06M30210
   10 CONTINUE                                                          06M30220
      RETURN                                                            06M30230
      END                                                               06M30240