Fortran语言
#LOCAL REAL Uds,Uqs,angles,vs,w,o,ang,ta,tb,t0,Ud1,Uq1,E
#LOCAL REAL tc,t,mt,pi
#LOCAL INTEGER sign,M,nl,n
t=TIME
tc=$TC
nl=DNINT(tc/delt)
pi=3.1415926
E=$Udc
w=2*pi*50
if(t.lt.delt) then
sign=0
M=0
Uds=0
Uqs=0
angles=0
vs=0
o=0
ang=0
ta=0
tb=0
t0=0
Ud1=0
Uq1=0
n=1
end if
if(mod(n,nl).lt.1E-7) then
!初始化变量
M=M+1
Uds=$Ud
Uqs=$Uq
angles=$angle
!变换到alphabeta坐标直角坐标
Ud1=cos(angles)*Uds+sin(angles)*Uqs
Uq1=sin(angles)*Uds-cos(angles)*Uqs
! Ud1=$Ud
! Uq1=$Uq
!生成矢量的模和角度
vs=sqrt(Ud1**2+Uq1**2)
if(Ud1.gt.0.and.Uq1.ge.0) then
ang=atan(Uq1/Ud1)
else
ang=atan(Uq1/Ud1)+2*pi
end if
if(Ud1.lt.0) then
ang=atan(Uq1/Ud1)+pi
end if
if(Ud1.eq.0) then
if(Uq1.ge.0) then
ang=pi/2
else
ang=3*pi/2
end if
end if
!判断矢量的所属,计算角度
if(ang.ge.0.and.ang.lt.pi/3) then
sign=1
o=ang
! o=w*t
end if
if(ang.ge.pi/3.and.ang.lt.2*pi/3) then
sign=2
o=ang-pi/3
! o=w*t-pi/3
end if
if(ang.ge.2*pi/3.and.ang.lt.3*pi/3) then
sign=3
o=ang-2*pi/3
! o=w*t-2*pi/3
if(ang.ge.3*pi/3.and.ang.lt.4*pi/3) then
sign=4
o=ang-3*pi/3
! o=w*t-3*pi/3
end if
if(ang.ge.4*pi/3.and.ang.lt.5*pi/3) then
sign=5
o=ang-4*pi/3
! o=w*t-4*pi/3
end if
if(ang.ge.5*pi/3.and.ang.lt.6*pi/3) then
sign=6
o=ang-5*pi/3
! o=w*t-5*pi/3
end if
!计算ta,tb,to
ta=sqrt(2.0)*vs*tc*sin(pi/3-o)/E
tb=sqrt(2.0)*vs*tc*sin(o)/E
if((ta+tb).ge.tc) then
ta=tc*ta/(ta+tb)
tb=tc*tb/(ta+tb)
end if
t0=tc-ta-tb
!减小误差
!end 111
end if
!23456
! 选择开关变量
if(sign.eq.0) then
p1=0
p2=0
p3=0
p4=0
p6=0
end if
if(sign.eq.1) then
!前半采样周期
if(t.lt.(t0/4+M*tc).and.t.ge.M*tc) then
!vo
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
if(t.ge.(t0/4+M*tc).and.t.lt.(t0/4+ta/2+M*tc)) then
!v1
p1=1
p2=1
p3=0
p4=0
p5=0
p6=1
end if
if(t.ge.(t0/4+ta/2+M*tc).and.t.lt.(t0/4+ta/2+tb/2+M*tc)) then ! %v2
p1=1
p2=1
p3=1
p4=0
p5=0
p6=0
end if
if(t.ge.(t0/4+ta/2+tb/2+M*tc).and.t.lt.(tc/2+M*tc)) then
!%v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
!%后半采样周期
if(t.ge.(tc/2+M*tc).and.t.lt.(tc/2+t0/4+M*tc)) then
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
if(t.ge.(tc/2+t0/4+M*tc).and.t.lt.(tc/2+t0/4+tb/2+M*tc)) then
! %V2
p1=1
p2=1
p3=1
p4=0
p5=0
p6=0
end if
if(t.ge.(tc/2+t0/4+tb/2+M*tc).and. t.lt.(tc/2+t0/4+tb/2+ta/2+M*tc)) then !%v1
p1=1
p2=1
p3=0
p4=0
p5=0
p6=1
end if
if(t.ge.(tc/2+t0/4+tb/2+ta/2+M*tc).and.t.lt.(tc+M*tc)) then
!%v0
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
end if
!23456
if(sign.eq.2) then
!前半采样周期
if(t.lt.(t0/4+M*tc).and.t.ge.M*tc) then
!vo
p1=0
p3=0
p4=1
p5=0
p6=1
end if
if(t.ge.(t0/4+M*tc).and.t.lt.(t0/4+ta/2+M*tc)) then
!v2 p1(n)=1;p2(n)=1;p3(n)=1;p4(n)=0;p5(n)=0;p6(n)=0;
p1=1
p2=1
p3=1
p4=0
p5=0
p6=0
end if
if(t.ge.(t0/4+ta/2+M*tc).and.t.lt.(t0/4+ta/2+tb/2+M*tc)) then ! %v3 p1(n)=0;p2(n)=1;p3(n)=1;p4(n)=1;p5(n)=0;p6(n)=0;v3
p1=0
p2=1
p3=1
p4=1
p5=0
p6=0
end if
if(t.ge.(t0/4+ta/2+tb/2+M*tc).and.t.lt.(tc/2+M*tc)) then !%v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
!%后半采样周期
if(t.ge.(tc/2+M*tc).and.t.lt.(tc/2+t0/4+M*tc)) then ! %v7
p1=1
p2=0
p3=1
p4=0
p6=0
end if
if(t.ge.(tc/2+t0/4+M*tc).and.t.lt.(tc/2+t0/4+tb/2+M*tc)) then
! %V3 p1(n)=0;p2(n)=1;p3(n)=1;p4(n)=1;p5(n)=0;p6(n)=0;v3
p1=0
p2=1
p3=1
p4=1
p5=0
p6=0
end if
if(t.ge.(tc/2+t0/4+tb/2+M*tc).and. t.lt.(tc/2+t0/4+tb/2+ta/2+M*tc)) then !%v2 p1(n)=1;p2(n)=1;p3(n)=1;p4(n)=0;p5(n)=0;p6(n)=0;
p1=1
p2=1
p3=1
p4=0
p5=0
p6=0
end if
if(t.ge.(tc/2+t0/4+tb/2+ta/2+M*tc).and.t<(tc+M*tc)) then
!%v0
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
end if
if(sign.eq.3) then
!前半采样周期
if(t.lt.(t0/4+M*tc).and.t.ge.M*tc) then
!vo
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
if(t.ge.(t0/4+M*tc).and.t.lt.(t0/4+ta/2+M*tc)) then
!v3 p1(n)=0;p2(n)=1;p3(n)=1;p4(n)=1;p5(n)=0;p6(n)=0;v3
p1=0
p2=1
p3=1
p4=1
p5=0
p6=0
end if
if(t.ge.(t0/4+ta/2+M*tc).and.t.lt.(t0/4+ta/2+tb/2+M*tc)) then
! %v4
p1=0
p2=0
p3=1
p4=1
p5=1
p6=0
end if
if(t.ge.(t0/4+ta/2+tb/2+M*tc).and.t.lt.(tc/2+M*tc)) then !%v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
!%后半采样周期
if(t.ge.(tc/2+M*tc).and.t.lt.(tc/2+t0/4+M*tc)) then ! %v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
if(t.ge.(tc/2+t0/4+M*tc).and.t.lt.(tc/2+t0/4+tb/2+M*tc)) then
! %V4
p1=0
p2=0
p3=1
p4=1
p5=1
p6=0
end if
if(t.ge.(tc/2+t0/4+tb/2+M*tc).and. t.lt.(tc/2+t0/4+tb/2+ta/2+M*tc)) then !%v3 p1(n)=0;p2(n)=1;p3(n)=1;p4(n)=1;p5(n)=0;p6(n)=0;v3
p1=0
p2=1
p3=1
p4=1
p5=0
p6=0
end if
if(t.ge.(tc/2+t0/4+tb/2+ta/2+M*tc).and.t.lt.(tc+M*tc)) then
!%v0
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
end if
if(sign.eq.4) then
!前半采样周期
if(t.lt.(t0/4+M*tc).and.t.ge.M*tc) then
!vo
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
if(t.ge.(t0/4+M*tc).and.t.lt.(t0/4+ta/2+M*tc)) then
!v4
p1=0
p2=0
p3=1
p4=1
p5=1
p6=0
end if
if(t.ge.(t0/4+ta/2+M*tc).and.t.lt.(t0/4+ta/2+tb/2+M*tc)) then
! %v5 p1(n)=0;p2(n)=0;p3(n)=0;p4(n)=1;p5(n)=1;p6(n)=1;v5
p1=0
p2=0
p3=0
p4=1
p5=1
p6=1
end if
if(t.ge.(t0/4+ta/2+tb/2+M*tc).and.t.lt.(tc/2+M*tc)) then
!%v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
!%后半采样周期
if(t.ge.(tc/2+M*tc).and.t.lt.(tc/2+t0/4+M*tc)) then
! %v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
if(t.ge.(tc/2+t0/4+M*tc).and.t.lt.(tc/2+t0/4+tb/2+M*tc)) then
! %V5
p1=0
p2=0
p3=0
p4=1
p5=1
p6=1
end if
if(t.ge.(tc/2+t0/4+tb/2+M*tc).and. t.lt.(tc/2+t0/4+tb/2+ta/2+M*tc)) then !%v4 p1(n)=0;p2(n)=0;p3(n)=1;p4(n)=1;p5(n)=1;p6(n)=0;v4
p1=0
p2=0
p4=1
p5=1
p6=0
end if
if(t.ge.(tc/2+t0/4+tb/2+ta/2+M*tc).and.t.lt.(tc+M*tc)) then !%v0
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
end if
if(sign.eq.5) then
!前半采样周期
if(t.lt.(t0/4+M*tc).and.t.ge.M*tc) then
!vo
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
if(t.ge.(t0/4+M*tc).and.t.lt.(t0/4+ta/2+M*tc)) then
!v5
p1=0
p2=0
p3=0
p4=1
p5=1
p6=1
end if
if(t.ge.(t0/4+ta/2+M*tc).and.t.lt.(t0/4+ta/2+tb/2+M*tc)) then ! %v6
p1=1
p3=0
p4=0
p5=1
p6=1
end if
if(t.ge.(t0/4+ta/2+tb/2+M*tc).and.t.lt.(tc/2+M*tc)) then
!%v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
!%后半采样周期
if(t.ge.(tc/2+M*tc).and.t.lt.(tc/2+t0/4+M*tc)) then
! %v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
if(t.ge.(tc/2+t0/4+M*tc).and.t.lt.(tc/2+t0/4+tb/2+M*tc)) then
! %V6
p1=1
p2=0
p3=0
p4=0
p5=1
p6=1
end if
if(t.ge.(tc/2+t0/4+tb/2+M*tc).and. t.lt.(tc/2+t0/4+tb/2+ta/2+M*tc)) then !%v5
p1=0
p2=0
p3=0
p4=1
p5=1
p6=1
end if
if(t.ge.(tc/2+t0/4+tb/2+ta/2+M*tc).and.t.lt.(tc+M*tc)) then
!%v0
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
end if
if(sign.eq.6) then
!前半采样周期
if(t.lt.(t0/4+M*tc).and.t.ge.M*tc) then
!vo
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
if(t.ge.(t0/4+M*tc).and.t.lt.(t0/4+ta/2+M*tc)) then
!v6
p1=1
p2=0
p3=0
p4=0
p5=1
p6=1
end if
if(t.ge.(t0/4+ta/2+M*tc).and.t.lt.(t0/4+ta/2+tb/2+M*tc)) then ! %v1p1(n)=1;p2(n)=1;p3(n)=0;p4(n)=0;p5(n)=0;p6(n)=1
p1=1
p2=1
p3=0
p4=0
p5=0
p6=1
end if
if(t.ge.(t0/4+ta/2+tb/2+M*tc).and.t.lt.(tc/2+M*tc)) then !%v7
p2=0
p3=1
p4=0
p5=1
p6=0
end if
!%后半采样周期
if(t.ge.(tc/2+M*tc).and.t.lt.(tc/2+t0/4+M*tc)) then
! %v7
p1=1
p2=0
p3=1
p4=0
p5=1
p6=0
end if
if(t.ge.(tc/2+t0/4+M*tc).and.t.lt.(tc/2+t0/4+tb/2+M*tc)) then
! %V1 p1=1;p2(n)=1;p3(n)=0;p4(n)=0;p5(n)=0;p6(n)=1
p1=1
p2=1
p3=0
p4=0
p5=0
p6=1
end if
if(t.ge.(tc/2+t0/4+tb/2+M*tc).and. t.lt.(tc/2+t0/4+tb/2+ta/2+M*tc)) then !v6 p1(n)=1;p2(n)=0;p3(n)=0;p4(n)=0;p5(n)=1;p6(n)=1;v6
p1=1
p2=0
p3=0
p4=0
p5=1
p6=1
end if
if(t.ge.(tc/2+t0/4+tb/2+ta/2+M*tc).and.t.lt.(tc+M*tc)) then
!%v0
p1=0
p2=1
p3=0
p4=1
p5=0
p6=1
end if
$pulse1=p1
$pulse2=p2
$pulse3=p3
$pulse4=p4
$pulse5=p5
$pulse6=p6
n=n+1
! $pulse6=tb
! $pulse4=t0
! $pulse5=sin(pi/3-o) ! $pulse3=o
! $pulse2=ang
! $pulse1=sign
(完整)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语法手册
1 FORTRAN77四则运算符 + - * / ** (其中**表示乘方) 在表达式中按优先级次序由低到高为：+或-→*或/→**→函数→() 2 FORTRAN77变量类型 2.1 隐含约定：I-N规则 凡是以字母I，J，K，L，M，N六个字母开头的，即认为是整型变量，其它为实型变量。 如IMPLICIT REAL (I,J) 三种定义的优先级别由低到高顺序为：I-N规则→IMPLICIT语句→类型说明语句，因此，在程序中IMPLICIT语句应放在类型说明语句之前。 2.4 数组的说明与使用 使用I-N规则时用DIMENSION说明数组，也可在定义变量类型同时说明数组，说明格式为：数组名(下标下界,下标上界），也可省略下标下界，此时默认为1，例： DIMENSION IA(0:9),ND(80:99),W(3,2),NUM(-1:0),A(0:2,0:1,0:3) REAL IA(10),ND(80:99)使用隐含DO循环进行数组输入输出操作：例如WRITE(*,10) ('I=',I,'A=',A(I),I=1,10,2) 10FORMAT(1X,5(A2,I2,1X,A2,I4)) 2.5 使用DATA语句给数组赋初值 变量表中可出现变量名,数组名,数组元素名,隐含DO循环，但不许出现任何形式的表达式：例如 DATA A,B,C/-1.0,-1.0,-1.0/ DATA A/-1.0/,B/-1.0/,C/-1.0/ DATA A,B,C/3*-1.0/CHARACTER*6 CHN(10)
DATA CHN/10*' '/INTEGER NUM(1000) DATA (NUM(I),I=1,500)/500*0/,(NUM(I),I=501,1000)/500*1/ 3 FORTRAN77程序书写规则 程序中的变量名，不分大小写； 变量名称是以字母开头再加上1到5位字母或数字构成，即变更名字串中只有前6位有效； 一行只能写一个语句； 程序的第一个语句固定为PROGRAM 程序名称字符串 某行的第１个字符至第５个字符位为标号区，只能书写语句标号或空着或注释内容； 某行的第１个字符为Ｃ或＊号时，则表示该行为注释行，其后面的内容为注释内容； 某行的第６个字符位为非空格和非０字符时，则该行为上一行的续行，一个语句最多可有19个续行； 某行的第７至７２字符位为语句区，语句区内可以任加空格以求美观； 某行的第７３至８０字符位为注释区，８０字符位以后不能有内容。 4 FORTRAN77关系运算符 .GT. 大于 .GE. 天于或等于 .LT. 小于 .LE. 小于或等于 .EQ. 等于 .NE. 不等于 .AND. 逻辑与 .OR. 逻辑或 .NOT. 逻辑非 .EQV. 逻辑等 .NEQV. 逻辑不等 运算符优先级由高到低顺序为:()→**→*或/→+或-→.GT.或.GE.或.LT. 或.LE.或.EQ.或.NE.→.NOT.→.AND.→.OR.→.EQV.或.NEQV 5 FORTRAN77语句
(完整版)FORTRAN90用法总结,推荐文档
FORTRAN90用法总结 1、数据类型 整型 INTEGER 实型 REAL 复型 COMPLEX 字符型 CHARACTER 逻辑型 LOGICAL 2、特殊函数 1>定值函数 PARAMETER(变量列表) 或类型说明符，PARAMETER::变量列表 2>IMPLICIT 类型说明符（字母表）类型定义函数（默认I-N规则，取消IMPLICIT NONE） 3>** 乘方 4>STOP[标号]终止执行 5>PAUSE[标号]暂停程序 6>字符串长度函数 LEN（String）、 LEN_TRIM(String) 7>除去字符串尾部空格 TRIM(String) 8>字符串位置函数 INDEX(String1,String2) 9>字符串比较函数 LGE(String1,String2) LGT(String1,String2) LLE(String1,String2) LLT(String1,String2) 10>字符串空格首位调整 ADJUSTL(String)、ADJUSTR(String) 11>字符串、ASSII码互相转化函数 ①化为ASSII码 ICHAR(String)、IACHAR(String) ②化为字符串CHAR(I)、ACHAR(I) 12>.LT.(<) .LE.(<=) .EQ.(=) .NE.(/=) .GT.(>) .GE.(>=) 13>.AND.与 .OR.或 .NOT.非 .EQV.等于 .NEQV.不等 .XOR.异或 3、IF函数 1> 单支IF IF(逻辑表达式) THEN 块语句 END IF 2> 双支IF IF(逻辑表达式) THEN 块语句1 ELSE 块语句2 END IF 3>多支IF IF(逻辑表达式1) THEN 块语句1 ELSE IF(逻辑表达式2) THEN
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=）的语法格式声明字符变量的长度。
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)
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.360docs.net/doc/636011130.html,\MSDEV98\dfdev.exe Microsoft Developer Studio的界面如下图所示： 文档窗口 工作空间窗口 输出窗口
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语言) (总分：36.00，做题时间：90分钟) 一、 (总题数：36，分数：36.00) 1.编译程序能将高级语言编写的源程序转换成( )。 A．解释程序 B．汇编程序 C．映象程序 D．目标程序 （分数：1.00） A. B. C. D. √ 解析： 2.一个完整的FORTRAN源程序( )。 A．至少包括一个主程序 B．至少包括一个主程序和一个子程序 C．由一个主程序与一个以上的子程序组成 D．由一个主程序与一个子程序组成 （分数：1.00） A. √ B. C. D. 解析： 3.语句函数定义语句在程序内合法的位置是( )。 A．在程序块开头语句之后，END语句之前 B．在程序块开头语句之后，可执行语句之前 C．在说明语句之后，END语句之前 D．在说明语句之后，可执行语句之前 （分数：1.00） A. B. C. D. √ 解析： 4.下列关于函数子程序虚实参数的错误说法是( )。 A．可以没有形参数 B．虚实结合的数组长度可以不同 C．实参表与虚参表类型可以不同 D．函数名可以作为虚参
（分数：1.00） A. B. C. √ D. 解析： 5.下列叙述中正确的是( )。 A．FORTRAN程序块中，无名公用语句只能有一个B．FORTRAN子程序中，至少应有一个RETURN语句C．FORTRAN程序块中，最后一行必须是END语句D．FORTRAN程序块中，必须有变量说明语句 （分数：1.00） A. B. C. √ D. 解析： 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．循环参数设置错误 （分数：1.00） A. B. C. D. √ 解析： 7.下列数据中，不符合FORTRAN常量表示法的是( )。 A．-25.6 B．2.758D3 C．'FOOT"=' D．TRUE （分数：1.00） A. B. C. D. √ 解析：
fortran语法手册
1F O R T R A N77四则运算符+ - */ ** (其中**表示乘方) 在表达式中按优先级次序由低到高为：+或-→*或/→**→函数→() 2 FORTRAN77变量类型 隐含约定：I-N规则 凡是以字母I，J，K，L，M，N六个字母开头的，即认为是整型变量，其它为实型变量。 如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) REAL IA(10),ND(80:99)使用隐含DO循环进行数组输入输出操作：例如 WRITE(*,10) ('I=',I,'A=',A(I),I=1,10,2) 10FORMAT(1X,5(A2,I2,1X,A2,I4)) 使用DATA语句给数组赋初值 变量表中可出现变量名,数组名,数组元素名,隐含DO循环，但不许出现任何形式的表达式：例如 DATA A,B,C/,, DATA A/,B/,C/ DATA A,B,C/3*CHARACTER*6 CHN(10) DATA CHN/10*''/INTEGER NUM(1000) DATA (NUM(I),I=1,500)/500*0/,(NUM(I),I=501,1000)/500*1/ 3 FORTRAN77程序书写规则
(完整版)Fortran之COMMON语句
Fortran之COMMON语句 FORTRAN程序中各程序单位之间的数据交换可以通过虚实结合来实现，还可以通过建立公用区的方式来完成。公用区有两种，一种是无名公用区，任何一个程序中只可能有一个无名公用区；一种是有名公用区，一个程序中可以根据需要由程序员开辟任意多个有名公用区。建立无名和有名公用区都通过COMMON语句来进行。 按照结构化程序设计的指导思想，通常主张通过虚实结合的途径来传送数据而不主张使用公用区。因为前者使程序有较好的可读性，而且容易跟踪数据的流向，便于对程序进行调试和维护。但是虚实结合的传送方式速度较慢，特别是在各程序单位之间有大量的数据需要传送时速度问题就更为突出。不同程序单位利用公用区交换数据的速度却比虚实结合的方式传送速度快得多。所以在很多应用程序中如果各程序单位之间有大量数据需要传送时，程序员一般仍然采用开辟公用区的方式。 11.2.1 无名公用区 开辟无名公用区COMMON语句的一般形式如下： COMMON a1,a2,… 其中a1,a2,…允许是普通变量名、数组名和数组说明符（注意：不是数组元素），它们之间用逗号隔开。例如： 在主程序中写：COMMON X，Y，I，Z（3） 在子程序中写：COMMON A，B，J，T（3） 于是，在无名公用区中变量X和A，Y和B，I和J分别分配在同一个存储单元中，数组Z和T同占三个存储单元。占同一个存储单元的那些变量在不同的程序单位中，它们的名字不需要相同。
FORTRAN编译程序在编译时为COMMON语句开辟一个无名公用区，当把不同的程序单位连接在一起的时候，不同程序单位在COMMON语句中的变量按其在语句中出现的先后顺序占用无名公用区连续的存储单元，也就是说每个程序单位中的COMMON语句总是给出现在语句中的第一个变量分配在无名公用区的第一个存储区单元。 COMMON语句开辟公用区的主要用途就是使不同程序单位的变量之间 进行数据传送，另一种用途是用来节省存储空间。下面的例子通过公用区进行交换数据： PROGRAM MAIN COMMON Z1,Z2 READ(*,*)A1,A2,A3 CALL QUAD(A1,A2,A3) WRITE(*,*)Z1,Z2 END SUBROUTINE QUAD(A,B,C) COMMON X1,X2 . . . END 或者: PROGRAM MAIN COMMON Z1,Z2,A1,A2,A3 READ(*,*)A1,A2,A3 CALL QUAD() WRITE(*,*)Z1,Z2 END
Simple算法_fortran语言编写
Cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc Cccc This computer program was copied from the graduate student course program Cccc of the University of Minnesota. Part of it was re-formulated to meet the Cccc personal computer environment. Some inappropriate expressions were also Cccc corrected. The program is used only for teaching purpose. No part of it Cccc may be published. You may use it as a frame to re-develop your own code Cccc for research purpose. XJTU Instructor, 1995.11 **************************************************************************** *----------------------------MAIN PROGRAM----------------------------------- **************************************************************************** LOGICAL LSTOP COMMON/CNTL/LSTOP **************************************************************************** OPEN(08,FILE='teresul') CALL SETUP0 CALL GRID CALL SETUP1 CALL START 10 CALL DENSE CALL BOUND CALL OUTPUT IF(.NOT.LSTOP) GO TO 15 CLOSE(08) STOP 15 CALL SETUP2 GO TO 10 END *--------------------------------------------------------------------------- SUBROUTINE DIFLOW **************************************************************************** COMMON/COEF/FLOW,DIFF,ACOF **************************************************************************** ACOF=DIFF IF(FLOW .EQ.0.0)RETURN TEMP=DIFF-ABS(FLOW)*0.1 ACOF=0. IF(TEMP .LE. 0. ) RETURN TEMP=TEMP/DIFF ACOF=DIFF*TEMP**5 RETURN END *-------------------------------------------------------------------------- SUBROUTINE SOLVE ****************************************************************************
fortran90-知识点整理
Fortran90关于变量说明的新功能： 1.在变量说明的同时，可以给变量赋初值。 INTEGER::I=5，J=126 REAL::X=7.2，Z，W=774.2 2.在说明变量的同时也可说明其种别 REAL(KIND=4)::X,Y 或REAL(4)::X,Y 3.在说明变量的同时，还可说明变量的属性 INTEGER，PARAMETER::I=5，J=123 Ps :PARAMETER 属性 （1）功能：用一个符号代表一个常量，称为符号常量（常数） （2）写法： Real，Parameter ：：G=9.8 &&说明类型时赋值 （3）位置：位于可执行语句之前 REAL，DIMENSION(1:10)::A 2.3.6 派生数据类型 根据需要而由基本数据类型定义新的数据类型。在一个派生类型中可包含多个基本类型。 如：TYPE STUDENT (定义开始) CHARACTER(LEN=20)::DEPARTMENT CHARACTER(LEN=10)::CLASS CHARACTER(LEN=15)::NAME INTEGER::NUMBER （成员定义） END TYPE STUDENT （定义结束）
派生类型变量的定义： TYPE（STUDENT）::PERSON 变量的赋值： PERSON=（”COMPUTER”，”92_2”，”LI LIN”，21） 成员的表示： PERSON%CLASS=“92_2” PERSON%NAME=“LI LIN” 主程序其他限制 主程序的可执行部分不能包含有RETURN语句或者ENTRY语句。 程序名对可执行程序是全局的，而且不得于该可执行程序中的任何其它程序单元名、外部过程名或公用块名相同，也不得于主程序内的任何局部名相同。 在主程序的作用范围内的说明不得包含OPTIONAL语句、INTENT语句、PUBLIC语句或它们的等价属性，在主程序内SAVE语句不起作用。 主程序内的任何内部过程的定义必须跟在CONTAINS语句之后。 2.3.2 子程序 子程序是可以完成某一独立算法或功能的程序单元，但它功能的具体体现要通过主程序（或子程序）的调用来实现。 按子程序与主程序的位置关系分为内部过程和外部过程。 2.4.3 模块 模块也是一种在程序单元之外独立编写的程序单元。它有独特的形式，即模块程序单元内没有可执行语句，除了说明语句外，最多包含内部过程。模块的主要作用是供其它程序单元引用。 程序单元引用模块，实际是将模块内的全部语句复制到本程序单元。因此模块起共享及复制的作用。 模块的引用采用语句： USE 模块1，模块2，…
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语言
#LOCAL REAL Uds,Uqs,angles,vs,w,o,ang,ta,tb,t0,Ud1,Uq1,E #LOCAL REAL tc,t,mt,pi #LOCAL INTEGER sign,M,nl,n t=TIME tc=$TC nl=DNINT(tc/delt) pi=3.1415926 E=$Udc w=2*pi*50 if(t.lt.delt) then sign=0 M=0 Uds=0 Uqs=0 angles=0 vs=0 o=0 ang=0 ta=0 tb=0 t0=0 Ud1=0 Uq1=0 n=1 end if if(mod(n,nl).lt.1E-7) then !初始化变量 M=M+1 Uds=$Ud Uqs=$Uq angles=$angle !变换到alphabeta坐标直角坐标
Ud1=cos(angles)*Uds+sin(angles)*Uqs Uq1=sin(angles)*Uds-cos(angles)*Uqs ! Ud1=$Ud ! Uq1=$Uq !生成矢量的模和角度 vs=sqrt(Ud1**2+Uq1**2) if(Ud1.gt.0.and.Uq1.ge.0) then ang=atan(Uq1/Ud1) else ang=atan(Uq1/Ud1)+2*pi end if if(Ud1.lt.0) then ang=atan(Uq1/Ud1)+pi end if if(Ud1.eq.0) then if(Uq1.ge.0) then ang=pi/2 else ang=3*pi/2 end if end if !判断矢量的所属,计算角度 if(ang.ge.0.and.ang.lt.pi/3) then sign=1 o=ang ! o=w*t end if if(ang.ge.pi/3.and.ang.lt.2*pi/3) then sign=2 o=ang-pi/3 ! o=w*t-pi/3 end if if(ang.ge.2*pi/3.and.ang.lt.3*pi/3) then sign=3 o=ang-2*pi/3 ! o=w*t-2*pi/3
Fortran90第六章、第七章习题答案
1.输入整数,判断能否被3或5整除,如能整除,则打印,否则不打印。编写程序实现之。解答： !判断可以被3或5整除的整数 PROGRAM exam61 INTEGER :: n PRINT*, '请输入一个整数:' READ*, n IF(mod(n,15)==0) THEN WRITE(*, "(1X,I5,'是一个可以被3和5整除的整数')") n ELSE IF(mod(n,3)==0) THEN WRITE(*, "(1X,I5,'是一个可以被3整除的整数')") n ELSE IF(mod(n,5)==0) THEN WRITE(*, "(1X,I5,'是一个可以被5整除的整数')") n ENDIF END 2.计算职工工资,工人每周工作40小时,超过40小时的部分应该按加班工资计算(为正常工资的2倍。输入工作时间和单位报酬,计算出该职工应得的工资,并打印输出。编写程序实现之。 解答： !计算职工工资 PROGRAM exam62 INTEGER t,p,pt PRINT*, '请输入工人工作的时间t和单位报酬pt' READ*, t,pt IF(t>40)then p=40*pt+(t-40)*2*pt ELSE p=t*pt ENDIF PRINT*, '该工人本周的应得工资为:',p,'元' END
1.已知:ｘ=0?、10?、20?、…、180?,输出x、sin(x)、cos(x)、tg(x)的值。编写程序实现之。 解答： ! 输出x、xin(x)、cos(x)、tan(x)的值 PROGRAM exam71 PARAMETER(pi=3.1415926) INTEGER i REAL :: x,sinx,cosx,tanx PRINT *,' x sin(x) cos(x) tg(x)' DO i=0,180,10 x=i/180.0*pi sinx=sin(x) cosx=cos(x) tanx=tan(x) PRINT '(I3,3(2X,F10.7))',i,sinx,cosx,tanx END DO END 2.已知:ｘ=1.0、1.1、1.2、…、2.9,输出x、x2、eｘ、ln(x)的值。编写程序实现之。解答： ! 计算x、x平方、e的x次方、ln(x)的值 PROGRAM exam72 INTEGER i REAL :: x,x2,ex,lnx PRINT*,' x x**2 exp(x) ln(x) ' DO i=10,29 x=i/10.0 x2=x*x expx=exp(x) lnx=log(x) PRINT '(F3.1,2X,F4.2,2X,F10.5,2X,F10.5)',x,x2,expx,lnx END DO END 3.输入10个整数,计算它们的和、积、平方和、和的平方。编写程序实现之。 解答： ! 计算和,积,平方和,和的平方 PROGRAM exam73 INTEGER i REAL :: s1=0.0,s2=1.0,s3=0.0,s4=0.0 PRINT*,'请输入十个实数(每行一个)' DO i=1,10 READ*,a s1=s1+a s2=s2*a s3=s3+a**2 END DO s4=s1**2 PRINT*,'10个数之和',s1 PRINT*,'10个数之积',s2
fortran初步学习资料
第一章: Fortran语言程序设计初步 Fortran语言的发展概况 本节介绍Fortran的起源与发展历史，讲述Fortran由产生到形成标准FortranIV、Fortran77，并进一步形成新标准Fortran90/95的发展历程。 1.1.1 Fortran的历史 a)a)FortranI FortranIV Fortran是目前国际上广泛流行的一种高级语言，适用于科学计算。Fortran 是英文FORmula TRANslatio n的缩写，意为“公式翻译”。它是为科学、工程问题中的那些能够用数学公式表达的问题而设计的语言，主要用于数值计算。这种语言简单易学，因为可以像抄写数学教科书里的公式一样书写数学公式，它比英文书写的自然语言更接近数学语言。Fortran语言是第一个真正推广的高级语言。至今它已有四十多年历史，但仍历久不衰，始终是数值计算领域所使用的主要语言。Fortran语言问世以来，根据需要几经发展，先后推出形成了很多版本。 第一代Fortran语言是在1954年提出来的，称为FortranI。它于1956年在IBM 704计算机上得以实现。在此之前编写计算机程序是极为繁琐的，程序员需要详细了解为之编写代码的计算机的指令、寄存器和中央处理器(CPU)等方面的知识。源程序本身是用数学符号(八进制码)编写的，后来采用了助记符，即所谓机器码或汇编码，这些编码由汇编程序转换为指令字。在50年代书写和调试一个程序要很长时间，因为用这种方式编写程序显然是很不方便的，尽管它能使CPU高效地工作。正是这些原因，促使由John Backus率领的IBM公司的一个小组研究开发最早的高级程序设计语言Fortran。其目的是开发一种容易理解、简单易学又能几乎像汇编一样高效运行的语言，他们取得了极大的成功。Fortran 语言作为第一种高级语言不仅是一次创新，也是一次革命。它使程序员摆脱了使用汇编语言的冗长乏味的负担，而且它使得不再只是计算机专家才能编写计算机程序，任何一名科学家或工程技术人员，只要稍加努力学习和使用Fortran，就能按自己的意图编写出用于科学计算的程序。 经过不断发展，FortranI形成了很多不同版本，其中最为流行的是1958年出现的FortranII，它对FortranI进行了很多扩充(如引进了子程序)，FortranII 在很多机器上得以实现。其后出现的FortranIII未在任何计算机上实现。1962年出现的FortranIV对原来的Fortran作了一些改变，使得FortranII源程序在FortranIV编译程序下不能全部直接使用，导致了语言不兼容的问题。这样就形成了当时同时使用FortranII和FortranIV两种程序设计语言的局面。
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