SUBROUTINE MODEL4 (NPROB, B, FY, NOB, NC, X, NVARX, NOBMAX, NCMAX,07M40010
1 KTOU) 07M10011
C 07M40020
C --- 5 OXIDES( - MIN VALUES) VS D AND A. B, G AND H ARE COMMON. 07M40030
C 07M40040
IMPLICIT REAL*8(A-H,O-Z) 07M40050
DIMENSION B(NCMAX), FY(NOBMAX), X(NVARX,NOBMAX) 07M40060
C 07M40070
RB = 1.0/B(1) 07M40080
DO 10 J=1,NOB 07M40090
X1 = (X(2,J) - 21.0) / 10.66 07M40100
X2 = (X(3,J) - 3.48) / 10.66 07M40110
X3 = (X(4,J) - 1.18) / 10.66 07M40120
X4 = (X(5,J) - 61.6) / 10.66 07M40130
X5 = ( 97.92 - (X(2,J) + X(3,J) + X(4,J) + X(5,J))) / 10.66 07M40140
X6 = X1*X3 07M40150
X7 = X1*X4 07M40160
X8 = X2*X4 07M40170
X9 = X3*X4 07M40180
X10 = X1*X2 07M40190
D = B( 2)*X1 + B( 3)*X2 + B( 4)*X3 + B( 5)*X4 + B( 6)*X5 07M40200
1 - B( 7)*X6 + B( 8)*X8 07M40210
A = B( 9)*X1 + B(10)*X2 + B(11)*X3 + B(12)*X4 + B(13)*X5 07M40220
1 + B(14)*X10- B(15)*X7 07M40230
C = B(16)*D - B(17) 07M40240
F = (A/C)**RB 07M40250
FY(J) = D - C + ( A / ( (X(1,J) + F )**B(1) ) ) 07M40260
10 CONTINUE 07M40270
RETURN 07M40280
END 07M40290