Fortran上机作业(六)

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

数组（2）

上机目的：练习数组的声明、存储、操作，以及数组参数、动态数组、数组函数的使用。

1、从键盘上输入10个整数，然后逆序打印出来。

program array

implicit none

integer,parameter::n=10

integer i

real A(n),temp

read*,A

DO i=1,n/2

temp=A(i)

A(i)=A(n+1-i)

A(n+1-i)=temp

end DO

print*,A

end program

2、输入任意n个数存放在数组中（如5个数1、2、8、2、10），请在屏幕上打印如下方阵

1 2 8 2 10

10 1 2 8 2

2 10 1 2 8

8 2 10 1 2

2 8 2 10 1

program array

implicit none

integer,parameter::n=5

integer i,j

real A(n),temp

print*,'input',n,'numbers'

read*,A

print*,'output the array'

print*,A

DO i=1,n-1

temp=A(n)

DO j=n,2,-1

A(j)=A(j-1)

end Do

A(1)=temp

print*,A

end DO

end program

3、用选择法进行排序：假定A数组中有n个数，设定两个变量P和J，J按顺序从数组开

头指向数组的各个位置（也就是应该存放数据的位置），P则指向剩余的数据中最小的数。

然后把J指向的元素和P指向的元素进行比较大小，若P指向的元素更小，则把二者对换。如此把J从1号位向后移动至数组末尾，则完成了排序。例如

第一步：

15 8 4 10 1

J P J指向1号位，P指向2~5号位的最小值5号位的1，二者对换。

第二步：

1 8 4 10 15

J P J指向２号位，P指向3~5号位的最小值3号位的4，二者对换。

第三步：

1 4 8 10 15

J P J指向3号位，P指向4~5号位的最小值4号位的10，二者不对换。第四步：

1 4 8 10 15

J P J指向4号位，P指向5号位的最小值号位的15，二者不对换。

根据以上四步，整个数组完成排序。要求在屏幕上显示排序法每步进行的情况。有兴趣的同学还可以比较选择法排序和“冒泡算法”进行排序的速度。

!选择法

program array

implicit none

integer,parameter::n=10 !数列元素个数

integer i,j

real A(n),temp

print*,'input',n,'numbers'

read*,A

print*,'output the array'

DO i=1,n DO j=i+1,n if (A(i)>A(j)) then temp=A(i) A(i)=A(j) A(j)=temp endif end Do enddo

print*,A end program

4、 高斯消去法求解线性代数方程组：对于一般的n 阶方程组，

??????

?=+++=+++=+++n

n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a 22112222212111212111.....................................................

高斯消去法步骤如下：第一步：若011≠a ，令n i a a l i i ,...,3,2,/1111==，用()1i l -乘第1个方程加到第i 个方程上()n i ,...,3,2=，得同解方程组

(1)(1)(1)(1)11112211(2)(2)(2)

22222(2)(2)(2)22...................................n n n n n nn n n a x a x a x b a x a x b a x a x b ?+++=?++=??

??++=?

其中 .,,3,2,,,,3,2,,,,,,3,2,1,11)

2(11)2(1)1(11)1(1n i b l b b n j i a l a a b b n j a a i i i j i ij ij j j =-==-==== 第二步：若,0)2(22≠a 令n i a a l i i ,,4,3,/)

2(22)2(22 ==，用)(2i l -乘第2个方程加到第i 个方程上

),,4,3(n i =，则将),,4,3()

2(2

n i a i =消去。一般，设第k-1步后方程组化为如下的同解方程组 (1)(1)(1)(1)

11112211(2)(2)(2)

22222

()()

()()()()

............................................................................n n n n k k k kk k kn n k k k k nk n nn n n a x a x a x b a x a x b a x a x b a x a x b ?+++=?++=????++=?

??

?++=?

则第k 步：若若,0)(≠k kk a 令n k k i a a l k kk k ik

ik ,,2,1,/)()( ++==，用)(ik l -乘第k 个方程加到第i 个方程上(,,,2,1n k k i ++=)，得到如下的同解方程组

????

??

???????=++=++=++=++=++++++++++++++++)

1()1(1)1(1,)

1(1)1(,11)1(1,1)()

()()

2(2

)2(22)2(22)1(1)1(12)1(121)1(11............................................................................k n n k nn k k k n k k n k n k k k k k k k n k kn k k kk n n n n b x a x a b x a x a b x a x a b x a x a b x a x a x a 其中(1)()()(1)()(),.1,,,,1,

,k k k ij ij ik kj k k k i

i ik k a a l a i j k n b b l b i k n

++?=-=+??=-=+?? 按上述做法，做完n-1步，原方程组化为同解的上三角形方程组

??

??

?

??????==++=++=+++)

()()()

()()

2(2

)2(22)2(22)

1(1)1(12)1(121)1(11....................................................................n n n n nn k k n k kn k k kk n n n n b x a b x a x a b x a x a b x a x a x a 最后，设0)

(≠n nn

a ，逐步代回得原方程组的解()()

()()1/,()/(1,2,,2,1).n n n n nn n

k k k

k k kj j kk

j k x b a x b a x a k n n =+?=??=-??

?=--?

∑

注意：上述公式中)()

(,k i k ij b a 的上标k ，是用来区别消去过程中第k 步利用的量。在用编程求解

时，可把)(k ij a 存在ij a 位置，)(k i b 存在i b 位置。

解方程组（1）?????-=++-=-+=++3

118743513213213

21321x x x x x x x x x （2）???=+=+-232121021215

x x x x program Gauss implicit none

integer,parameter::n=3 !数列每维元素个数 integer i,j,k,f

real ::A(n,n+1),temp=0,X(n)

!data ((A(i,j),i=1,3),j=1,4) /5,1,7,1,2,8,-3,3,11,-4,1,-3/

print*,'input',n,'*',n+1,'数列，注意：按列输入' read*,A

DO i=1,n-1 DO j=i+1,n if (A(i,i)/=0) then

A(j,:)=-A(j,i)/A(i,i)*A(i,:)+A(j,:) !right

end if enddo

enddo

此方程病态，提示：列主元消去法，第k 步

消去过程选取第k 行绝对值最大元素a kq ，交换k 和q 列，然后继续消去过程

X(n)=A(n,n+1)/A(n,n)

DO k=n-1,1,-1

DO j=k+1,n

temp=temp+A(k,j)*x(j)

enddo

X(k)=(A(k,n+1)-temp)/A(k,k)

temp=0

enddo

Do i=1,n

print*,'X',i,'= ',x(i)

enddo

end program

FORTRAN 90 程序编程规范

FORTRAN 90 程序编程规范 Fortran 90 编程规范，使程序代码高度组织化，更加易读、易懂、易于维护，程序更加高效。使编出的程序更易懂、易于维护。 1 语言选择 数值预报创新系统软件开发应避免使用Fortran77 的某些过时特征以Fortran 90不一致的特征。选择Fortran 90 作为开发语言，并采用Fortran 90 的新功能，如动态内存的分配(dynamic memory allocation)、递归（recursion ）, 模块(modules)、POINTER 、长变量名、自由格式等。 Fortran 77其中某些只是一些冗余的功能，这些功能已经过时，另外，还有一些在Fortran90 中被证明是不好的用法，建议不要使用。 2 Fortran 90 的新特性 2.1.1 建议使用的Fortran 90 新特性 建议使用Fortran 90 提供的模块（module ），并用Use ONLY 指定module 中哪些变量或派生类型定义可用于调用程序。 尽量使用数组下标三元组，这样可优化并减少所需的代码行数。为提高可读性，要在括号内表明数组的维数，例如： 1dArrayA(:) = 1dArrayB(:) + 1dArrayC(:) 2dArray(: , :) = scalar * Another2dArray(: , :) 当访问数组的子集时，例如在有限差分等式中，可以通过使用下标三元组实现。例如：2dArray(: , 2:len2) = scalar *（ & Another2dArray(:, 1:len2 -1) & - Another2dArray(:, 2:len2) & ) 对程序单元（program units ）命名，并使用End program ，End subroutine ，End interface ，End module 等结构再次指定“program unit ”的名称。 在逻辑表达式中使用>、 >=、 ==、 <、 <=、 /=，它们分别代 替.gt.、.ge.、.eq.、.lt.、.le.、.ne. 。新的表示方法更接近标准的数学符号 在变量定义中始终使用“：：”；始终用“DIMENSION ”定义数组形状；始终用（len=）的语法格式声明字符变量的长度。

Fortran95习题答案

第四章 1.program main implicit none write(*,*) "Have a good time." write(*,*) "That's not bad." write(*,*) '"Mary" isn''t my name.' end program 2.program main real, parameter :: PI=3 implicit none.14159 real radius write(*,*) "请输入半径长" read(*,*) radius write(*,"(' 面积='f8. 3)") radius*radius*PI end program 3.program main implicit none real grades write(*,*) "请输入成绩" read(*,*) grades write(*,"(' 调整后成绩为'f8.3)") SQRT(grades)*10.0 end program 4.integer a,b real ra,rb a=2 b=3 ra=2.0 rb=3.0 write(*,*) b/a ! 输出1, 因为使用整数计算, 小数部分会无条件舍去write(*,*) rb/ra ! 输出1.5 5.program main implicit none type distance real meter, inch, cm end type type(distance) :: d write(*,*) "请输入长度:" read(*,*) d%meter d%cm = d%meter*100 d%inch = d%cm/2.54 write(*,"(f8.3'米='f8.3'厘米='f8.3'英寸')") d%meter, d%cm, d%inch end program 第五章 1.program main

fortran常见问题解决

楼主为了减少重复回答问题，特编此帖，并不定期添加和更新内容。 错误难免，欢迎讨论，仅供参考。 很多人问哪里可以找到Fortran编译器，有不少热心学友提供网址，特汇集在这里。虽然俺检验过这些链接，但是它们不一定总有效。 Fortran编译器下载： CVF？ FTN95(License:Freeforpersonaluse) 以下操作，如无特别说明，都是以为例。 1.如何加大Stacksize？ 选Project=>Settings=>Link=>Category:Output=>？ Stackallocations Reserve:这里填新值（默认为1M，若需要10M，则填） 2.如何用Fortran批量生成文件？ 设要生成4000个文件，文件名为AA1-AA4000，如何写循环生成文件,而不用写4000次write 命令呢？ 用内部文件： character(len=80)::filename,form integer::i doi=1,4000 selectcase(i) case(1:9) write(form,'(i1)')i case(10:99) write(form,'(i2)')i case(100:999) write(form,'(i3)')i case(1000:9999) write(form,'(i4)')i endselect write(filename,*)"AA",trim(form),".TXT" open(10,file=filename) write(10,*)i close(10)

enddo？ stop end 3.如何用Fortran动态生成输出格式？ 设有一个数组data(100)，输出时，希望每行输出num个数，而num由用户输入，如何实现？ 用内部文件： character(len=80)::form real::data(100) integer::i,num data=(/(i,i=1,100)/)/ read(*,*)num write(form,*)"(",num,"" write(*,form)data stop end 4.MS是不是很垃圾？ 是垃圾，其中Bug太多，多到不可用的地步！ 在这个主题里，换了CVF后问题就没了的人已有相当的数目。 如果你用，遇到莫名其妙的错误，建议换，这是一个比较成熟的编译器。 5.如何用F90/95生成随机数？ 注意： 现在计算机产生的随机数都是伪随机数。 random_number(x)产生一个0到1之间的随机数（x可以是向量），但是每次总是那几个数。用了random_seed()后，系统根据日期和时间随机地提供种子，使得随机数更随机了。 programrandom implicitnone real::x callrandom_seed()!系统根据日期和时间随机地提供种子 callrandom_number(x)!每次的随机数就都不一样了 write(*,*)x stop endprogramrandom 6.函数/子程序超载的例子

fortran程序案例题汇编(14道)

1．Fibonacci数列定义如下： F1=1 F2=1 F n=F n-1+F n-2 (n>2) 求Fibonacci数列的前30项。 integer f(30),i f(1)=1 f(2)=2 do i=3,30 f(i)=f(i-1)+f(i-2) enddo print*,f end 2．输入10个学生的总分，求每个学生的名次integer s(10),a(10),i,j do i=1,10 read*,s(i) enddo do i=1,10 a(i)=1 do j=1,10 if(s(i)

3．给定一组数，按照从小到大的顺序输出。integer a(10) integer p do i=1,10 read *,a(i) enddo do j=1,9 p=j do i=j+1,10 if (a(i)

4．输入若干名学生的学号和三门课程（语数英）的成绩，要求从键盘输入一个学生的学号，能打印出该学生的三门功课成绩和总分。 character*6,dimension(:),allocatable::xue integer,dimension(:,:),allocatable::g integer,dimension(:),allocatable::zong integer i,j,n character*6,xh print *,"请输入学生的个数" read *,n allocate(xue(n)) allocate(g(n,3)) allocate(zong(n)) do i=1,n read *,xue(i),(g(i,j),j=1,3) enddo do i=1,n zong(i)=0 do j=1,3 zong(i)=zong(i)+g(i,j) enddo enddo print *,"请输入你要打印的学生的学号" read *,xh do i=1,n if(xue(i)==xh)then print *,(g(i,j),j=1,3),zong(i) exit endif enddo end

几种Fortran 编译器简介

几种Fortran 编译器 --------------------------------------------------------------------------------------- 1.CVF Compaq Visual Fortran （CVF）, 当今PC平台上功能相当强大与完整的Fortran程序开发工具，还用于Abaqus的开发。 1997年，微软将Fortran PowerStation卖给DEC之后，微软就不再出版Fortran编译器了。后来DEC并入了Compaq，再后来Compaq又和HP合并了。现在最新的版本是HP出的Fortran for Windows v6.6，现在HP/Compaq已经不再开发Fortran了，CVF 6.6是最终的版本了，Compaq的Fortran开发小组已经投入Intel旗下，目前Intel已经有Intel Visual Fortran 11.0。Compaq Visual Fortran 6.6官方的单价也相当昂贵。 Compaq Visual Fortran 6.6 下载： https://www.wendangku.net/doc/e010156287.html,/SoftDown.asp?ID=11937 Compaq Visual Fortran 6.6 绿色版下载： https://www.wendangku.net/doc/e010156287.html,/down/10915.html Compaq Visual Fortran 6.5 下载： https://www.wendangku.net/doc/e010156287.html,/soft/fortran6.5.rar ftp://2006:2006@https://www.wendangku.net/doc/e010156287.html,/36/https://www.wendangku.net/doc/e010156287.html,-002124.rar --------------------------------------------------------------------------------------- 2. IVF Intel Visual Fortran (IVF)将Compaq Visual Fortran* (CVF) 语言的丰富功能与英特尔代码生成及优化技术结合在一起。目前Intel已经有Intel Visual Fortran 11.0。 下载： Intel官方网站 ftp://166.111.26.159/software/science/ ftp://202.112.85.101/pub2/Windows/Scientific_Tools/Fortran/https://www.wendangku.net/doc/e010156287.html,p iler/ 从https://www.wendangku.net/doc/e010156287.html,/maths/下载。IVF10可以用IVF9的licence https://www.wendangku.net/doc/e010156287.html,/irc_nas/730/w_cc_p_10.0.025_ia3 2.exe 集成VS2005的IVF10.1：w_fc_p_10.1.021.exe （302MB） Intel(R) Visual Fortran for IA-32 and Intel(R) 64 (with Microsoft Visual Studio 2005

fortran程序实例

1)实例3—求多个半径下的圆周长 ! z3.f90 --Fortran95 ! FUNCTIONS: ! z3 - Entry point of console application. !************************************************************************* ! PROGRAM: z3 ! PURPOSE: Entry point for the console application. !************************************************************************ program z3 ! 求多个半径下的圆周长 ! 主程序 ! PROGRAM Z3 PRINT *, 'R=',1.2,'C=',C(1.2) PRINT *, 'R=',3.4,'C=',C(3.4) PRINT *, 'R=',15.6,'C=',C(15.6) PRINT *, 'R=',567.3,'C=',C(567.3) END program z3 !子程序 FUNCTION C(R) PI=3.1415926 C=2*PI*R RETURN ! Body of z3 end

2)实例4—键盘与显示器输入/输出 a)Fortran 基本操作 b)程序指令 ! ZXZ_I_O.f90 ! FUNCTIONS: ! ZXZ_I_O - Entry point of console application. ! PROGRAM: ZXZ_I_O ! PURPOSE: Entry point for the console application. !***************输入、输出样式种种************************** program ZXZ_I_O implicit none !变量声明的位置 INTEGER(2) i; INTEGER(4) j; INTEGER(4) m; REAL n INTEGER A,B ! Variables PRINT*,'输入整数A'; READ*, A PRINT*,'输入整数B'; READ*, B B=A+B PRINT*,'B=A+B=',B WRITE(*,*) 'A*B=',A*B PRINT* ,'以上为计算机的计算结果，注意B的值'

Fortran95 简介-全文版

Fortran95简介－全文版 By陈鲸太FORTRAN的演进 FORTRAN的起源，要追溯到1954年IBM公司的一项计划。由JOHN BACKUS 领导的一个小组，尝试着在IBM 704计算机上面发展一套程序，它可以把使用接近数学语言的文字，翻译成机械语言。这个计划在刚开始并不被大家看好，但他们在1957年交出了成果，也就是第一套FORTRAN编译器，FORTRAN语言也就因此诞生了。FORTRAN语言的执行效率普遍的令各界满意，它证明了这项计划的可行性，也成为第一个被广泛使用的高级语言。FORTRAN的名字来自于英文的FORMULA TRANSLATOR这两个字，而这两个字恰是数学公式翻译器的意思。 旧版的FORTRAN77是在1978年由美国国家标准局(ANSI)所正式公布的，之后改版有1992年提出的FORTRAN90以及1997年的FORTRAN95，本文是为了FORTRAN 95所撰写。 编译器简介 1、VISUAL FORTRAN VISUAL FORTRAN一开始是起源于MICROSOFT的FORTRAN POWERSTATION 4.0，这套工具后来卖给DIGITAL公司来继续发展，下一个版本称为DIGITAL VISUAL FORTRAN 5.0，DIGITAL后来被COMPAQ合并，所以接下来的6.0及6.5版就称为COMPAQ VISUAL FORTRAN。而COMPAQ目前又跟HP合并，也许下一个版本会称为HP VISUAL FORTRAN。 VISUAL FORTRAN被整合在一个叫作MICROSOFT VISUAL STUDIO的图形接口开发环境中，VISUAL STUDIO提供一个统一的使用接口，这个接口包括文书编辑功能，PROJECT的管理、除错工具等等，所以在使用上其实跟上学期的VISUAL C++满类似的，同学们上课用过VISUAL C++，对 VISUAL FORTRAN应该不会陌生。 VISUAL FORTRAN 6.5除了完全支持FORTRAN 95的语法外，扩充功能方面提供完整的WINDOWS程序开发工具，专业版还内含IMSL数值链接库。 另外它还可以和VISUAL C++直接互相连结使用，也就是把FORTRAN和C 语言的程序代码混合编译成同一执行档案。 2、在工作站使用FORTRAN 学校计中工作站也提供FORTRAN COMPILER，使用方式很简单，只需要在存放FORTRAN档案的目录下面输入下面叙述即可：

中国海洋大学Fortran程序考试典型题

1、找出100-150之间和400-450之间能被9整除的数 2、找出100-999之间的回文数 3、找出水仙花数

4、输出“*”图形 5、输入n个数，找出大于平均值的数和最小数

6、把一个数值型数组的相同数删除到只剩一个 7、形成一个5*5矩阵，对角线元素为“i”，其余为“j”program juzhen implicit none integer I,J integer,parameter::size=5 integer::a(size,size) forall(I=1:size,J=1:size,I>J) a(I,J)=j forall(I=1:size,J=1:size,I==J) a(I,J)=i forall(I=1:size,J=1:size,I

integer,allocatable::c(:) write(*,*) "输入数组A的数据个数" read(*,*) n write(*,*) "输入数组B的数据个数" read(*,*) m l=m+n allocate (a(n)) allocate (b(m)) allocate (c(l)) write(*,*) "从小到大输入A的元素" do i=1,n read(*,*) a(i) end do write(*,*) "从小到大输入B的元素" do i=1,m read(*,*) b(i) end do do i=1,n c(i)=a(i) end do i=1 do while(i=c(n))then n=n+1 c(n)=b(i) goto 10 else if(b(i)>c(j).and.b(i)

FORTRAN程序设计复习题及答案

FORTRAN程序设计复习题 一、选择题 B (1)下列各FORTRAN表达式中合法的是 A) S+T*2P >= B) .NOT. (A*B+C) C) A2+B2/(C+D) <= D) (A+B).NOT．A*B．GT．(.NOT.只跟一个表达式) C (2)数学式(3/5)ex+y的FORTRAN表达式是 A) 3*EXP(X+Y)/5 B) 3*E* *(X+Y)/ C) (3/5)*EXP(X+Y)D) EXP(X+Y) D (3)下列FORTRAN77表达式中不合法的是 A) A．GT．B．EQV．C．GT．D B) A．AND．B．AND．C．AND．D C) .NOT．(X．LE．D) A．LT．B．LT．C．LT．D D(4)下列叙述中不正确的是 A) FORTRAN子程序可以单独编译 B) 对一个FORTRAN源程序进行编译和连接无误后可生成可执行文件 C) 即使编译和连接都正确无误，FORTRAN程序运行时仍可能出错 D) FORTRAN连接的主要任务是把函数库中的函数翻译成机器指令(正确描述：主要任务为连接目标文件) B (5)在下列FORTRAN77运算符中，优先级最高的是 A) .AND. B) .NOT. C) .OR. D) .EQ. B (6)FORTRAN表达式"6/5+9/2**3/2"的值为 A) 33 B) 1 C) 5 D) 3 A (7)下列FORTRAN77表达式中，合法的是： A) .AND.. B) 10.0 C) D) 提示：A）相当于 .AND.（.NOT.（）） D (8)关于编译一个FORTRAN源程序文件，下列说法中错误的是 A) 允许编译只有一个主程序而没有子程序的源文件 B) 允许编译有多个子程序的源文件 C) 允许编译只有一个子程序而没有主程序的源文件 D) 允许编译有多个主程序的源文件 C (9)在FORTRAN77源程序中，续行标志符必须放在 A) 第1列 B) 第1-6列C) 第6列D) 第5列 D (10)下列关于"SUBROUTIN E MAP(X，Y)"语句行的叙述中，不正确的是 A) 这是子程序的第一个语句 B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参D) 子程序执行后，MAP将返回整型数据 提示：子程序无返回值，自定义函数才有） A (11)FORTRAN表达式"2/4+"的值是 A) B) 1 C) D) 0 提示：2/4默认等于整型，=》 D (12)FORTRAN表达式"MOD，"的值是 A) B)0.0 C) D) A (13下列FORTRAN运算符中，优先级最低的是 A）逻辑运算符.AND. B）算术运算符*

Fortran编译器及下载地址

Fortran 编译器 --------------------------------------------------------------------------------------- 1. Fortran Powerstation 4.0 Fortran Powerstation 4.0是最老的版本 下载ftp://2006:2006@https://www.wendangku.net/doc/e010156287.html,/36/https://www.wendangku.net/doc/e010156287.html,-002123.ZIP --------------------------------------------------------------------------------------- 2.CVF Compaq Visual Fortran （CVF）, 当今PC平台上功能相当强大与完整的Fortran程序开发工具，还用于Abaqus的开发。 1997年，微软将Fortran PowerStation卖给DEC之后，微软就不再出版Fortran编译器了。后来DEC并入了Compaq，再后来Compaq又和HP合并了。现在最新的版本是HP出的Fortran for Windows v6.6，现在HP/Compaq已经不再开发Fortran了，CVF 6.6是最终的版本了，Compaq的Fortran开发小组已经投入Intel旗下，目前Intel 已经有Intel Visual Fortran 11.0。Compaq Visual Fortran 6.6官方的单价也相当昂贵。 Compaq Visual Fortran 6.6 下载： https://www.wendangku.net/doc/e010156287.html,/SoftDown.asp?ID=11937 Compaq Visual Fortran 6.6 绿色版下载： https://www.wendangku.net/doc/e010156287.html,/down/10915.html Compaq Visual Fortran 6.5 下载： https://www.wendangku.net/doc/e010156287.html,/soft/fortran6.5.rar ftp://2006:2006@https://www.wendangku.net/doc/e010156287.html,/36/https://www.wendangku.net/doc/e010156287.html,-002124.rar ---------------------------------------------------------------------------------------

几个fortran程序

计算圆周率 REAL R,R1,R2,PI ISEED=RTC() N0=0 N=300000 DO I=1,N R1=RAN(ISEED) R2=RAN(ISEED) R=SQRT(R1*R1+R2*R2) IF(R<1.0)N0=N0+1 END DO PI=4.0*N0/N WRITE(*,*)PI END 一）蒙特卡洛计算生日问题 假设有N个人在一起，各自的生日为365天之一，根据概率理论，与很多人的直觉相反，只需23个人便有大于50％的几率人群中至少有2个人生日相同。 INTEGER M(1:10000), NUMBER1(0:364), NUMBER2 REAL X,Y ISEED=RTC() DO J=1, 10000 NUMBER1=0 X=RAN(ISEED) NUMBER1(0)=INT(365*X+1) JJJ=1 DO I=1,364 Y=RAN(ISEED) NUMBER2=INT(365*Y+1) ETR=COUNT(NUMBER1.EQ.NUMBER2) IF (ETR= =1) THEN EXIT ELSE JJJ=JJJ+1 M(J)=JJJ NUMBER1(I)=NUMBER2 END IF END DO END DO DO I=1,10000 IF(M(I).LE.23) SUM=SUM+1 END DO PRINT *,SUM/10000 END 二）MONTE CARLO SIMULATION OF ONE DIMENSIONAL DIFFUSION 蒙特卡罗计算一维扩散问题

INTEGER X,XX(1:1000,1:1000) REAL XXM(1:1000) ! X:INSTANTANEOUS POSITION OF ATOM ! XX(J,I)：X*X ,J:第几天实验,I:第几步跳跃 ! XXM(I): THE MEAN OF XX WRITE(*,*) "实验天数JMAX，实验次数IMAX" READ(*,*) JMAX,IMAX ISEED=RTC() DO J=1,JMAX !第几天实验 X=0 !!! DO I=1,IMAX !第几步跳跃 RN=RAN(ISEED) IF(RN<0.5)THEN X=X+1 ELSE X=X-1 END IF XX(J,I)=X*X END DO END DO OPEN(1,FILE="C:\DIF1.DAT") DO I=1,IMAX XXM=0.0 XXM(I)=1.0*SUM(XX(1:JMAX,I))/JMAX !! WRITE(1,*) I, XXM(I) END DO CLOSE(1) END 三维的！ 三）通过该程序了解FORTRAN语言如何画图（通过像素画图）USE MSFLIB INTEGER XR，YR !在的区域中画一个圆 PARAMETER XR=400,YR=400 INTEGER R, S(1:XR,1:YR) X0=XR/2 ! 圆心位置X0，YO Y0=YR/2 R=MIN(X0-10,Y0-10) !圆半径 S=0 !像素的初始状态（颜色） DO I=1,XR DO J=1,YR IF((I-X0)**2+(J-Y0)**2<=R**2)S(I,J)=10 IER=SETCOLOR(S(I,J)) IER=SETPIXEL(I,J) END DO END DO END

Fortran95简介-全文版

Fortran95簡介－全文版 By陳鯨太FORTRAN的演進 FORTRAN的起源，要追溯到1954年IBM公司的一項計畫。由JOHN BACKUS 領導的一個小組，嘗試著在IBM 704電腦上面發展一套程式，它可以把使用接近數學語言的文字，翻譯成機械語言。這個計畫在剛開始並不被大家看好，但他們在1957年交出了成果，也就是第一套FORTRAN編譯器，FORTRAN語言也就因此誕生了。FORTRAN語言的執行效率普遍的令各界滿意，它證明了這項計畫的可行性，也成為第一個被廣泛使用的高階語言。FORTRAN的名字來自於英文的FORMULA TRANSLATOR這兩個字，而這兩個字恰是數學公式翻譯器的意思。 舊版的FORTRAN77是在1978年由美國國家標準局(ANSI)所正式公布的，之後改版有1992年提出的FORTRAN90以及1997年的FORTRAN95，本文是為了FORTRAN 95所撰寫。 編譯器簡介 1、VISUAL FORTRAN VISUAL FORTRAN一開始是起源於MICROSOFT的FORTRAN POWERSTATION 4.0，這套工具後來賣給DIGITAL公司來繼續發展，下一個版本稱為DIGITAL VISUAL FORTRAN 5.0，DIGITAL後來被COMPAQ合併，所以接下來的6.0及6.5版就稱為COMPAQ VISUAL FORTRAN。而COMPAQ目前又跟HP合併，也許下一個版本會稱為HP VISUAL FORTRAN。 VISUAL FORTRAN被整合在一個叫作MICROSOFT VISUAL STUDIO的圖形介面開發環境中，VISUAL STUDIO提供一個統一的使用介面，這個介面包括文書編輯功能，PROJECT的管理、除錯工具等等，所以在使用上其實跟上學期的VISUAL C++滿類似的，同學們上課用過VISUAL C++，對VISUAL FORTRAN應該不會陌生。 VISUAL FORTRAN 6.5除了完全支援FORTRAN 95的語法外，擴充功能方面提供完整的WINDOWS程式開發工具，專業版還內含IMSL數值程式庫。另外它還可以和VISUAL C++直接互相連結使用，也就是把FORTRAN和C語言的程式碼混合編譯成同一執行檔案。 2、在工作站使用FORTRAN 學校計中工作站也提供FORTRAN COMPILER，使用方式很簡單，只需要在存放FORTRAN檔案的目錄下面輸入下面敘述即可： ccsun33 [u8623033/fortran]% f77 filename.for

fortran程序30个

程序1 PROGRAM chaper1272 INTEGER :: A(10)=(/5,7,4,8,12,2,10,3,9,11/) INTEGER :: sum=0,ave !打开一数据文件,设置一个有格式直接存取文件,将10个数分2个记录写入文件。 !数据文件生成2个记录,每个记录行长度相同,记录长度为25 OPEN(1,FILE='input21.dat',FORM='FORMATTED',ACCESS='D IRECT',RECL=25) WRITE(1,"(5I5)",REC=1)(A(I)+10,I=1,5) !按格式说明将头5个数写入第1个记录 WRITE(1,"(5I5)",REC=2)(A(I)+10,I=6,10)!按格式说明将后5个数写入第2个记录 READ(1,"(5I5)",REC=2)(A(I),I=6,10) !按格式说明从第2个记录中读取后5个数 READ(1,"(5I5)",REC=1)(A(I),I=1,5) !按格式说明从第1个记录中读取头5个数 DO I=1,10 sum=sum+A(I) ENDDO ave=sum/10 !打开一个最大记录长度为22的有格式顺序存取文件

OPEN(2,FILE='input22.dat',FORM='FORMATTED',ACCESS='D IRECT',RECL=22) WRITE(2,"('10个数之和为：',I5)",REC=1) sum !输出1记录行,记录长度为22 WRITE(2,"('10个数平均值为：',I5)",REC=2)ave !输出1记录行,记录长度为22 WRITE(2,"(A)",REC=3) '程序运行正常结束。' !输出一个记录行,记录长度为22 END 程序2 PROGRAM average REAL sum,ave INTEGER n OPEN(1,file='score.dat') PRINT*,'正在统计平均成绩,请等待。' sum=0.0;n=0 DO READ(1,*,END=100) S sum=sum+s n=n+1 ENDDO 100 ave=sum/n

FORTRAN 95 语法速查

FORTRAN 95 语法速查 ----------白云、李学哲、陈国新、贾波编著《FORTRAN95程序设计》读书笔记 目录：一、应用程序的创建与运行/FORTRAN 95所用的字符/ 变量类型及其声明，常量声明/表达式与运算符 二、输入与输出：表控、有格式 三、选择语句与结构：IF语句、CASE结构 四、DO循环结构 五、数组：数组的声明，数组的引用，数组的算术运算，数组的输入、输出，给数组赋初值， 动态数组，WHERE、FORALL语句 六、子程序：语句函数，内部子程序，调用子程序时的虚实结合：形参为数组、非定界数组、子 程序名、星号，递归子程序，外部子程序，纯子程序，逐元子程序 七、派生数据类型与结构体 八、指针与动态链表 九、文件：存取方式，基本操作语句，各类文件的读写操作 十、接口、模块 十一、公用区、存储关联、数据块子程序 十二、绘图：坐标系、设置图形颜色、创建图形程序/ 常用过程：设置线型、绘一像素点、设置当前位置、绘直线、绘弧线、绘矩形、绘多边形、绘制扇形（圆、椭圆）/ 文字信息的显示 附/录：标准函数与标准子例行程序 一、基础部份 1-1 FORTRAN 95 应用程序的创建与运行 创建或运行FORTRAN 95程序必须在Microsoft Developer Studio平台上进行。尽管程序文本及相关文件的编辑可以在任一文本编辑器上进行，然后再拷到Studio的文档窗口中。但最好还是一开始就进入Studio环境。创建FORTRAN 95 程序的步骤大致如下： 1）启动Microsoft Developer Studio 可以通过不同方式运行dfdev.exe程序以启动Microsoft Developer Studio [开始] \ Compaq Visual Fortran 6 \ Developer Studio \ dfdev.exe：或 ……\CVF66 \https://www.wendangku.net/doc/e010156287.html,\MSDEV98\dfdev.exe Microsoft Developer Studio的界面如下图所示： 文档窗口 工作空间窗口 输出窗口

Fortran基本用法小结档

Fortran基本用法小结 目录： 一、说明 二、概述 三、数据类型及基本输入输出 四、流程控制 五、循环 六、数组 七、函数 八、文件 一、说明 本文多数内容是我读彭国伦《Fortran 95 程序设计》的笔记。只读到第九章，主要是3~9章，都是些最基本的用法（原书共16章）。这里主要摘录了我看书过程中总结的一些Fortran和C不同的地方，主要是语法方面。希望这份笔记能够给学过C但没有接触过Fortran的同学带去一些帮助。要想弄得更清楚些，推荐看一下原书，觉得作者真的写得很好，很清楚；如果有C语言的基础，看完前九章应该很快的，花一两天就行了。觉得如果耐心看完本文，基本功能应该也可以顺利用起来了。 另外，由于我之前没有用过Fortran，这次为了赶文档看书又看得很粗浅，大多数东西看过之后都没来得及仔细想，只是按着作者的意思去理解。所以这份笔记还处于纸上谈兵的层次。如果有不妥的地方，希望大家指正。谢谢！ 文中蓝色的部分是程序代码，!后面的内容为注释。 二、概述 1、名词解释 Fortran=Formula Translator/Translation 一看就知道有什么特色了：可以把接近数学语言的文本翻译成机械语言。的确，从一开始，IBM设计它的时候就是为了方便数值计算和科学数据处理。设计强大的数组操作就是为了实现这一目标。Fortran奠定了高级语言发展的基础。现在Fortran在科研和机械方面应用很广。 2、Fortran的主要版本及差别 按其发展历史，Fortran编译器的版本其实很多。现在在广泛使用的是Fortran 77和Fortran90。Fortran 90在Fortran 77基础上添加了不少使用的功能，并且改良了77编程的版面格式，所以编程时推荐使用90。鉴于很多现成的程序只有77版本，有必要知道77的一些基本常识，至少保证能够看懂77程序。以下是77和90的一些格式上的区别。

fortran课后习题答案

第一章 FORTRAN程序设计基础第15页 1、2 1．简述程序设计的步骤。 “程序设计”：反映了利用计算机解决问题的全过程，通常要经过以下四个基本步骤：(1)分析问题，确定数学模型或方法；(2)设计算法，画出流程图；(3)选择编程工具，编写程序；(4)调试程序，分析输出结果。 2. 什么是算法？它有何特征？如何描述算法？ 解决问题的方法和步骤称为算法。 算法的五个特征：(1) 有穷性。 (2) 确定性。 (3) 有效性。 (4) 要有数据输入。(5) 要有结果输出。 算法的描述有许多方法，常用的有：自然语言、一般流程图、N－S图等。 第二章顺序结构程序设计 第29页 1、2、3、4、5、6、7、8、9 1．简述符号常量与变量的区别？ 符号常量在程序运行过程中其值不能改变。变量在程序运行过程中其值可以改变。 2. 下列符号中为合法的FORTRAN 90标识符的有哪些? (1) A123B (2) M%10 (3) X_C2 (4) 5YZ (5) X+Y (6) F(X) (7) COS(X) (8) A.2 (9) ‘A’ONE (10) U.S.S.R. (11) min*2 (12) PRINT 3. 下列数据中哪一些是合法的FORTRAN常量？ (1) 9,87 (2) .0 (3) 25.82(4) -356231 (5) 3.57*E2 (6) 3.57E2.1 (7) 3.57E+2(8) 3,57E-2 4. 已知A=2,B=3,C=5(REAL)；且I=2,J=3(INTEGER),求下列表达式的值: (1) A*B+C 表达式的值: 11 (2) A*(B+C) 表达式的值: 16 (3) B/C*A 表达式的值: 1.2 (4) B/(C*A) 表达式的值: 0.3 (5) A/I/J 表达式的值: 0.33 (6) I/J/A 表达式的值: 0 (7) A*B**I/A**J*2 表达式的值: 4.5 (8) C+(B/A)**3/B*2. 表达式的值: 7.25 (9) A**B**I 表达式的值: 512 5. 将下列数学表达式写成相应的FORTRAN表达式： (1) 1E-2 （2）(-B+SQRT(B*B-4*A*C)/(2*A) (3) 1+X+X*X/2+X**3/2/3 (4) COS(ATAN((A**3+B**3)**(1.0/3)/(C*C+1))) (5) EXP(A*X**2+B*X+C) (6) COS(X*Y/SQRT(X*X+Y*Y))**3 6. 用FORTRAN语句完成下列操作： (1) 将变量I的值增加1。I=I+1 (2) I的立方加上J,并将结果保存到I中。 I=I**3+J (3) 将E和F中大者存储到G中。G=Max(E,F) (4) 将两位自然数N的个位与十位互换，得到一个新的数存储到M中(不考虑个位为0的情况) M=MOD(N,10)*10+N/10 第三章选择结构程序设计第43页 1、2、3、5、6、7、9 1．分析下列程序运行结果 (1) LOGICAL P INTEGER I,I1,I2,I3 P=.FALSE. READ*,I I1=MOD(I,10) I2=MOD(I/10,10) I3=I/100

fortran习题参考答案

工程分析程序设计上机作业（五） 模块化编程 上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。 1、有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面 积 2 a b c s ++ =，a、b、c为三个边长）。要求用内部函 数来计算每个三角形的面积。 =21 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 function end 2、编写一个子例程子程序，SUM(S，T，N1，N2)。把整型数N1到N2进行求和，并把求 和的结果放置到S，把整型数N1到N2进行求积，并把求积的结果放置到T。并用这个子程序来计算： y=++++++++++- (1234)(345678)(3*4*5*6)(1*2*3) 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)=i end 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，这时的除数就是两个数的最大公约数。