fortran习题参考答案

工程分析程序设计 上机作业（五）模块化编程上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。1、 有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面积=，其中，a、b、c为三个边长）。要求用内部函数来计算每个三角形的面积。program main implicit none real : A(9)=(/10,30,16,13,21,14,20,36,28/) integer i real s s=area(A(3),A(4),A(7)+area(A(2),A(7),A(8)+ & area(A(1),A(8),A(9)+area(A(5),A(6),A(9) print*,The area is,s contains function area(x,y,z) implicit none real x,y,z real p,area p=(x+y+z)/2 area=sqrt(p*(p-x)*(p-y)*(p-z) end functionend2、 编写一个子例程子程序，SUM(S，T，N1，N2)。把整型数N1到N2进行求和，并把求和的结果放置到S，把整型数N1到N2进行求积，并把求积的结果放置到T。并用这个子程序来计算：program main implicit none integer A(8) integer i,b,c,d,e,f,g,m,n,y b=0;d=0;g=1;n=1 do i=1,8 A(i)=iend do call sum(b,c,1,4) call sum(d,e,3,8) call sum(f,g,3,6) call sum(m,n,1,3) y=b+d+g-n print*,y=(1+2+3+4)+(3+4+5+6+7+8)+(3*4*5*6)-(1*2*3)=,y contains subroutine sum(s,t,N1,N2) integer s,t,N1,N2,i do i=N1,N2 s=s+A(i) t=t*A(i) end do end subroutine end 3、 编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。调用此函数，求1260，198，72三个数的最大公约数。program main implicit none integer a,b,c,m,n,temp print*,请输入三个数： print*,a b c read*,a,b,c if (ab) then temp=a a=b b=temp end if call GDC(a,b,m) if (c1.0e-6) y=y+(-1)*(i-1)/2)*k*i/(fact(i) i=i+2 end do print*,sin(x) = ,y contains function fact(m) integer m,j,factfact=1do j=1,m fact=fact*jend do end function factend subroutine sub6、 编写程序，用欧拉法数值求解一阶微分方程。数值解法的基本思想：用差分方程代替微分方程，然后在若干个离散点上逐点求解差分方程，得到各离散点x0, x1, x2, 等处函数的近似值y0, y1, y2, 。其中各离散点x0, x1, x2, 之间的距离称为步长h。欧拉法的基本思想是，一阶微分可用向前差分代替，即。带入微分方程，可得因此，其中。当给定x0, y0和步长h后，即可按下面步骤求得数值解：请按上述方法编写程序用Euler法求解微分方程，当x0=0时，y0=1.0。取h=0.1，试求出当xi=0.1,0.2,0.3,0.4.,1.0时的yi值。program main implicit none real x(0:10) real y(0:10) real p(0:10) integer(2) i,j y(0)=1.0; x(0)=0.0 p(0)=f(x(0),y(0) do i=1,10 x(i)=x(i-1)+0.1 end do do j=1,10 y(j)=y(j-1)+0.1*p(j-1) p(j)=f(x(j),y(j) prin