资源描述：
《西安交大fortran上机作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、工程分析程序设计上机作业(三)模块化编程(2)上机目的：进一步掌握内部例程、外部例程、接口块、模块等功能的使用方法。练习例程重载、例程递归的使用方法。1、利用例程重载编写一个子程序求余数ACR(A,B)。要求，对两个整型数和两个实型数都有效(实型数相除的余数：两个实型数相除后的商仍然収一个整数，但剩余值小于除数)[不能用Fortran的内部惭数MOD(x,y)]oMODULEModIMPLICITNONEINTERFACEACRMODULEPROCEDUREACRREALSENDINTERFACECONTAIN
2、SSUBROUTINEACRREALS(A,B,C)IMPLICITNONEREALA,BINTEGERCDOWHILE(A>B)A=A-BENDDOC=AENDSUBROUTINEENDMODULEPROGRAMMAINUSEModIMPLICITNONEREALA,BINTEGERCPRINT弩请输入A,B的值READ*,A,BCALLACR(A,B,C)PRINT*,A和B的余数是:;CENDPROGRAM2、编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中大的那个数作为被除数
3、，两数相除得到一个余数。把余数去除除数得到新一轮的余数。不断重复这一过程直到余数为0,这时的除数就是两个数的最大公约数。调用此函数，求1260,198,72三个数的最大公约数。programmainimplicitnoneintegera,b,c,i,j,kprint*,'请输入要求的数:a,b,c'read*,a,b,ci=gdc(a,b)j=gdc(b,c)k=gdc(i,j)print*,kcontainsfunctiongdc(a,b)implicitnoneintegera,b,c,d,gdcif(b
4、>a)thend=aa=bb=dendifc=mod(a,b)dowhile(c>0)a=bb=cc=mod(a9b)enddogdc=bendfunctionendca*C:Progra>FilesBicrosoftVisualStudioCoB*onISDEV98lYPROJECTS...诂祈入要求的I12601987218Presskeytocontinue3、编写一个函数子程序求sinh(x),以下分三个步骤完成：(1)用递归算法，求出尢"兀Zx，、宀才一心d亠I.X2X3X4xn要求计算—
5、=；(2)自仃编稈求出幺匕1+x11=〉—/?!n-2!3!4!,?=0/?!精度是—<10_6[^能用Fortran的内部函数EXP(x)]；(3)计算sinh(x)=nPROGRAMMAINIMPLICITNONEINTEGERN,XREALMPRINT*；请输入X值「READ*,XM=SINH(X)PRINT*,'SINH(X)二;MCONTAINSFUNCTIONSINH(X)IMPLICITNONEINTEGERI,X,NREALA,B,C,SINH1=1A=X/IDOWHILE(A>1.0E
6、-6)C=l+A1=1+1A=A*(X/I)ENDDOSINH=(C-(l/C))/2ENDFUNCTIONEND4、编写一程序用Euler法求解微分方程^-=y2-x2,当x=0时，y二1.0。试求出dxx=0.1,0.2,0.3,0.4....,1.0时的y值。算法如下：Euler法求解y'=/（兀,y（Q），定解条件：x=x0,y=y（）。取向前差分，令丁=；（儿+1一儿），h=xM+1-xz,h得yw+i=儿+力*/（兀，儿）+。（胪）/(兀0,儿)n/'/(Xl’yjn歹2…nynPROGRAMMAI
7、NIMPLICITNONEREALX,YINTEGERIY=1.0X=0DO1=1,10Y=Y+0」*(Y*Y-X*X)X=X+0.1PRINT*；X=,,X；Y=YENDDOEND
