工程分析程序设计上机作业Fortran部分1.docVIP

工程分析程序设计 上机作业（二） 模块化编程（1） 上机目的：掌握内部例程、主程序、外部例程、接口块、模块等功能的使用方法。 有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面积=，其中，a、b、c为三个边长）。要求用内部函数来实现。 PROGRAM MAIN IMPLICIT NONE REAL ::L1=10,L2=20,L3=18,L4=15,L5=21,L6=14,L7=30,L8=36,L9=28 REAL F F=AREA(L5,L6,L9)+AREA(L1,L8,L9)+AREA(L2,L7,L8)+AREA(L3,L4,L7) PRINT*,F CONTAINS FUNCTION AREA(A,B,C) REAL S,A,B,C,AREA S=(A+B+C)/2 AREA=SQRT(S*(S-A)*(S-B)*(S-C)) END FUNCTION END 使用牛顿法求解以下方程的根：（要求用外部子程序） PROGRAM MAIN IMPLICIT NONE REAL F,DF EXTERNAL GEN CALL GEN(F,DF) END SUBROUTINE GEN(F,DF) IMPLICIT NONE REAL ::X=0 REAL F,DF INTEGER ::N=0 INTEGER ::MaxN=100