Trailing-Edge
-
PDP-10 Archives
-
decuslib10-06
-
43,50402/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