Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-06 - 43,50402/model5.pgm
There are 2 other files named model5.pgm in the archive. Click here to see a list.
      SUBROUTINE MODEL5 (NPROB, B, FY, NOB, NC, X, NVARX, NOBMAX, NCMAX,08M50010
     1 KTOU)                                                            08M10011
C          OXIDES MEASURED IN MOL PERCENTS.    D' = D1 TO D5            08M50020
C --- 5 OXIDES( - MIN VALUES) VS D AND A.  B, G AND H ARE COMMON.       08M50030
C                                                                       08M50040
      IMPLICIT REAL*8(A-H,O-Z)                                          08M50050
      DIMENSION  B(NCMAX), FY(NOBMAX), X(NVARX,NOBMAX)                  08M50060
      RB  = 1.0/B(1)                                                    08M50070
      DO 10  J=1,NOB                                                    08M50080
      X1 = (X(2,J) - 20.9) / 7.43                                       08M50090
      X2 = (X(3,J) - 2.03) / 7.43                                       08M50100
      X3 = (X(4,J) - 0.44) / 7.43                                       08M50110
      X4 = (X(5,J) - 66.1) / 7.43                                       08M50120
      X5 = (X(6,J) - 3.10) / 7.43                                       08M50130
      X12 = X1*X2                                                       08M50140
      X13 = X1*X3                                                       08M50150
      X14 = X1*X4                                                       08M50160
      X24 = X2*X4                                                       08M50170
      X23 = X2*X3                                                       08M50180
      X34 = X3*X4                                                       08M50190
      A =  B( 4)*X1  - B( 5)*X2  + B( 6)*X3  + B( 7)*X4  + B( 8)*X5     08M50200
     1   + B( 9)*X12 - B(10)*X13 - B(11)*X14 - B(12)*X23 + B(13)*X24    08M50210
     2   - B(14)*X34                                                    08M50220
      D =  B(15)*X1  + B(16)*X2  + B(17)*X3  + B(18)*X4  + B(19)*X5     08M50230
      C =  B( 2)*D - B( 3)                                              08M50240
      F = (A/C)**RB                                                     08M50250
      FY(J) = D - C + ( A / ( (X(1,J) + F )**B(1) )    )                08M50260
   10 CONTINUE                                                          08M50270
      RETURN                                                            08M50280
      END                                                               08M50290