没有合适的资源？快使用搜索试试~ 我知道了~

文库首页信息化管理其它工程分析程序设计上机作业Fortran部分1.doc

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

1.该资源内容由用户上传，如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款（资源遇到问题，请及时私信上传者）

版权申诉

文档资料

0 下载量 107 浏览量 2022-06-17 05:12:55 上传评论收藏 785KB DOC 举报

温馨提示

试读

64页

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

资源推荐

资源详情

资源评论

工程分析程序设计上机作业（二）

模块化编程（1）

上机目的：掌握内部例程、主程序、外部例程、接口块、模块等功能的使用方法。

1、有一个六边形，求其面积。为求面积，作了 3 条辅助线。如图所示：（提示，三角形面积 =

，其中，a、b、c 为三个边长）。要求用内部函数来实现。

PROGRAM MAIN

IMPLICIT NONE

REAL ::L1=10,L2=20,L3=18,L4=15,L5=21,L6=14,L7=30,L8=36,L9=28

REAL F

F=AREA(L5,L6,L9)+AREA(L1,L8,L9)+AREA(L2,L7,L8)+AREA(L3,L4,L7)

PRINT*,F

CONTAINS

FUNCTION AREA(A,B,C)

REAL S,A,B,C,AREA

S=(A+B+C)/2

AREA=SQRT(S*(S-A)*(S-B)*(S-C))

END FUNCTION

END

2、使用牛顿法求解以下方程的根：（要求用外部子程序）

PROGRAM MAIN

IMPLICIT NONE

REAL F,DF

EXTERNAL GEN

CALL GEN(F,DF)

END

SUBROUTINE GEN(F,DF)

IMPLICIT NONE

REAL ::X=0

REAL F,DF

INTEGER ::N=0

INTEGER ::MaxN=100

REAL ::Eps=1.0E-6

F=X**2+4*X+1

DF=2*X+4

DO WHILE (ABS(F)>Eps.AND.N<MaxN)

X=X-F/DF

N=N+1

F=X**2+4*X+1

DF=2*X+4

END DO

IF (F<=Eps) THEN

PRINT*,'ANSWER IS ',X

ELSE

PRINT*,'NO ANSWER'

END IF

END SUBROUTINE

PROGRAM MAIN

IMPLICIT NONE

REAL F,DF

EXTERNAL GEN

CALL GEN(F,DF)

END

SUBROUTINE GEN(F,DF)

IMPLICIT NONE

REAL ::X=0

REAL F,DF

INTEGER ::N=0

INTEGER ::MaxN=100

REAL ::Eps=1.0E-6

F=7*X**4+6*X**3-5*X**2+4*X+3

DF=28*X**3+18*X**2-10*X+4

DO WHILE (ABS(F)>Eps.AND.N<MaxN)

X=X-F/DF

N=N+1

F=7*X**4+6*X**3-5*X**2+4*X+3

DF=28*X**3+18*X**2-10*X+4

END DO

IF (F<=Eps) THEN

PRINT*,'ANSWER IS ',X

ELSE

PRINT*,'NO ANSWER'

END IF

END SUBROUTINE

3、编写一个子例程子程序，SUM(S，T，N1，N2)。把整型数 N1 到 N2 进行求和，并把求和的结果放

置到 S，把整型数 N1 到 N2 进行求积，并把求积的结果放置到 T。并用这个子程序来计算

PROGRAM MAIN

IMPLICIT NONE

INTEGER N1,N2,S,T,A1,A2,A3,A4,Y

PRINT*,'请输入 N1,N2:'

READ*,N1,N2

CALL SUM(S,T,N1,N2)

A1=S

PRINT*,'请输入 N1,N2:'

READ*,N1,N2

CALL SUM(S,T,N1,N2)

A2=S

PRINT*,'请输入 N1,N2:'

READ*,N1,N2

CALL SUM(S,T,N1,N2)

A3=T

PRINT*,'请输入 N1,N2:'

READ*,N1,N2

CALL SUM(S,T,N1,N2)

A4=T

Y=A1+A2+A3-A4

PRINT*,'Y=',Y

END

SUBROUTINE SUM(S,T,N1,N2)

IMPLICIT NONE

INTEGER S,T,N1,N2,I

S=0

T=1

DO I=N1,N2

S=S+I

T=T*I

END DO

END

4、编写一个模块程序，提供以下服务：定义出常量



、e。定义出子程序，实现求和、求阶乘 n!。

并在主程序中计算如下结果：从键盘上输入整数 n、实型数 A、R、R

，求（实型）和

（实型）

5、 MODULE CX

6、 IMPLICIT NONE

7、 REAL,PARAMETER::PI=3.1415926,E=2.718

8、 CONTAINS

9、 FUNCTION ADD(N)

10、 IMPLICIT NONE

11、INTEGER I,N,ADD

12、 ADD=1

13、 DO I=2,N

14、 ADD=ADD+I

15、 END DO

16、 END FUNCTION

17、 FUNCTION CJ(N)

18、 IMPLICIT NONE

19、 INTEGER I,N,CJ

20、 CJ=1

21、 DO I=1,N

22、 CJ=CJ*I**2

23、 END DO

24、 END FUNCTION

25、 END MODULE

26、 PROGRAM MAIN

27、 USE CX

28、 IMPLICIT NONE

29、 INTEGER N

30、 REAL A,R,R0,D,F

31、 PRINT*,'请输入 n,A,R,R0'

32、 read*,N,A,R,R0

33、 PRINT*,ADD(N)

34、 PRINT*,CJ(N)

35、 D=CJ(N)/ADD(N)

36、 F=A*N/(2*PI*R**2)*(R/R0)**N*E**(-1*(R/R0)**N)

37、 PRINT*,D,F

38、 END

剩余63页未读，继续阅读

评论收藏

内容反馈

1.该资源内容由用户上传，如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款（资源遇到问题，请及时私信上传者）

版权申诉

资源评论

资源反馈

评论星级较低，若资源使用遇到问题可联系上传者，3个工作日内问题未解决可申请退款~

老帽爬新坡

粉丝: 81
资源: 2万+

下载权益

C知道特权

VIP文章

课程特权

VIP享7折，此内容立减4.47元

开通VIP

上传资源快速赚钱

我的内容管理展开

我的资源快来上传第一个资源

我的收益

登录查看自己的收益

我的积分登录查看自己的积分

我的C币登录后查看C币余额

我的收藏

我的下载

下载帮助

前往需求广场，查看用户热搜

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

西安交通大学工程分析程序设计Fortran上机作业参考答案.pdf

哈尔滨工程大学数值分析大作业2014-附fortran程序文件.doc

1999年4月 全国计算机等级考试二级笔试试题 基础部分与Fortran程序设计.doc

FORTRAN模块详解.doc

数值分析大作业-附fortran程序.doc

Fortran基本用法小结.doc

Intel.Visual.Fortran.Composer.XE.2011.5.221.ISO

fortran语言语法.doc

fortran函数总结.doc

Fortran中文教程.pdf

空间桁架结构程序的设计(Fortran).doc

Intel.Visual.Fortran.Composer.XE.2011.5.221.ISO-TBE.iso

2021年计算机技术基础Fortran试题库.doc

Fortran程序设计初步fortran学习课程.pptx

《程序设计基础—fortran90程序实验》.doc

有限单元法FORTRAN程序设计作业.doc

全国计算机等级考试二级Python真题及解析.docx

1000份ppt模版，PPT模板优秀PPT

导入证书可以解决”无法建立到信任根颁发机构的证书链"问题。

matlab批量读取excel表格数据并处理画图

OpenCv车辆识别训练模型

代码随想录知识星球精华-大厂面试八股文第二版v1.2.pdf

数学建模对乙醇偶合制备C4烯烃的问题研究

Vue-Element UI集成ECharts实现数据统计分析页代码部分(如果帮助到你，感谢关注点赞)

STM32F103C8T6中文数据手册

（头歌）计算机组成原理存储系统设计（HUST）1-7关答案

MATLAB深度学习入门实例（果树病虫害识别VGG19版）

（中文）半导体器件物理第三版课后习题答案，施敏著，

网络安全攻防演练防守方方案.docx

最新资源

1999年4月全国计算机等级考试二级笔试试题基础部分与Fortran程序设计.doc