R=DSQRT(1.DO-QE/QT) Regression equation is (5分) F=N-M-】)LU/OEM FORMAT(/4X,'Roots of mormal equation ,5X,'A(i,m+1),i=1,m'& T1gp=2077. +6.081 *T+110.2 /1x,601H- FORMAT(1X,10G12.4) R= .999999 RETURN Fs 686329. END 6.08138 B= 1407.22 SUBROUTINE GS(A.N.M,EPS) (5分) = -110.160 (10分) DOUBLE PRECISION EPS,A DIMENSION A(6,7) 1.计算原理（化学原理和计算方法）(15分) DO 50K=LN 放射性衰变系列 P=0.0 D0101=K,N 2U4→20Th点→2Ra—→… IFOPGE ABS(AOLKINGOTO I0 P=ABS(A(LK)) 衰变方式服从连串反应动力学， L=I 其动力学方程式为一个常微分方程组： CONTINUE IF(P.LT.EPS)STOPNO SOLUTION IF(L.EQ.K)GOTO 30 d凹=-k,U叮 DO20 EK M (10分) TMP=A(L,J) A(LJA(KJ) d=k心]-kT A(K JETMP CONTINUE P=L0A(KK) 利用RngK法解此方程组。可符2过U和20Th这两种核素在10a内每释2a DO40 K+LM A(KJ)=A(K.J)*P 的质量变化 (5分) DO40 BK+IN A(1,JFA(IJ)-A(L,KA(K,J) 2.程序框图(20分》 CONTINUE 50 CONTINUE 开始 DOMK=2N I=M-K DO 60 J=I+1,N 输入：[U,Th初值234,230 A(L,M)=A(L,M)-A(LJ)*A(J,M) 时间间隔H, 计算精度PS (5分) 60 CONTINUE RETURN 计算：积分步长H0(5分) END 4.运行结果。(15分) 调用变步长Runge--Kuta法子程序解动力学方程组计算时刻t的U, [Th值 (其中调用计算动力学方程组的子程序计算DY值)(5分) 3 R=DSQRT(1.D0-QE/QT) F=(N-M-1)*U/QE/M 22 FORMAT(/4X,'Roots of mormal equation',5X,'A(i,m+1),i=1,m' & /1x,60(1H-)) 33 FORMAT(1X,10G12.4) RETURN END SUBROUTINE GS(A,N,M,EPS) （5 分） DOUBLE PRECISION EPS,A DIMENSION A(6,7) DO 50 K=1,N P=0.0 DO 10 I=K,N IF(P.GE.ABS(A(I,K)))GOTO 10 P=ABS(A(I,K)) L=I 10 CONTINUE IF(P.LT.EPS)STOP'NO SOLUTION' IF(L.EQ.K)GOTO 30 DO 20 J=K,M TMP=A(L,J) A(L,J)=A(K,J) A(K,J)=TMP 20 CONTINUE 30 P=1.0/A(K,K) DO 40 J=K+1,M A(K,J)=A(K,J)*P DO 40 I=K+1,N A(I,J)=A(I,J)-A(I,K)*A(K,J) 40 CONTINUE 50 CONTINUE DO 60 K=2,N I=M-K DO 60 J=I+1,N A(I,M)=A(I,M)-A(I,J)*A(J,M) 60 CONTINUE RETURN END 4.运行结果。（15 分） Regression equation is: （5 分） ------------------------------------------------------------ Tlgp=-2077. + 6.081 *T+ 110.2 *lgp R= .999999 F= 686329. A= 6.08138 B= 1407.22 C= -110.160 （10 分） B 1.计算原理（化学原理和计算方法）（15 分） 放射性衰变系列 234 92 U⎯⎯→k1 23090Th⎯⎯→k2 22684 Ra ⎯⎯→…… 衰变方式服从连串反应动力学。 其动力学方程式为一个常微分方程组： [U] [Th] d d[Th] [U] dd[U] 1 2 1 k k t k t = − = − （10 分） 利用 Runge-Kutta 法解此方程组，可得 U 23492 和 Th 23090 这两种核素在 10 a 内每隔 2 a 的质量变化。 （5 分） 2.程序框图（20 分） 开始 输入：[U]，[Th]初值 234，230 时间间隔 H， 计算精度 EPS （5 分） 计算：积分步长 H0 （5 分） 调用变步长Runge-Kutta法子程序解动力学方程组计算时刻t的[U]， [Th]值 （其中调用计算动力学方程组的子程序计算 DY 值） （5 分）