fortran中epsilon函数

一个语言可以从很多方面来介绍，但由于你没有提你想了解哪些方面，所以下面从①性质和②特点及应用场景 两方面来简单介绍。①从性质上来说，fortran是一种编译型语言。然后它是一种面向过程的编译型语言。编译型语言即将代码直接编译和汇编成机器码的语言，而非送到解释器逐行执行的解释型语言。编译型语言总共也没几种，比较知名的只有C、PACAL、Fortran。而编译型语言和解释型语言的区别及优劣我在这里就不展开讲了，你可以自己再百度一下。②从应用场景来说，虽然Fortran是最早的几种高级语言之一，但Fortran现在已经完全和商业化无关了，它目前只被应用在科学计算领域，或者更确切点说它现在只被应用在科研用的数值模式、数值计算领域。在科研领域，我们总是需要手动设计并编写一些新的而且是复杂的数值模式（而商业软件只是集成一些确定的可行的东西），在这种情况下，Fortran与C++效率相当但书写更方便，因此目前大部分复杂的数值模式还使用Fortran编写。

Fortran在工程软件比如Aspen以及大型计算软件比如DFT软件占据统治地位。

fortran是计算机领域最早的高级语言，其标准从fortran66、fortan77、fortran90、95，一直到最新的2000、2003。目前fortran2008标准已经在一些编译器上部分实现。从它的不断发展来看，fortran肯定没有被淘汰。相反，它的不断发展也在继续适应各方面的需求。从2003标注开始，fortran已经跟C++等一样，成为了面向对象的现代计算机编程语言。fortran在长期应用中积累了大量的算法资源，尤其在数值运算、数组处理方面。因此，它在航天、航空、气象、数值模拟、并行计算等方面具有很好的应用和前景。

如果你今后想从事的是软件工程师的工作，做专业的软件开发，那么毫无疑问应该学习C语言。但如果仅仅是为了你自己的本专业中相关场合的计算问题，可以考虑学习fortran语言，它的计算能力要强于其他编程语言（因为这是专门为科学计算而发明的语言），但通用性很有限，并不适合做其他工作。

http://baike.baidu.com/view/36402.html?wtp=ttFORTRAN 开放分类： 编程、语言、计算机、程序设计、科学计算 目录u2022 FORTRAN简介 u2022 FORTRAN开发历史 u2022 Fortran的特性 u2022 Fortran语言的Hello World程序 u2022 Fortran编译器 u2022 Fortran程序包

FORTRAN语言中文叫公式翻译语言，是IBM公司发明的计算机高级语言。是我国六七十年代流行的两大编程语言之一,用于较大的机种。七十年代末起，用于PC。FORTRAN语言逻辑性强，程序结构清晰，语法语义简捷好懂，特别适合用于科学计算，数据采集处理，调用绘图库（例如GKS,DISPLA等）可以绘图 。大型 MainFrame 计算机，DEC 计算机 等都用 FORTRAN。从FORTRAN90开始，加入了可视化。现在叫Compag Visual FORTRAN。 可以在 MS Visual Studio 里工作。凡C,C++ 能干的活，FORTRAN都能做。近十年来 FORTRAN 不如 C,C++ 流行。主要因为FORTRAN编译器价格贵。不像C,C++ 有免费的版本。unix 和 工作站流行推动了C的传播，同时MainFrame，DEC的退役，老互联网BITNET消亡，使FORTRAN失去了主要的硬件基础。微软视窗统治了PC世界，与微软视窗紧密结合的VC++比FORTRAN好用。FORTRAN能干的活，C,C++ 都能干。从目前势头看，尽管FORTRAN还能并存一段时期，但终将走向消亡。FORTRAN程序例子：write(*,*) "Hello World !"end上面是一个完整的FORTRAN程序，打印：Hello World !

FORTRAN，亦译为福传，是英文“FORmula TRANslator”的缩写，译为“公式翻译器”，它是世界上最早出现的计算机高级程序设计语言，广泛应用于科学和工程计算领域。FORTRAN语言以其特有的功能在数值、科学和工程计算领域发挥着重要作用。 Fortran的特性　　Fortran语言的最大特性是接近数学公式的自然描述，在计算机里具有很高的执行效率。 　　易学，语法严谨。 　　可以直接对矩阵和复数进行运算，这一点类似matlab。 　　自诞生以来广泛地应用于数值计算领域，积累了大量高效而可靠的源程序。 　　很多专用的大型数值运算计算机针对Fortran做了优化。 　　广泛地应用于并行计算和高性能计算领域。 　　Fortran90，Fortran95，Fortran2003的相继推出使Fortran语言具备了现代高级编程语言的一些特性。 Fortran的将来　　Fortran语言是一种极具发展潜力的语言，在全球范围内流行过程中，Fortran语言的标准化不断吸收现代化编程语言的新特性，并且在工程计算领域仍然占有重要地位。 　　不可否认，Fortran语言与目前流行的JAVA,C#等高级语言相比，它缺乏创造力。但是，工科学生必须注意到，由于Fortran在工程计算领域长期处于统治地位，很多优秀的工程计算软件都是运用fortran语言编写，例如ANSYS,Marc,为了能够使用这些商业软件的高级功能，必须先学会fortran语言，才能编写应用程序接口。 　　在数值计算中，Fortran语言仍然不可替代。Fortran90标准引入了数组计算等非常利于矩阵运算的功能。在数组运算时，Fortran能够自动进行并行运算，这是很多编程语言不具备的。运用fortran语言，你能够运用很多现成的函数软件包，所以非常便利。（目前流行的Matlab的早期版本，主要就是为两个著名的Fortran函数包提供程序接口）

1、打开visualstudio2013，文件-新建-项目。2、在intelvisualfortran下面，点击emptyproject，并根据自己的设置，将命名以及存储位置设置好。3、创建好，在右侧就会出来一个项目，然后点击添加-新建项。4、此时即可以输入相调试的fortran程序，直接输入。5、进行调试(F5)或直接运行(Ctrl+F5)。

有免安装的，直接用就可以了，如fps4

明知故问。鄙视楼主Program www_fcode_cn implicit none real a(5,5) write(*,*) size(a)End Program www_fcode_cn这个程序会输出 25，你懂的

matlab是解释性语言，速度非常慢，因此大量计算适合用fortran。matlab的工具箱非常多，简单研究和作图比较合适。

你已的编译器已经能够完全正常运行，而且程序也可以运行了呢，只是你不知道如何运行它而已！在unix或linux下编译fortran，一般系统没有自带fortran的编译器，可以自行安装g77或gfortran编译器，它们是gcc一个系列的，也可以下载单独的g95编译器。假定你的fortran源程序为hello.for则输入如下的命令来编译它（g77、gfortran、g95）g77hello.for这样就会在当前目录生成a.out的可执行文件，这样既可执行./a.out如果需要指定生成的文件名，可以加个参数-o来制定文件名，像这样g77-ohellohello.for./hello明白了吗？unix类操作系统的执行档，并不需要.exe这样的扩展名的。

不是因为功能上的原因，而是语言本身的性质决定的。FORTRAN本身跟机器语言十分贴近，也就是说其翻译成机器语言的效率极高。另外，FORTRAN设计之初就是为数值计算而设计的，它的很多语法跟机器的数值交互要比C或者CPP效率高得多，有说法就会FORTRAN是原生支持数值计算的高级语言。由于长期以来在数值计算领域FORTRAN占据了统治地位，很多知名的程序的代码都是用FORTRAN写成的，另外由于F90的出现和F2003的发展，FORTRAN和特性进一步发展，交互性进一步增强，其擅长数值计算的特点能够进一步加强。总体来说，当下的数值计算的程序语言还是首选FORTRAN。

第一种方法： 1. 建立工作目录/ ]"2. 将ABAQUS安装目录6.4-pr11site下的aba_param_dp.inc 或 aba_param_sp.inc拷贝到工作目录，并改名为aba_param.inc；3. 将编译的fortran程序拷贝到工作目录；4. 将.obj文件拷贝到工作目录；5. 建立好输入文件.inp；6. 运行abaqus job=inp_name user=fortran name即可。第二种方法：在Job模块里，创建工作，在Edit Job对话框中选择General选项卡，在User subroutine file中点击Select按钮，从弹出对话框中选择你要调用的子程序文件（后缀为.for或.f）。

FORTRAN=FORmula TRANslator

Fortran 主要的编译有1. gfortran 免费的编译器，Linux 和 Windows 下都可以用，GCC 众多编译器的一个，我主要用的就是这个。2. G95 免费的编译器，安装文件很小的编译器，相比gfortran 优势不明显。3. Watcom 的 Fortran 编译器， 现在是免费的，Watcom 跟 GCC 一样有很多语言的编译器，Watcom 的 C 编译器不错，Fortran 没感觉多好。4. Intel Visual Fortran（简称 IVF），原来叫 Compaq Visual Fortran（简称 CVF） 收费软件，不过可以先下载下来试用 30 天，编译出来的二进制文件针对 Intel CPU 优化过，所以会快一些，但不是很明显。5. PGI Visual Fortran (PVF), 是 Portland Group Inc. 出的，收费软件。另外还有 Salford FTN95，simply fortran 等。 推荐使用 gfortran 和 Intel Visual Fortran.

A

FORTRAN语言是世界上第一个被正式推广使用的高级语言。它是1954年被提出来的，1956年开始正式使用，至今已有五十多年的历史，但仍历久不衰，它始终是数值计算领域所使用的主要语言。

实现并行的方式有很多种从硬件上分为 CPU 并行和GPU 并行CPU 并行目前主要有：1. 进程级适用于分布内存的MPI 并行(参考书 百度 都志辉) 2. 线程级别的openmp （参考书：百度 fortran openmp） 需要注意的是 MPI 相对底层一点，openmp 则非常简单 基本上1个月能够将你的程序并行化。 但前者的效率要高于后者 当然我是说你程序编写合理的情况下。GPU 并行：目前的显卡基本上都支持GPU并行，只要你的显卡不是很老，都可以用来做GPU 并行计算。 GPU计算说通俗一点就是借助显卡的处理核心来做科学计算，有效发挥显卡的多核计算能力如果你要进行GPU 编程，那么建议学习CUDA 对应于C 有cuda C 对应于fortran 有cuda fortran 当然 目前cuda fortran 的编译器只有pgi， 你可以很方便在网上下载到试用版本。就这么多吧，仅供参考。有什么问题可以在群里 或者在论坛里面提问，大家来讨论。最后还补充一下，如果你有正版PGI的话，也可以用openacc 也是一种简化的cuda 大体上就是你在代码中加入一些编译指导语句，然后编译器帮你把代码cuda 并行化。相当于cpu 并行里的openmp。

用.GT.等表示比较（关系运算)是FORTRAN77的用法，两边的点号与字母是一个整体，在Fortran90中可以用新的符号来表示，建议以后都用新方法。对照如下，括号里是新符号：.LT. (< ) .LE.(<=) .EQ. (=) .NE. (/=，不等于) .GE.(>=) .GT.(>)最近编写F77到F90的转换程序，发现有的将 .GT. 中的点号与字母分开，甚至有的在中间换行到续行中，为了处理这种情况，费了不少脑筋。

没有

（1）在求矩阵的范数之前，我们首先要清楚我们要求得是那一类矩阵范数，通常我们常用的矩阵范数可以分为：1范数，2范数，无穷范数，和Frobenius范数。具体的范数表示形式如下图所示：（2）上面介绍了几种常用的范数表示形式了，那么下面来看下怎么求具体的范数值。当然，我们可以根据定义来求每个范数的值，这样只针对于矩阵维度较小的矩阵适用，下面我们来看下当矩阵维数较大时我们怎么通过matlab来求矩阵的不同范数。（3）首先，我们来看下矩阵的1范数怎么通过matlab来求。先在matlab命令窗口中生成一个希尔伯特矩阵a=hilb(4)，然后，在命令窗口中输入nm1=norm(a,1) ，其中norm就是求矩阵范数的函数，1表示的是1范数。程序运行结果如下图所示，显然红色圈中部分就是所求的结果对应的列。（4）其次，看下怎么求矩阵的2范数。先在matlab命令窗口中生成一个希尔伯特矩阵a=hilb(4)，然后，在命令窗口中输入nm2=norm(a,2) ，其中norm就是求矩阵范数的函数，2表示的是2范数。程序运行结果如下图所示，当然这里不能向1范数那样，一眼看出结果。（5）下面看下怎么求矩阵的无穷范数。（相信聪明的同学已经想到了）先在matlab命令窗口中生成一个希尔伯特矩阵a=hilb(4)，然后，在命令窗口中输入nm3=norm(a,inf) ，其中norm就是求矩阵范数的函数，inf表示的是无穷范数。程序运行结果如下图所示。步骤阅读6（6）最后我们看下怎么求矩阵的Frobenius范数。先在matlab命令窗口中生成一个希尔伯特矩阵a=hilb(4)，然后，在命令窗口中输入nm4=norm(a,"for") ，其中norm就是求矩阵范数的函数，for表示的是Frobenius范数,就是前三个字母嘛。程序运行结果如下图所示。至此几种常用的矩阵范数都求出来了，大家可以试试了！

Fortran自带MOD函数可以解决你的问题MOD（x,y）就是x/y的余数

有的。fortran最新标准是2012年发布的fortran2008。fortran最新的编译器目前还在更新并作商业化发布，比如 Intel fortran。在水力计算、气象数值预报、航天航空领域有着大量应用。

real*8 :: a(n)real*8 :: tempinteger :: i,jdo i=1,n-1 do j=i+1,n if (a(i) .lt. a(j)) then temp = a(i) a(i) = a(j) a(j) = temp endif enddoenddo

DO I=1,m DO J=1,n sum=sum+A(I,J) B(I)=sum ENDDO ENDDO这个应该是求每行的和放到B数组里吧，你这用的临时量sum在每行求完后没有清零

在fortran语言中，abs可以输入的参数类型为 integer, real, or complex；iabs只能输入 INTEGER(4) 类型 ，返回类型 INTEGER(4)http://zhidao.baidu.com/question/327752376.html

fortran是一种程序设计语言

解决方法很多!第一，V（5000）=0的意思是，第5000个元素为0，而前4999个没有进行赋值。如果写为：V(1:5000)=0，这就表示5000个全部为0了。第二，V（0:5000）=1,2,3这样写是非法的。第三，V(5000)也可以写为V(1:5000)表示V数组里有5000个数。如果省略了冒号，则默认从1开始如果写成V(0:5000)则表示V数组有5001个数。请采纳答案，支持我一下。

Fortran Coder 论坛：fcode.cn含工具，资源，电子书，代码共享，算法研讨等。bbs.pfan.cn 访问老有问题，现在人气不如以前了。此外，仿真论坛，气象论坛也有Fortran版块。地址可自己百度。

double randomz (int ia, int ib){double x; //返回值static int initial[15]={0};static double iz,iy[15]={0.0};//使用static类型，为下次调用保留值，不然每次都要从头开始switch(ia){ case 1: //ia参数为1，从键盘输入种子初始化随机数生成器 iz=100000001.0 read(5,"(5i8)") initial 这句话直接翻成C很难，//意思就是从键盘(5号设备)连续读入5个整数（难道不是15？// 我怀疑源码写错了，应该是15i8），存入整型数组 initial ；//而且Fortran的格式描述符i8强制每个整数都是8位（千万位），//如果输入不足8位就切换到下一个数进行输入，则Fortran会//自动在不足8位的数的右侧补上足够的0，以放大到千万位 iy=dble(initial) 这句话直接翻成C很难，//本句用到Fortran90/95特色的数组整体操作，C/C++要用循环；//是把整型数组initial逐项复制给double数组iy，转换成double型 x = iy[ib]*1.0E-8 ;case 2: //ia参数为2，继续使用已经初始化好了的生成器 iy[ib] = 329.0 * iy[ib] % iz ; x = iy[ib]*1.0E-8 ;case 3: //ia参数为3，重新初始化，但还使用原先的种子 iy=dble(initial) 这句话直接翻成C很难，用循环完成。}//switch结束return x;}=================================算法的主要思想就是“线性同余法”，linear-congruential method其基本迭代公式为X[n+1] = ( A * X[n] + B ）% CX的初始值随便取在你给的源代码里：A= 329.0B= 0C= 100000001.0源代码中的关键是 iy[ib] = 329.0 * iy[ib] % iz另外 x = iy[ib]*1.0E-8 是为了将结果归一化到0~1之间再返回你可以自己手工算几个数，就能看出这个算法的奥妙了。另外需要指出的是——源代码里用static就是为了每次case 2时候的调用，都是在对上一次的结果进行迭代。而ib参数的用处就是保持有几组不同的独立迭代序列可用，防止不同用途的几处“生成伪随机数”调用互相干扰。

自然对数是 log(x)十为第的对数是 log10(x)其他底的对数，用换底公式。

复数运算的速度程序参数与字串内存的动态管理多维阵列的处理函数调用与参数传递

FORTRAN是Formula Translation的缩写，意为“公式翻译”。它是为科学、工程问题或企事业管理中的那些能够用数学公式表达的问题而设计的，其数值计算的功能较强。FORTRAN语言问世以来，根据需要几经发展，先后推出了不同的版本，其中最流行的是1958年出现的FORTRANⅡ和1962年出现的FORTRANⅣ。1966年美国标准化协会(ANSI)公布了两个美国标准文本：·标准FORTRAN（X3.9-1966）。大致相当于FORTRANⅣ。·标准基本FORTRAN（X3.10-1966）大致相当FORTRANⅡ。1972年国际标准化组织(ISO)接受了美国标准，在稍加修改后公布了ISO FORTRAN标准，即《程序设计语言FORTRAN ISO 1539-1972》，它分为三级，即：· 完全的(一级)FORTRAN，相当于FORTRANⅣ。· 中间的(二级)FORTRAN，介于FORTRANⅡ和FORTRANⅣ之间。· 基本的(三级)FORTRAN，相当于FORTRANⅡ。FORTRANⅣ(即FORTRAN 66)流行了十几年，几乎统治了所有的数值计算领域，许多应用程序和程序库都是用FORTRANⅣ语言编写的。美国标准化协会(ANSI)在1976年对ANSI FORTRAN(X3.9-1966)进行了修订，预定在1977年通过，为了区别于FORTRAN 66，新标准定名为FORTRAN 77。实际上到1978年4月才由ANSI正式公布作为新的美国国家标准。即FORTRAN(X3.9-1978)。1980年，FORTRAN 77被接受为国际标准，即《程序设计语言FORTRAN ISO 1539-1980》，该标准分为全集和子集。中国制订的FORTRAN标准，基本上采用了国际标准，于1983年5月公布执行，标准号为GB3057-82。FORTRAN 77标准完成后，新版本的修订工作也在同一时间开始进行。这个版本进行了15年，最后在1992年正式由国际标准组织ISO公布，它就是FORTRAN 90。FORTRAN 90对以往的FORTRAN语言标准作了大量的改动，使之成为一种功能强大、具有现代语言特征的计算机语言。其主要特色是加入了面向对象的概念及工具、提供了指针、加强了数组的功能、改良了旧式FORTRAN语法中的编写“版面”格式。FORTRAN 95标准在1997年同样由ISO公布，它可以视为是FORTRAN 90的修正版，主要加强了FORTRAN在并行运算方面的支持。同时一些公司纷纷推出Visual Fortran，这为工程技术界进行科学计算和编写面向对象的工程实用软件的用户提供了极大的方便。熟悉VB或VC的读者可以很容易地掌握Visual Fortran的使用，进一步开发出自己专业领域的Windows下的界面友好的工程应用软件。

科学计算。fortran是数值计算领域所使用的主要语言，主要应用于科学计算，广泛应用于科学及工程计算应用程序的开发。FORTRAN语言是FormulaTranslation的缩写，意为“公式翻译”，它是为科学、工程问题或企事业管理中的那些能够用数学公式表达的问题而设计的，其数值计算的功能较强。

当然有必要了。Fortran本身就是用来做科学计算的。目前依然存在除了最起码的原因历史上有庞大的优秀的fortran代码，可以直接拿来使用，如果再也C等语言在写一遍会超级复杂之外，最根本的原因，就是他执行效率高。它原生支持矩阵计算，这在科学计算高维计算中尤其重要和方便，原生支持复数运算。FORTRAN语言是Formula Translation的缩写，意为“公式翻译”。它是为科学、工程问题或企事业管理中的那些能够用数学公式表达的问题而设计的，其数值计算的功能较强。FORTRAN语言是世界上第一个被正式推广使用的高级语言。它是1954年被提出来的，1956年开始正式使用，直到2021年已有六十年的历史，但仍历久不衰，它始终是数值计算领域所使用的主要语言。

FORTRAN语言是世界上第一个被正式推广使用的高级语言。它是1954年被提出来的，1956年开始正式使用，Fortran语言的发明人John W. Backus。

FORTRAN语言问世以来，根据需要几经发展，先后推出了不同的版本，其中最流行的是1958年出现的FORTRANⅡ和1962年出现的FORTRANⅣ。1966年美国标准化协会(ANSI)公布了两个美国标准文本：·标准FORTRAN（X3.9-1966）。大致相当于FORTRANⅣ。·标准基本FORTRAN（X3.10-1966）大致相当FORTRANⅡ。1972年国际标准化组织(ISO)接受了美国标准，在稍加修改后公布了ISO FORTRAN标准，即《程序设计语言FORTRAN ISO 1539-1972》，它分为三级，即：· 完全的(一级)FORTRAN，相当于FORTRANⅣ。· 中间的(二级)FORTRAN，介于FORTRANⅡ和FORTRANⅣ之间。· 基本的(三级)FORTRAN，相当于FORTRANⅡ。FORTRANⅣ(即FORTRAN 66)流行了十几年，几乎统治了所有的数值计算领域，许多应用程序和程序库都是用FORTRANⅣ语言编写的。美国标准化协会(ANSI)在1976年对ANSI FORTRAN(X3.9-1966)进行了修订，预定在1977年通过，为了区别于FORTRAN 66，新标准定名为FORTRAN 77。实际上到1978年4月才由ANSI正式公布作为新的美国国家标准。即FORTRAN(X3.9-1978)。1980年，FORTRAN 77被接受为国际标准，即《程序设计语言FORTRAN ISO 1539-1980》，该标准分为全集和子集。中国制订的FORTRAN标准，基本上采用了国际标准，于1983年5月公布执行，标准号为GB3057-82。FORTRAN 77标准完成后，新版本的修订工作也在同一时间开始进行。这个版本进行了15年，最后在1992年正式由国际标准组织ISO公布，它就是FORTRAN 90。FORTRAN 90对以往的FORTRAN语言标准作了大量的改动，使之成为一种功能强大、具有现代语言特征的计算机语言。其主要特色是加入了面向对象的概念及工具、提供了指针、加强了数组的功能、改良了旧式FORTRAN语法中的编写“版面”格式。FORTRAN 95标准在1997年同样由ISO公布，它可以视为是FORTRAN 90的修正版，主要加强了FORTRAN在并行运算方面的支持。同时一些公司纷纷推出Visual Fortran，这为工程技术界进行科学计算和编写面向对象的工程实用软件的用户提供了极大的方便。熟悉VB或VC的读者可以很容易地掌握Visual Fortran的使用，进一步开发出自己专业领域的Windows下的界面友好的工程应用软件。

找到你的源代码对应的 fortran 版本（for77、for90、for95）、对应平台（win、unix、...）的编译器、连接程序；用它们将你的源代码转换为可执行文件；再运行这个可执行文件。

%是FORTRAN90新增加的字符，表示结构中成员的一般形式：结构名％成员名。即在所属结构名后写一个百分号（％）而后写出成员本身名称。这样的成员可以像访问变量一样被访问，包括赋值、打印、引用等。

看你虚参的定义。虚参定义从哪儿开始，就从哪儿开始。

分类: 电脑/网络 >> 程序设计 >> 其他编程语言 问题描述: 如0.25*(x2)-2 >=如何用fortran表示写出来?谢谢 解析: FORTRAN中大于等于用 (.ge.)表示，例如X>=Y可以表示为 （X.GE.Y）

像是fortran77的程序，又不完全是。改为fortran77的程序如下，编译、连接均通过。 PROGRAM JOURNAL DIMENSION H(61,21),P(61,21) DATA B,R,C0,AN,EDA,EPSON,aSUM1/6E-2,25E-3,5E-5,6E4,0.05,0.7,0/ OPEN(9,FILE="PRESSURE.DAT",STATUS="new") OPEN(8,FILE="FILM.DAT",STATUS="new") PI=3.1415926 N=61 M=21 DX=2.0*PI/FLOAT(N-1) DY=1./FLOAT(M-1) OMEGA=AN*2.0*PI/60.0 U=OMEGA*R ALFA=(R/B*DX/DY)**2 CALL SUBH(N,M,DX,EPSON,H) CALL SUBP(N,M,DX,EPSON,ALFA,H,P) CALL SUBM(N,M,DX,DY,aSUM1,PI,P) CALL OUTPUT(N,M,DX,DY,H,P) STOP END SUBROUTINE SUBH(N,M,DX,EPSON,H) DIMENSION H(N,M) DO 1 I=1,N SETA=(I-1.0)*DX DO 2 J=1,M H(I,J)=1.0+EPSON*COS(SETA)2 continue1 continue RETURN END SUBROUTINE SUBP(N,M,DX,EPSON,ALFA,H,P) DIMENSION H(N,M),P(N,M) DO 3 I=1,N DO 4 J=2,M-1 P(I,J)=0.54 continue 3 continue DO 5 J=1,M P(1,J)=0.0 P(N,J)=0.05 continue DO 6 I=1,N P(I,1)=0.0 P(I,M)=0.06 continue IK=010 C1=0.0 ALOAD=0.0 DO 12 I=2,N-1 I1=I-1 I2=I+1 DO 13 J=2,M-1 PD=P(I,J) J1=J-1 J2=J+1 A1=(0.5*(H(I2,J)+H(I,J)))**3 A2=(0.5*(H(I,J)+H(I1,J)))**3 A3=ALFA*(0.5*(H(I,J2)+H(I,J)))**3 A4=ALFA*(0.5*(H(I,J)+H(I,J1)))**3 P(I,J)=(-DX*(H(I2,J)-H(I1,J))+A1*P(I2,J)+A2*P(I1,J)+A3*P(I,J2)+ 1 A4*P(I,J1))/(A1+A2+A3+A4) P(I,J)=0.7*PD+0.3*P(I,J) IF(P(I,J).LT.0.0)P(I,J)=0.0 C1=C1+ABS(P(I,J)-PD) ALOAD=ALOAD+P(I,J)20 CONTINUE13 continue12 continue IK=IK+1 C1=C1/ALOAD WRITE(*,*)IK,C1,ALOAD IF(C1.GT.1.E-7)GOTO 10 RETURN END SUBROUTINE SUBM(N,M,DX,DY,aSUM1,PI,P) DIMENSION P(N,M) PX=0 PY=0 TEMP=PI/60 DO 30 I=1,N AI=(I-1)*TEMP DO 35 J=1,M PX=PX-P(I,J)*COS(AI)*DX*DY PY=PY+P(I,J)*SIN(AI)*DX*DY35 continue30 continue aSUM1=SQRT(PX*PX+PY*PY) WRITE(*,*)aSUM1 RETURN END SUBROUTINE OUTPUT(N,M,DX,DY,H,P) DIMENSION Y(21),H(N,M),P(N,M) DO 36 J=1,M Y(J)=(J-1.)*DY-0.536 continue WRITE(8,40)Y(1),(Y(J),J=1,M) WRITE(9,40)Y(1),(Y(J),J=1,M) DO 37 I=1,N AX=(I-1.0)*360.0/(N-1.0) WRITE(8,40)AX,(H(I,J),J=1,M) WRITE(9,40)AX,(P(I,J),J=1,M)37 continue40 FORMAT(22(E12.6,1X)) STOP END

用.GT.等表示比较（关系运算)是FORTRAN77的用法，两边的点号与字母是一个整体，在Fortran90中可以用新的符号来表示，建议以后都用新方法。对照如下，括号里是新符号：.LT. (< ) .LE.(<=) .EQ. (=) .NE. (/=，不等于) .GE.(>=) .GT.(>)最近编写F77到F90的转换程序，发现有的将 .GT. 中的点号与字母分开，甚至有的在中间换行到续行中，为了处理这种情况，费了不少脑筋。

这是C++的用法看++在哪个位置是加法运算符的缩写譬如i=i+1可以写成i++可以是先赋值再运算也可以是先运算在赋值再看看别人怎么说的。

Fortran 主要的编译有1. gfortran 免费的编译器，Linux 和 Windows 下都可以用，GCC 众多编译器的一个，我主要用的就是这个。2. G95 免费的编译器，安装文件很小的编译器，相比gfortran 优势不明显。3. Watcom 的 Fortran 编译器， 现在是免费的，Watcom 跟 GCC 一样有很多语言的编译器，Watcom 的 C 编译器不错，Fortran 没感觉多好。4. Intel Visual Fortran（简称 IVF），原来叫 Compaq Visual Fortran（简称 CVF） 收费软件，不过可以先下载下来试用 30 天，编译出来的二进制文件针对 Intel CPU 优化过，所以会快一些，但不是很明显。5. PGI Visual Fortran (PVF), 是 Portland Group Inc. 出的，收费软件。另外还有 Salford FTN95，simply fortran 等。 推荐使用 gfortran 和 Intel Visual Fortran.

dowhile表示一直做某事/=表示不等于意思是一直做下面的代码知道age/=k为止或者理解为当age/=k就执行代码

FORTRAN语言是世界上第一个被正式推广使用的高级语言。它是1954年被提出来的，1956年开始正式使用，至今已有三十多年的历史，但仍历久不衰，它始终是数值计算领域所使用的...

fortran的强项在于有长期的积累，长于数值计算、并行计算；弱项是其界面功能、交互功能、图形功能稍弱。

第一个问题，虚参有9个，实参却有10个。第二个问题，应该是第一个问题引起的。解决了第一个，第二个应该就没了。我猜，虚参里是否少写了一个 DSOL ？