工程分析程序设计 上机作业(二).docx

工程分析程序设计 上机作业（二） 模块化编程（1） 上机目的：掌握内部例程、主程序、外部例程、接口块、模块等功能的使用方法。 有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面积=，其中，a、b、c为三个边长）。要求用内部函数来实现。 PROGRAM MAIN REAL::L1=10,L2=20,L3=18,L4=15,L5=21,L6=14,L7=30,L8=36,L9=28 PRINT*,S(L5,L6,L9)+S(L1,L8,L9)+S(L2,L8,L7)+S(L3,L7,L4) CONTAINS FUNCTION S(A,B,C) REAL P,A,B,C,S P=(A+B+C)/2 S=SQRT(P*(P-A)*(P-B)*(P-C)) END FUNCTION S END 使用牛顿法求解以下方程的根：（要求用外部子程序） (1) PROGRAM Main IMPLICIT NONE REAL A,B,C COMPLEX X EXTERNAL NEWTON PRINT*,请输入一元二次方程的系数A,B,C: READ*,A,B,C CALL NEWTON(A,B,C,X) PRINT*,X=,X END PROGRAM SUBROUTINE NEWTON(A,B,C,X) IMPLICIT NONE REAL A,B,C COMPLEX X REAL F,F1 X=0.0 F=A*X**2+B*X+C F1=2*A*X+B DO WHILE (ABS(F)=1E-6) F1=2*A*X+B F=A*X**2+B*X+C X=X-F/(F1) END DO END SUBROUTINE (2) PROGRAM Main IMPLICIT NONE REAL A,B,C,D,E COMPLEX X EXTERNAL NEWTON PRINT*,请输入一元四次方程的系数A,B,C,D,E: READ*,A,B,C,D,E CALL NEWTON(A,B,C,D,E,X) PRINT*,X=,X END PROGRAM SUBROUTINE NEWTON(A,B,C,D,E,X) IMPLICIT NONE REAL A,B,C,D,E COMPLEX X REAL F,F1 X=0.0 F=A*X**4+B*X**3+C*X**2+D*X+E F1=4*A*X**3+3*B*X**2+2*C+D DO WHILE (ABS(F)=1E-6) F=A*X**4+B*X**3+C*X**2+D*X+E F1=4*A*X**3+3*B*X**2+2*C+D X=X-F/(F1) END DO END SUBROUTINE 编写一个子例程子程序，SUM(S，T，N1，N2)。把整型数N1到N2进行求和，并把求和的结果放置到S，把整型数N1到N2进行求积，并把求积的结果放置到T。并用这个子程序来计算 PROGRAM MAIN IMPLICIT NONE INTEGER N1,N2 INTEGER::M1=1,M2=4,M3=3,M4=8,M5=3,M6=6,M7=1,M8=3 PRINT*,y=,S(M1,M2)+S(M3,M4)+T(M5,M6)-T(M7,M8) CONTAINS FUNCTION S(N1,N2) IMPLICIT NONE INTEGER S,N1,N2 S=N1 DO WHILE (N1N2) N1=N1+1 S=S+N1 END DO END FUNCTION FUNCTION T(N1,N2) IMPLICIT NONE INTEGER T,N1,N2 T=N1 DO WHILE (N1N2) N1=N1+1 T=T*N1 END DO END FUNCTION END PROGRAM MAIN 编写一个模块程序，提供以下服务：定义出常量?、e。定义出子程序，实现求和、求阶乘n!。并??主程序中计算如下结果：从键盘上输入整数n、实型数A、R、R0，求（实型）和（实型） MODULE Myutils IMPLICIT NONE REAL,PARAMETER::PI=3.1415926,e=2CONTAINS FUNCTION X(n) INTEGER n,X,I I=1 DO WHILE (I=n) X=X+I**2 I=I+