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