西安交大fortran上机作业.docxVIP

工程分析程序设计上机作业（三）模块化编程（2）上机目的：进一步掌握内部例程、外部例程、接口块、模块等功能的使用方法。练习例程重载、例程递归的使用方法。利用例程重载编写一个子程序求余数ACR(A,B)。要求，对两个整型数和两个实型数都有效（实型数相除的余数：两个实型数相除后的商仍然取一个整数，但剩余值小于除数）[不能用Fortran的内部函数MOD(x,y)]。MODULE Mod IMPLICIT NONE INTERFACE ACR MODULE PROCEDURE ACRREALSEND INTERFACECONTAINS SUBROUTINE ACRREALS(A,B,C) IMPLICIT NONE REAL A,B INTEGER C DO WHILE (AB) A=A-B END DO C=A END SUBROUTINEEND MODULEPROGRAM MAIN USE Mod IMPLICIT NONE REAL A,B INTEGER C PRINT*,请输入A,B的值: READ*,A,B CALL ACR(A,B,C) PRINT*,A和B的余数是:,C END PROGRAM编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。调用此函数，求1260，198，72三个数的最大公约数。program mainimplicit noneintegera,b,c,i,j,k print*,请输入要求的数:a,b,cread*,a,b,ci=gdc(a,b) j=gdc(b,c) k=gdc(i,j)print*,kcontainsfunctiongdc(a,b)implicit noneintegera,b,c,d,gdcif (ba) then d=a a=b b=dend if c=mod(a,b)do while (c0) a=b b=c c=mod(a,b)end dogdc=bend functionend编写一个函数子程序求sinh(x)，以下分三个步骤完成：（1）用递归算法，求出；（2）自行编程求出，要求计算精度是[不能用Fortran的内部函数EXP(x)]；（3）计算。 PROGRAM MAINIMPLICIT NONEINTEGER N,XREAL M PRINT*,请输入X值: READ*,X M=SINH(X) PRINT*,SINH(X)=,M CONTAINS FUNCTION SINH(X) IMPLICIT NONE INTEGER I,X,N REAL A,B,C,SINH I=1 A=X/I DO WHILE (A1.0E-6) C=1+A I=I+1 A=A*(X/I) END DO SINH=(C-(1/C))/2 END FUNCTION END编写一程序用Euler法求解微分方程，当x=0时，y=1.0。试求出x=0.1,0.2,0.3,0.4….,1.0时的y值。算法如下：Euler法求解，定解条件：。取向前差分，令，得PROGRAM MAINIMPLICIT NONEREAL X,YINTEGER I Y=1.0 X=0 DO I=1,10 Y=Y+0.1*(Y*Y-X*X) X=X+0.1 PRINT*,X=,X,Y=,Y END DOEND