当前位置：文档库 › fortran习题集参考答案解析

fortran习题集参考答案解析

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

模块化编程

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

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

积

a b c

=，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，这时的除数就是两个数的最大公约数。

调用此函数，求1260，198，72三个数的最大公约数。

program main

implicit none

integer a,b,c,m,n,temp

print*,'请输入三个数：'

print*,'a b c'

read*,a,b,c

if (a

temp=a

a=b

b=temp

end if

call GDC(a,b,m)

if (c

temp=c

c=m

m=temp

end if

call GDC(m,c,n)

print*,'最大公约数为',n

contains

subroutine GDC(x,y,z) integer x,y,z,k

do while(mod(x,y)/=0) k=mod(x,y)

x=y

y=k

end do

z=y

end subroutine end

4、编写一个模块程序，提供以下服务：定义出常量π、e 。定义出子程序，实现求和21n

i n =∑、

求阶乘n !。并在主程序中计算如下结果：从键盘上输入整数n 、实型数A 、R 、R 0，求

21!

n i n n =?? ???∑（实型）和0202n n R R An R e R R π??- ????? ???（实型）。

module mytask

implicit none

real,parameter :: PI=3.1415926

real,parameter :: e =2.7182818

contains

function sub(n)

integer i,n,sub

sub=0

do i=1,n

sub=sub+i*i end do end function sub

function fact(n)

integer i,n

real fact

fact=1

do i=1,n

fact=fact*i

end do

end function fact

end module mytask

program main

use mytask

implicit none

integer n

real A,R,R0

real x,y

print*,'请输入依次n,A,R,R0 :'

read*,n,A,R,R0

x=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)=',x

print*,'(A*n/(2*PI*R*R))*(R/R0)**n*e**(-(R/R0)**n)',y end program

5、编写外部函数，计算...!

7!5!3sin 7

53+-+-=x x x x x 的值，直到最后一项的绝对值小于10-6为止。再编写主程序，从键盘读入x ，调用该外部函数并输出sinx 的计算结果。注意不能用Fortran 的标准函数SIN(X)。

program main

implicit none

real x,p

external sub

print*,'请输入角度 x :'

read*,x

p=x/180*3.14

call sub(p)

end program

subroutine sub(k)

implicit none

real k,y

integer i

y=0

i=1

do while (abs(k**i/(fact(i)))>1.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,fact

fact=1

do j=1,m

fact=fact*j

end do

end function fact

end subroutine sub

6、编写程序，用欧拉法数值求解一阶微分方程。

数值解法的基本思想：用差分方程代替微分方程，然后在若干个离散点上逐点求解差分方程，得到各离散点x 0, x 1, x 2, … 等处函数(),()dy f x y x dx

=的近似值y 0, y 1, y 2, … 。其中各离散点x 0, x 1, x 2, …之间的距离称为步长h 。

欧拉法的基本思想是，一阶微分

dy dx 可用向前差分代替，即11'()n n y y y h

+≈-。带入微分方程，可得 1(,)n n n n y y f x y h

+-≈ 因此，1*(,)n n n n y y h f x y +≈+，其中n n x x h -=+1。

当给定x 0, y 0和步长h 后，即可按下面步骤求得数值解：

0000(,)x f x y y ?????1111,(,)y f x y x ????()2222,n y f x y y x ???????

请按上述方法编写程序用Euler 法求解微分方程22dy y x dx

=-，当x 0=0时，y 0=1.0。取h =0.1，试求出当x i =0.1,0.2,0.3,0.4….,1.0时的y i 值。

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))

print*,'x(',j,')=',x(j),'y(',j,')=',y(j) end do

contains

function f(a,b)

real a,b,f

f=b*b-a*a

end function

end

fortran调试经验

FORTRAN调试程序的时候注意的问题调试程序的时候注意的问题。程序编好，能够直接运行而且完全正确，基本不可能，这就有调试程序的问题。主要有一下几个方面：其一，在每个子程序被调用的时候留个心眼，写个输出语句，表示程序已经运行到了这里。这样一个小提示会给调试带来巨大的方便，如果程序运行出错，至少你可以知道它是在运行到哪里出的错，这样，直接去检查那个程序就可以了。不必重头开始检查。其二，注意对中间计算结果的输出。有时候，而且是很多的时候，程序编译成功，运行也没有问题，就是结果不对，这肯定是计算原理有问题，此时，输入一些重要步骤的中间结果，往往可以检查出问题所在。同时，就算查出了问题所在，也可以不删除这段输出中间计算结果的代码，有可能后面还会有用处，此时，在每行输出语句前加一个感叹号，把代码变成注释的绿体字就可以了。等到再次需要输出的时候，直接删除一个“！”比再写一遍输出代码，当然要简单的多。其三，对WATCH功能的应用，FORTRAN提供的这个功能很实用，可以查很多问题，尤其是程序中间计算值，这个和上述的中间结果的输出有点相似。但两者的不同是前者可以进行中间结果的输出控制，就是只有符合了某个条件的才能被输出，这样可以便捷程序的调试，同时对中间结果输出后可以直接用STOP停止程序的运行，这样对于大型程序来说，节省了很多后面继续计算的时间——因为前面的结果已经不对了，后面的算也是白算。其四，对中间结算结果输出形式的控制，一般来说，FORTRAN计算结果可以输出到文件里面和计算界面两个地方。对较大的计算结果，可以输出到文件里面，反之较少的结果可以直接输出到屏幕上，为了增强数据的可读性，最好进行有格式的数据输出，以利于相同性质的数据的比较。输出到屏幕上的结果直接用WRITE(6,*)就可以（无格式），对于输出到文件里面的数据，可以省些事情，直接用WRITE(X,*)就可以，其中X是一个任意的正整数，最好大于10，也不用事先对这个X设备进行说明，程序会将结果输出到一个FORT.X的文件里面，例如10，就是FORT.10，此时，用NOTEPAD或者ULTRA-EDIT都可以把它打开——FORT.10实质上就是一个.DAT的文件，你可以把它重命名。 3.对数据计算时的误差控制。以前觉得小数点后的误差不是那回事，没有太在意，可经过实战，终于明白了小小的误差完全可以改变整个计算的结果。因此，如果程序能够输入结果而不正确时，除了寻找算法的问题，不要忽略的误差。一般认为，FORTRAN的REAL变量小数点后8位数字误差定义已经足够，而事实上,这个精度可能在一些情况下不满足，这个时候，需要用更精确的变量类型——REAL(8),同理，当要判断两个数是否相等的时候，一定要慎用相等判断（.EQ.）这个比较运算符，因为任何数据，别看着在现实中它们一定相等，在程序中就不一定了。一旦经过了计算，就不可避免的产生了舍入误差，对于整数和有限几位循环的有理数都问题不大，可一旦是一个无理数或者无限循环的小数，只有在判断了小数点后的每一位都相等的时候，程序才判断为相等成立。这个相等的标准是非常苛刻的，所以一般情况下，可行的方法是将

Fortran95习题答案

(完整)Fortran经典编程语言笔记(你值得拥有)

FORTRAN笔记 2014.10.07 目录第七讲_FORTRAN的基本知识.ppt (2) FORTRAN语言程序设计初步 (2) FORTRAN源程序的书写格式（以77为例） (2) 变量 (2) 变量类型 (2) 算术运算符和运算优先级 (3) 赋值语句 (3) 参数语句（PARAMETER语句） (3) END语句 (3) PAUSE语句 (3) 逻辑运算和选择结构 (4) 关系表达式 (4) FORTRAN中数组的定义及使用 (4) 其他 (5) 1. fortran语言定义CHARACTER*6 TTL(14,3),CNAM(400)是什么意思? (5) 2. fortran里character*10 是什么意思 (5) 3. Fortran中kind是什么函数？ (5)

第七讲_FORTRAN的基本知识.ppt FORTRAN语言程序设计初步 FORTRAN是Formula Translation的缩写，意为“公式翻译”，它是为科学、工程问题或企事业管理中的那些能够用数学公式表达的问题而设计的，其数值计算的功能较强。常用的是FORTRAN77和FORTRAN90两种标准。 1、一个程序由若干个程序单位组成。主程序和每一个子程序分别是一个独立的程序单位。 2、每一个程序单位都是以“END”结束的。 3、一个程序单位包括若干行。 1）语句行。由一个FORTRAN语句组成。 2）非语句行，即注释行。 4、FORTRAN程序中的语句可以没有标号，也可以有标号，根据需要而定。标号的作用是标志一个语句以便被其他语句引用。 5、一个程序单位中各类语句的位置是有一定规定的。 6、FORTRAN源程序必须按一定的格式书写。 FORTRAN源程序的书写格式（以77为例）每一行有80列，分别如下： 1、第1-5列为标号区。一行中第一列为“C”或“*”，该行即被认为是注释行。 2、第6列为“续行标志区”，如果在一行的第6列上写一个非空格和非零的字符，则该行作为其上一行的续行。 3、第7-72列为语句区。 4、第73-80列，注释区。变量变量名：一个变量需要用一个名字（变量名）来识别。在同一个程序单位中不能用同一个变量名代表两个不同的变量。 FORTRAN的变量名按以下规则选定： 1）第一个字符必须是字母，即变量名必须以字母开头； 2）在一个字母后面可以跟1-5为数字或字母。如果选定的变量名超过6个字符，则只有前面6个字符有效。注：在变量名中大写与小写字母是等价的。变量类型整型变量Integer、实型变量Real、双精度变量Double Precision、复型变量Complex、逻辑型变量Logical和字符型变量Character。 1、隐含约定（I-N规则） FORTRAN规定：在程序中的变量名，凡以字母I，J，K，L，M，N六个字母开头的，即认为该变量为整型变量。在程序中，凡是变量名以字母I,J,K,L,M,N,i,j,k,l,m,n开头的变量被默认为整型变量，以其他字母开头的变量被默认为实型变量。 2、用类型说明语句确定变量类型 1）INTEGER语句（整型说明语句） 2）REAL语句（实型说明语句） 3）DOUBLE PRECISION语句（双精度说明语句） 4）COMPLEX语句（复型说明语句） 5）LOGICAL语句（逻辑型说明语句）

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语言语法

FORTRAN是世界上最早出现的高级编程语言，是工程界最常用的编程语言，它在科学计算中（如航空航天、地质勘探、天气预报和建筑工程等领域）发挥着极其重要的作用。经过40多年的发展，伴随着FORTRAN语言多次版本的更新及相应开发系统的出现，其功能不断完善，最新版本的开发系统几乎具备了VC、VB的所有特点，如图形界面编程、数据库等。目前，工科院校开设的计算机编程语言课首选仍然是FORTRAN :< 说实话，从科技发展的趋势来说这不是好事。您可以设想一下，如果需要用鹅毛笔抄写大量的古籍是什么感受！强烈建议阅读《发掘C#特性赋予科学计算项目以威力》 1 FORTRAN77四则运算符 + - * / ** (其中**表示乘方) 在表达式中按优先级次序由低到高为： +或-→*或/→**→函数→() 2 FORTRAN77变量类型隐含约定：I-N规则凡是以字母I，J，K，L，M，N六个字母开头的，即认为是整型变量，其它为实型变量。用类型说明语句确定变量类型：可以改变I-N规则

用IMPLICIT语句将某一字母开头的全部变量指定为所需类型如IMPLICIT REAL (I,J) 三种定义的优先级别由低到高顺序为：I-N规则→IMPLICIT语句→类型说明语句，因此，在程序中IMPLICIT语句应放在类型说明语句之前。数组的说明与使用使用I-N规则时用DIMENSION说明数组，也可在定义变量类型同时说明数组，说明格式为：数组名(下标下界,下标上界），也可省略下标下界，此时默认为1，例：DIMENSION IA(0:9),ND(80:99),W(3,2),NUM(-1:0),A(0:2,0:1,0:3)

Fortran95程序设计课后习题答案(word版方便)

第四章 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.p rogram 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 implicit none integer money real tax write(*,*) "请输入月收入" read(*,*) money if ( money<1000 ) then tax = 0.03 else if ( money<5000) then tax = 0.1 else tax = 0.15 end if write(*,"(' 税金为'I8)") nint(money*tax) end program 2.program main implicit none integer day character(len=20) :: tv write(*,*) "请输入星期几" read(*,*) day select case(day) case(1,4) tv = "新闻" case(2,5) tv = "电视剧" case(3,6) tv = "卡通" case(7) tv = "电影" case default write(*,*) "错误的输入" stop end select write(*,*) tv end program 3.program main implicit none integer age, money real tax write(*,*) "请输入年龄" read(*,*) age write(*,*) "请输入月收入" read(*,*) money if ( age<50 ) then if ( money<1000 ) then tax = 0.03 else if ( money<5000 )then tax = 0.10 else tax = 0.15 end if else if ( money<1000 ) then tax = 0.5 else if ( money<5000 )then tax = 0.7 else tax = 0.10 end if end if write(*,"(' 税金为'I8)") nint(money*tax) end program 4.program main implicit none integer year, days logical mod_4, mod_100, mod_400 write(*,*) "请输入年份" read(*,*) year mod_4 = ( MOD(year,4) == 0 ) mod_100 = ( MOD(year,100) == 0 ) mod_400 = ( MOD(year,400) == 0 ) if ( (mod_4 .NEQV. mod_100) .or. mod_400 ) then days = 366 else days = 365 end if write(*,"('这一年有'I3'天')") days stop end program 第六章 1.program main implicit none integer i do i=1,5 write(*,*) "Fortran" end do stop end program 2.program main implicit none integer i,sum sum = 0 do i=1,99,2 sum = sum+i end do write(*,*) sum stop end program 3.program main implicit none integer, parameter :: answer = 45 integer, parameter :: max = 5 integer weight, i do i=1,max write(*,*) "请输入体重" read(*,*) weight if ( weight==answer ) exit end do if ( i<=max ) then write(*,*) "猜对了" else write(*,*) "猜错了" end if stop end program 4.program main implicit none integer, parameter :: max=10 integer i real item real ans ans = 1.0 item = 1.0 do i=2,max item = item/real(i) ans = ans+item

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

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常见错误

FAQ之常见错误 2014-02-02 13:45:35 来源：Fcode研讨团队评论：2点击：4419 本文从编译错误，链接错误，运行时错误，计算结果错误等四个方面介绍了常见的错误及解决思路。适合初学者阅读。首先应该明确：错误有哪几种？我们当前遇到的是何种错误？阐述这些问题前，我们先讨论一下常规的应用程序开发的过程： 1>>编写代码，使用一个或多个源代码文件。 2>>对第一步的每一个源代码文件执行编译操作。得到一个或若干个目标代码。 3>>将目标代码，运行时库（Run-time Library）和其他使用到的函数库链接起来。得到一个可执行文件（EXE 或其他） 4>>编写程序的说明书，必要的（输入）数据文件 5>>将上述得到的结果发布给用户。（发布的方式可以是刻录成光盘，销售，放在网站上供别人下载，或者其他） 6>>用户得到程序后，运行，输入数据，得到计算结果。对于很多 Fortran 程序员来说，可能用户就是自己，也可能仅仅是自己教研室的同事同学。所以第4，5，6步骤很多时候不明显。而如果使用集成开发环境（IDE）进行开发，第1，2，3步骤又可以一键完成。因此，很多初学者就认为，写程序就是：输入代码，运行，得到结果。这样的理解太狭义。不管我们面对什么使用者来写代码，程序开发应该是上述的过程。我们的编译器，编译环境，也是为这个过程而设计的。于是，我们将错误分为四种：一. 编译错误（发生在第2步）编译错误，一般是源代码书写格式不正确，不符合语法要求。二. 链接错误（发生在第3步）链接错误，一般是源代码结构不完整，运行时库或函数库使用不合理。三. 运行时错误（发生在第6步）运行时错误，一般是执行代码时，遇到了事先未料及的错误。比如内存不足了，磁盘空间不够了，输入文件格式不对了，输出文件写入失败了等等。四. 计算结果不符合预期（程序代码不规范，或不符合你的设想）计算结果不符合预期，可能性就很多了。语法与你的想法不一致，超出函数库的适用范围，执行流程控制不当等等。这四种错误，其排查难度依次增大。也就是，编译错误最容易排查和修改，而计算结果不正确，最让人头疼。

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中批处理命令的实现函数：利用systemqq命令（需要调用DFLIB 数据库） ********************************************* 例1： USE DFLIB character*100 CMD LOGICAL(4) res CMD="dir/a-d/b/s "//trim(fPath)//" >"//trim(outPut) res=SYSTEMQQ(CMD) 例2： USE DFLIB LOGICAL(4) result result = SYSTEMQQ('copy e:\dir.txt e:\test\dir.txt') !将e:\dir.txt 复制到e:\test\dir.txt文件中。！****************实例3:复制文件************************* 例3： programmain_pro USE DFLIB implicit none integer,parameter::sta_num=123 character(5),dimension(sta_num)::sta_ID character(500)::filein,fileout character(5000)::cmd logical(4)::judge

integer::status,is open(1,file='山东.txt',status='old',action='read',iostat=status) read(1,*) do is=1,sta_num read(1,*) sta_ID(is) filein='Z:\data\降水逐小时数据-戴至修\precip_data\'//sta_ID(is)//'_precip.txt' open(2,file= filein,status='old',action='read',iostat=status) if(status/=0) goto 1000 fileout='Z:\data\降水逐小时数据-戴至修\山东省-降水数据\'//sta_ID(is)//'_precip.txt' cmd='copy '//filein//' '//fileout judge=SYSTEMQQ( cmd) 1000 continue enddo end program

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语言内部函数

附录 FORTRAN 90标准函数符号约定： ●I代表整型;R代表实型;C代表复型;CH代表字符型;S代表字符串;L代表逻辑型;A代表数组;P代表指针;T代表派生类型;AT为任意类型。 ●s:P表示s类型为P类型(任意kind值)。s:P(k)表示s类型为P类型(kind值=k)。 ●[…]表示可选参数。 ●*表示常用函数。表1 数值和类型转换函数函数名说明 ABS(x)*求x的绝对值∣x∣。x:I、R,结果类型同x; x:C,结果:R AIMAG(x)求x的实部。x:C,结果:R AINT(x[,kind])*对x取整,并转换为实数(kind)。x:R, kind:I,结果:R(kind) AMAX0(x1,x2,x3,…)*求x1,x2,x3,…中最大值。x I:I,结果:R AMIN0(x1,x2,x3,…)*求x1,x2,x3,…中最小值。x I:I,结果:R ANINT(x[,kind])*对x四舍五入取整,并转换为实数(kind)。x:R, kind:I,结果:R(kind) CEILING(x)*求大于等于x的最小整数。x:R,结果:I CMPLX(x[,y][,kind]))将参数转换为x、(x,0.0)或(x,y)。x:I、R、C, y:I、R,kind:I,结果:C(kind) CONJG(x)求x的共轭复数。x:C,结果:C DBLE(x)*将x转换为双精度实数。x:I、R、C,结果:R(8) DCMPLX(x[,y])将参数转换为x、(x,0.0)或(x,y)。x:I、R、C, y:I、R,结果:C(8) DFLOAT(x)将x转换为双精度实数。x:I,结果:R(8) DIM(x,y)*求x-y和0中最大值,即MAX(x-y,0)。x:I、R, y的类型同x,结果类型同x DPROD(x,y)求x和y的乘积,并转换为双精度实数。x:R, y:R,结果:R(8)

计算机程序设计语言(FORTRAN语言)

计算机程序设计语言(FORTRAN语言) 1. 编译程序能将高级语言编写的源程序转换成。 A.解释程序 B.汇编程序 C.映象程序 D.目标程序答案：D 2. 一个完整的FORTRAN源程序。 A.至少包括一个主程序 B.至少包括一个主程序和一个子程序 C.由一个主程序与一个以上的子程序组成 D.由一个主程序与一个子程序组成答案：A 3. 语句函数定义语句在程序内合法的位置是。 A.在程序块开头语句之后，END语句之前 B.在程序块开头语句之后，可执行语句之前 C.在说明语句之后，END语句之前 D.在说明语句之后，可执行语句之前答案：D 4. 下列关于函数子程序虚实参数的错误说法是。 A.可以没有形参数 B.虚实结合的数组长度可以不同

C.实参表与虚参表类型可以不同 D.函数名可以作为虚参答案：C 5. 下列叙述中正确的是。 A.FORTRAN程序块中，无名公用语句只能有一个 B.FORTRAN子程序中，至少应有一个RETURN语句 C.FORTRAN程序块中，最后一行必须是END语句 D.FORTRAN程序块中，必须有变量说明语句答案：C 6. 运行下面的程序时得不到所需的结果，其主要原因是。 INTEGER X(11) DATA X/9，8，7，6，5，4，3，2，1，0，-1/ DO 10 1=1，X(1) ，-1 WRITE(*，*)1.0/SQRT(25.0-REAL(X(1) )* * 2) 10 CONTINUE END A.没有给X(11) 赋初值 B.发生除以零的情况 C.发生负数开平方的情况 D.循环参数设置错误答案：D 7. 下列数据中，不符合FORTRAN常量表示法的是。 A.-25.6 B.2.758D3

fortran问题

imsl7.0中用use linear_operators这句话是会出错的，当时intel论坛上也有人问，后来intel 给出了个X64 imsl的补丁，但是32位的没有。要使用包含在linear_operators这个库中的函数时，要用use+原函数例如：上面那个例子把use linear_operators改成use operation_xt即可另外，imsl7.0引用函数和以前版本不一样的补充一点： linear_operators这个文件是这样子的（一看就知道原因了）： modulelinear_operators usecond_int usedet_int usediag_int usediagonals_int usefft_int useifft_int useeye_int uselin_eig_self_int uselin_sol_self_int usenorm_int useoperation_i useoperation_ix useoperation_t useoperation_h useoperation_tx useoperation_hx useoperation_x useoperation_xi useoperation_xt useoperation_xh useorth_int userand_int userank_int usesvd_int useunit_int useeig_int usechol_int useisnan_int end module 1. 如何加大Stack size? 选Project => Settings => Link => Category: Output =>

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/fc2728019.html,\MSDEV98\dfdev.exe Microsoft Developer Studio的界面如下图所示：文档窗口工作空间窗口输出窗口