“…0OO74348D0/ , C6/ . 00079824D0/ DATA C7/-.00029166D0/ IF (X.GT.3.0D0) GO TO 10 X2=X*X/9.D0 DBJ1=X* ( ( ( ( ( (A7*X2+A6) *X2+A5) *X2+A4) X2+A3) *X2+A2) *X2+A1) RETURN 10 X1=3.0D0/X Fl=:( ( ( ( (B7*X1+B6) *X1+B5) *X1+B4) X1+B3) *X1+B2) *X1+B1 THETA=X+ ( ( ( ( (C7*X1+C6) *X1+C5) *X1+C4) X1+C3) *X1+C2) Xl+Cl FORMAT (2X, 'CALCULATION OF THE SPREADING RESISTANCE FROM THE 1 RESISTIVITY ',/,2X, 'USING THE 33-POINT GAUSS-LAGUERRE SCHEME') GO TO (15,20,25), IBBC WRITE (6,17) SUBTH+N+DELTH , N*DELTH . C Qi^i^l^ili^l7^l]^i^l^l^l^t^l^l^liii^l^^^ilii**^l*************************** ************************** DO 100 1=1,33 90 W (I) =W (I) * (DB Jl (KA (I) THE 500 LOOP IS THE MAJOR LOOP WHICH CONTAINS THE CALCULATION OF THE C CORRECTION FACTOR INTEGRAL FOR THE N-POINT LAYERED STRUCTURE STARTING C FROM THE DEEPEST POINT TO THE TOP SURFACE.…”