1 : C FBG2 SOURCE CHAT 05/01/12 23:57:00 5004 2 : 3 : 4 : FUNCTION FBG2(T) 5 : IMPLICIT INTEGER(I-N) 6 : IMPLICIT REAL*8 (A-H,O-Z) 7 : 8 : t1 = t**2 9 : t2 = t1**2 10 : t4 = t2**2 11 : t5 = t4**2 12 : t6 = t5**2 13 : t7 = t6**2 14 : t9 = 0.2126494641379152D-89*t7*t2*t1+1.D0 15 : t10 = 1.D0/t9 16 : t11 = t2*t 17 : t12 = 0.1316872427983539D0*t11+1.D0 18 : t13 = 1.D0/t12 19 : t14 = -0.156710172162D4-3020.D0*t10-0.595089827838D4*t13-0.1146775 20 : #430405D2*t1+0.458710172162D3*t 21 : t16 = t9**2 22 : t19 = 1.D0/t16*t7*t11 23 : t20 = t12**2 24 : t22 = 1.D0/t20*t2 25 : t23 = 0.4495409671875527D-84*t19+0.3918286932266667D4*t22-0.229355 26 : #086081D2*t+0.458710172162D3 27 : t25 = 1.D0/t14 28 : t26 = sign(1.D0,t14)*t23*t25 29 : t27 = abs(t14) 30 : t28 = t14**2 31 : t31 = t27/t28*t23 32 : t33 = 0.458710172162D4+3020.D0*t10+0.595089827838D4*t13+0.11467754 33 : #30405D2*t1-0.458710172162D3*t 34 : t35 = t27*t25 35 : t36 = 1.D0/2.D0+t35/2.D0 36 : t37 = -0.4495409671875527D-84*t19-0.3918286932266667D4*t22+0.22935 37 : #5086081D2*t-0.458710172162D3 38 : t40 = 0.156710172162D4+3020.D0*t10+0.595089827838D4*t13+0.11467754 39 : #30405D2*t1-0.458710172162D3*t 40 : t42 = 1.D0/2.D0-t35/2.D0 41 : s1 = 31536000.D0 42 : s3 = -0.8765881258141834D0*t26+0.8765881258141834D0*t31+0.58052193 43 : #7625287D-3*(t26/2.D0-t31/2.D0)*t33+0.580521937625287D-3*t36*t37+0. 44 : #3152398523064585D-3*(-t26/2.D0+t31/2.D0)*t40+0.3152398523064585D-3 45 : #*t42*t37 46 : s4 = exp(-0.2514499507768026D2-0.8765881258141834D0*t35+0.58052193 47 : #7625287D-3*t36*t33+0.3152398523064585D-3*t42*t40) 48 : s2 = s3*s4 49 : t50 = s1*s2 50 : FBG2=T50 51 : RETURN 52 : END 53 : 54 :
© Cast3M 2003 - All rights reserved.
Disclaimer