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

else

pay=hours*rate

END IF

print *,"rate=", rate

print *, "hours=", hours

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

T=A

A=B

B=T

ELSE

END IF

IF (B.GT.C) THEN

T=B

B=C

PROGRAM EQUATION

READ (*,*) X

IF (X .GE. 0.0) Y=3*X+6

IF (X .LT. 0.0) Y=-X**2+2*X-8

PRINT *,"X=",X, "Y=",Y

END

6.CONTINUE语句。为了使循环的起止范围清晰，使用终端语句

例：求5！

program ex1

FACT=1

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

10continue

print*,FACT

end

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

READ *,N

FACT=1

do 10, I=1,N,1

FACT=FACT*I

SUM=SUM+FACT

10continue

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

10continue

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

end

8.打印九九乘法表

do 10,I=1,9,1

do 20,J=1,I

PRINT*,I,"*",J,"=",I*J

20CONTINUE

10CONTINUE

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

30CONTINUE

20CONTINUE

10CONTINUE

END

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

(1)DO WHILE

DO 标号，WHILE(逻辑表达式)

循环体

CONTINUE 标号

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

INTEGER SCORE,COUNT

READ*,SCORE

SUM=0

COUNT=0

DO 10, WHILE(SCORE .GT. 0.0)

SUM=SUM+SCORE

COUNT=COUNT+1

READ *,SCORE

10CONTINUE

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)

10CONTINUE

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

10CONTINUE

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

10IF(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

10CONTINUE

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

10CONTINUE

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

15.求一元二次方程的根20Ax Bx C ++= 。

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

30CONTINUE

10CONTINUE

END

17.输出数组中的数值。

INTEGER A(1:10)

DO 10,I=1,10,1

READ *,A(I)

10CONTINUE

PRINT *,A

END

18.利用DO循环对数组进行输入和输出。从键盘上输入十个数，按顺序逐个放入到A(1)到A(10)中。INTEGER A(1:10)

DO 10,I=1,10,1

READ *,A(I)

10CONTINUE

DO 20,I=1,9,2

PRINT *,"A(",I,")",A(I)

20CONTINUE

END

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

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

DO 20, J=1,2,1

READ *,W(I,J)

20CONTINUE

10CONTINUE

PRINT*,W

END

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

10CONTINUE

DO 20,I=0,10,1

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

20CONTINUE

END

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

INTEGER A(1:7)

DATA A/55,62,10,58,22,33,66/

READ *,X

DO 10, I=1,7,1

IF(A(I) .EQ. X) GOTO 20

10CONTINUE

20IF (I>7)THEN

PRINT *,X,"it is not here"

ELSE

PRINT *,X,"it is in the A the",I,"position"

END IF

END