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

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

ID:50960769

大小:644.50 KB

页数:41页

时间:2020-03-08

工程分析程序设计上机作业Fortran部分.doc_第1页
工程分析程序设计上机作业Fortran部分.doc_第2页
工程分析程序设计上机作业Fortran部分.doc_第3页
工程分析程序设计上机作业Fortran部分.doc_第4页
工程分析程序设计上机作业Fortran部分.doc_第5页
资源描述:

《工程分析程序设计上机作业Fortran部分.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、工程分析程序设计上机作业(二)模块化编程(1)上机目的:掌握内部例程、主程序、外部例程、接口块、模块等功能的使用方法。1、有一个六边形,求其面积。为求面积,作了3条辅助线。如图所示:(提示,三角形面积=,其中,a、b、c为三个边长)。要求用内部函数来实现。PROGRAMMAINIMPLICITNONEREAL::L1=10,L2=20,L3=18,L4=15,L5=21,L6=14,L7=30,L8=36,L9=28REALFF=AREA(L5,L6,L9)+AREA(L1,L8,L9)+AREA(L2

2、,L7,L8)+AREA(L3,L4,L7)PRINT*,FCONTAINSFUNCTIONAREA(A,B,C)REALS,A,B,C,AREAS=(A+B+C)/2AREA=SQRT(S*(S-A)*(S-B)*(S-C))ENDFUNCTIONEND2、使用牛顿法求解以下方程的根:(要求用外部子程序)PROGRAMMAINIMPLICITNONEREALF,DFEXTERNALGENCALLGEN(F,DF)ENDSUBROUTINEGEN(F,DF)IMPLICITNONEREAL::X=0RE

3、ALF,DFINTEGER::N=0INTEGER::MaxN=100REAL::Eps=1.0E-6F=X**2+4*X+1DF=2*X+4DOWHILE(ABS(F)>Eps.AND.N

4、ALLGEN(F,DF)ENDSUBROUTINEGEN(F,DF)IMPLICITNONEREAL::X=0REALF,DFINTEGER::N=0INTEGER::MaxN=100REAL::Eps=1.0E-6F=7*X**4+6*X**3-5*X**2+4*X+3DF=28*X**3+18*X**2-10*X+4DOWHILE(ABS(F)>Eps.AND.N

5、X+4ENDDOIF(F<=Eps)THENPRINT*,'ANSWERIS',XELSEPRINT*,'NOANSWER'ENDIFENDSUBROUTINE1、编写一个子例程子程序,SUM(S,T,N1,N2)。把整型数N1到N2进行求和,并把求和的结果放置到S,把整型数N1到N2进行求积,并把求积的结果放置到T。并用这个子程序来计算PROGRAMMAINIMPLICITNONEINTEGERN1,N2,S,T,A1,A2,A3,A4,YPRINT*,'请输入N1,N2:'READ*,N1,N2CA

6、LLSUM(S,T,N1,N2)A1=SPRINT*,'请输入N1,N2:'READ*,N1,N2CALLSUM(S,T,N1,N2)A2=SPRINT*,'请输入N1,N2:'READ*,N1,N2CALLSUM(S,T,N1,N2)A3=TPRINT*,'请输入N1,N2:'READ*,N1,N2CALLSUM(S,T,N1,N2)A4=TY=A1+A2+A3-A4PRINT*,'Y=',YENDSUBROUTINESUM(S,T,N1,N2)IMPLICITNONEINTEGERS,T,N1,N2

7、,IS=0T=1DOI=N1,N2S=S+IT=T*IENDDOEND1、编写一个模块程序,提供以下服务:定义出常量p、e。定义出子程序,实现求和、求阶乘n!。并在主程序中计算如下结果:从键盘上输入整数n、实型数A、R、R0,求(实型)和(实型)2、MODULECX3、IMPLICITNONE4、REAL,PARAMETER::PI=3.1415926,E=2.7185、CONTAINS6、FUNCTIONADD(N)7、IMPLICITNONE8、INTEGERI,N,ADD9、ADD=110、DOI

8、=2,N11、ADD=ADD+I12、ENDDO13、ENDFUNCTION14、FUNCTIONCJ(N)15、IMPLICITNONE16、INTEGERI,N,CJ17、CJ=118、DOI=1,N19、CJ=CJ*I**220、ENDDO21、ENDFUNCTION22、ENDMODULE23、PROGRAMMAIN24、USECX1、IMPLICITNONE2、INTEGERN3、REALA,R,R0,D,F4、PRINT*,'请

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。