fortran习题集参考答案解析.doc

《fortran习题集参考答案解析.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、工程分析程序设计上机作业（五）模块化编程上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。1、有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面积=，其中，a、b、c为三个边长）。要求用内部函数来计算每个三角形的面积。programmainimplicitnonereal::A(9)=(/10,30,16,13,21,14,20,36,28/)integerirealss=area(A(3),A(4),A(7))+area(A(2),A(7),A(8))+&area(A(1),A(8),A(9))+ar

2、ea(A(5),A(6),A(9))print*,'Theareais',scontainsfunctionarea(x,y,z)implicitnonerealx,y,zrealp,areap=(x+y+z)/2area=sqrt(p*(p-x)*(p-y)*(p-z))endfunctionend1、编写一个子例程子程序，SUM(S，T，N1，N2)。把整型数N1到N2进行求和，并把求和的结果放置到S，把整型数N1到N2进行求积，并把求积的结果放置到T。并用这个子程序来计算：programmainimplicitnoneintegerA(

3、8)integeri,b,c,d,e,f,g,m,n,yb=0;d=0;g=1;n=1doi=1,8A(i)=ienddocallsum(b,c,1,4)callsum(d,e,3,8)callsum(f,g,3,6)callsum(m,n,1,3)y=b+d+g-nprint*,'y=(1+2+3+4)+(3+4+5+6+7+8)+(3*4*5*6)-(1*2*3)=',ycontainssubroutinesum(s,t,N1,N2)integers,t,N1,N2,idoi=N1,N2s=s+A(i)t=t*A(i)enddoendsu

4、broutineend2、编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。调用此函数，求1260，198，72三个数的最大公约数。programmainimplicitnoneintegera,b,c,m,n,tempprint*,'请输入三个数：'print*,'abc'read*,a,b,cif(a

5、,m)if(c

6、itnonereal,parameter::PI=3.1415926real,parameter::e=2.7182818containsfunctionsub(n)integeri,n,subsub=0doi=1,nsub=sub+i*ienddoendfunctionsubfunctionfact(n)integeri,nrealfactfact=1doi=1,nfact=fact*ienddoendfunctionfactendmodulemytaskprogrammainusemytaskimplicitnoneintegernrea

7、lA,R,R0realx,yprint*,'请输入依次n,A,R,R0:'read*,n,A,R,R0x=fact(n)/sub(n)y=(A*n/(2*PI*R*R))*(R/R0)**n*e**(-(R/R0)**n)print*,'n!/∑n^2(i=1,n)=',xprint*,'(A*n/(2*PI*R*R))*(R/R0)**n*e**(-(R/R0)**n)',yendprogram1、编写外部函数，计算的值，直到最后一项的绝对值小于10-6为止。再编写主程序，从键盘读入x，调用该外部函数并输出sinx的计算结果。注意不能用Fo

8、rtran的标准函数SIN(X)。programmainimplicitnonerealx,pexternalsubprint*,'请输入角度x:'read*,xp