FORTRAN经典入门程序20例

对于FORTRAN 的初学者。这些例子可作为小练习。

1.例题：计算工资问题。每小时工资为RATE, 如果工作超过40小时，加班呢部分工资是正常时间工资的1.5倍。

C Payroll with overtime program payroll real rate, hours, pay read (*,*) rate, hours if (hours>40) then

pay=40*rate+(hours-40)*1.5*rate

print *, "pay=" ,pay end

2.学生成绩问题。大于80为A级。大于60小于80为B级。小于60为C级。IF的嵌套。注意空格可以看清楚else if ,end if,print 的内容.

PROGRAM GRADESTUDENT

READ *,GRADE

IF (GRADE .GE. 80) THEN

PRINT *,GRADE, "GRADE=>A"

ELSE

IF (GRADE.LT.60) THEN

PRINT *,GRADE, "GRADE=>C"

ELSE

PRINT *,GRADE, "GRADE=>B"

END IF

END IF

END

3.三个数按从小到大排序。

PROGRAM MAXMIN

REAL A,B,C,T

READ *,A,B,C

IF (A.GT.B) THEN

END IF

IF (B.GT.C) THEN

T=B

B=C

PRINT *,A,B,C

END

4.运用EISE IF 语句。重做例子2

PROGRAM ex2

READ (*,*) GRADE

IF (GRADE .GE. 80.0) THEN

PRINT *, GRADE, "=>A"

ELSE IF(GRADE .GE. 70.0) THEN

PRINT *, GRADE, "=>B"

ELSE IF(GARDE .GE. 60.0) THEN

PRINT *, GRADE, "=>C"

ELSE

PRINT *, GARDE, "=>D"

END IF

END 5. 计算y 3x26,x 0

x2 2x 8,x 0

END

6.CONTINUE 语句。为了使循环的起止范围清晰，使用终端语句例：求5！program ex1

FACT=1

do 10, N=1,5,1 ##DO 语句标号，循环不变量=初值，终值，步长FACT=FACT*N 10 continue

print*,FACT

end

求1！+2！+3！+···+N！

READ *,N

FACT=1

sum=0.0

do 10, I=1,N,1

FACT=FACT*I

SUM=SUM+FACT

10 continue

print*,SUM

end

7.求E^X=1+X+X^2/2!+X^3/3!+...+X^N/N!

READ *,X,N

sum=1

TERM=1

do 10, I=1,N,1

TERM=TERM*X/I

SUM=SUM+TERM

10 continue

print*, "EXP(" ,X, ")" ,SUM end

8.打印九九乘法表

do 10,I=1,9,1

do 20,J=1,I 20 CONTINUE 10 CONTINUE

end

9.百钱买百鸡

公鸡5元每只。母鸡3元每只，小鸡1元三只。求100 元买一百只鸡

INTEGER X,Y,Z

DO 10,X=0,100

DO 20,Y=0,100

DO 30,Z=0,100

IF(X+Y+Z .EQ. 100) THEN

IF(15*X+9*Y+Z .EQ. 300) PRINT *,X,Y,Z

END IF

30 CONTINUE

20 CONTINUE

10 CONTINUE

END

10.DO 循环用来处理循环次数确定的情况比较方便。在处理无法事先知道循环次数时，使

用当型循环。

(1)DO WHILE DO 标号，WH ILE（逻辑表达式）

CONTINUE 标号

例统计若干同学的得分并计算平均值。INTEGER SCORE,COUNT

READ *,SCORE

SUM=0

COUNT=0

DO 10, WHILE (SCORE .GT. 0.0)

PRINT *,I, "*"

SUM=SUM+SCORE

COUNT=COUNT+1

READ *,SCORE

10 CONTINUE

AVER=SUM/COUNT

PRINT *,AVER

END

11.辗转相除法求最大公约数。如26 和6 的最大公约数。先用27 除以6，余数是3，再用6 除以3，余数为0，则3 就是最大公约数。命令MOD() 求余数。

INTEGER M,N,R,T

READ *,M,N

IF (M .LT. N) THEN

T=M

M=N

N=T

END IF

R=MOD (M,N)

DO 10, WHILE (R .NE. 0)

M=N

N=R

R=MOD (M,N)

10 CONTINUE

PRINT *, "zui da gong yue shu",N

END

12.判断一个数是否是素数。SQRT()开根号

READ *,N

MAX= sqrt(real(N))

I=2

DO 10, WHILE (I .LE. MAX .AND. MOD (N,I) .NE. 0)

I=I+1

10 CONTINUE

IF(I .GT. MAX) THEN

PRINT *,N, "IS A PRIME NUMBER"

ELSE

PRINT *,N, "IS NOT A PRIME NUMVER"

END IF

END

13.求100和200之间能被3或7整除的自然数。

I=100 10 IF(I .LE. 200) THEN

IF (MOD (I,3) .EQ. 0 .OR. MOD (I,7) .EQ. 0) THEN

PRINT *,I, "ABLE"

ELSE

PRINT *,I, "UNABLE" END IF

I=I+1

GOTO 10

END IF

END

14.直到型循环，UNTIL 语句实现直到型循环。

格式“ DO S1,UNTIL( 条件)

循环体

标号S1 CONTINUE

例14 用直到型循环求5 ！FORTRAN 不支持UNTIL 语句。

INTEGER I,PRODUCT

I=1

PRODUCT=1

DO 10, UNTIL (I .GT. 5)

PRODUCT=PRODUCT*I

I=I+1

10 CONTINUE

PRINT *, "5!=" ,PRODUCT

END

FORTRAN 不支持UNTIL 语句。让我们改成DO WHILE 格式。INTEGER I,PRODUCT I=1

PRODUCT=1

DO 10, WHILE (I .LT. 5)

PRODUCT=PRODUCT*I

I=I+1

10 CONTINUE

PRINT *, "5!=" ,PRODUCT

END

15.

INTEGER R,I,J,S,T,SUM

DO 10, I=100,999,1

T=I

SUM=0

DO 20 ,J=1,3,1

R=MOD (T,10)

SUM=SUM+R**3

T=T/10 20 CONTINUE

IF (SUM .EQ. I) PRINT *,I, "is a daffodil NUMBER"

10 CONTINUE

END

Ax2Bx C 0

DOUBLE PRECISION A, B, C

COMPLEX X1,X2 READ *,A,B,C IF(B**2-4*A*C .GE.0) THEN

X1=CMPLX ((-B+SQRT(B**2-4*A*C))/(2*A),0.0)

X2=CMPLX ((-B- SQRT(B**2-4*A*C))/(2*A),0.0) ELSE

X1=CMPLX ((-B/(2*A)), SQRT(-(B**2-4*A*C))/(2*A))

X2=CMPLX ((-B/(2*A)),- SQRT(-(B**2-4*A*C))/(2*A)) END IF PRINT *,X1,X2

END

16. 打印此图

CHARACTER *9 LINE

DO 10, I=-4,4,1

K=1

DO 40 J=1,ABS(I),1

LINE(K:K)= " "

K=K+1

40 CONTINUE

DO 20, J=1, 2*(4-ABS(I))+1,1

LINE(K:K)= "*" K=K+1

20 CONTINUE

PRINT *,LINE

DO 30, J=1,9,1 LINE(J:J)= " " 30 CONTINUE 10 CONTINUE END

17. 输出数组中的数值。

INTEGER A(1:10) DO 10,I=1,10,1 READ *,A(I)

10 CONTINUE PRINT *,A END

18.

利用 DO 循环对数组进行输入和输出。从键盘上输入十个数，按顺序逐个放入到

19.

双重循环将从终端读入的数据赋值给 W 数组中的各元素。

REAL W(1:3,1:2) ##表示一个三行两列的数组，当前面数一样时，可以写出 W(3,2) DO 10,I=1,3,1 DO 20, J=1,2,1

20。输入 10名学生的成绩，统计各个分数段的人数。

INTEGER GRADE(1:10),COUNTER(0:10)

READ *,(GRADE(I),I=1,10,1)

DO 10, I=1,10,1 COUNTER(GRADE(I)/10)=COUNTER(GRADE(I)/10)+1

10 CONTINUE

DO 20,I=0,10,1

PRINT *,COUNTER(I), "=>" ,I*10,(I+1)*10-1

20 CONTINUE

END

21. 顺序查找。一组数据从第一个开始找，折半查找，顺序查找适用数据少的情况，折半查找用于数 据多的情况，但数据必须按顺序排列。

INTEGER A(1:7)

A(1)到A(10) 中

DATA A/55,62,10,58,22,33,66/ READ *,X

DO 10, I=1,7,1