FORTRAN經(jīng)典入門程序20例(總8頁(yè)

1、對(duì)于FORTRAN的初學(xué)者。這些例子可作為小練習(xí)。1.例題：計(jì)算工資問(wèn)題。每小時(shí)工資為RATE,如果工作超過(guò)40小時(shí)，加班呢部分工資是正常時(shí)間工資的1.5倍。C Payroll with overtime program payroll real rate, hours, pay read (*,*) rate, hours if (hours40) then pay=40*rate+(hours-40)*1.5*rate else pay=hours*rate END IF print *,rate=, rate print *, hours=, hours print *,pay=,pay

2、 end2.學(xué)生成績(jī)問(wèn)題。大于80為A級(jí)。大于60小于80為B級(jí)。小于60為C級(jí)。 IF的嵌套。注意空格可以看清楚else if ,end if,print的內(nèi)容. 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 END3.三個(gè)數(shù)按從小到大排序。PROGRAM MAXMIN REAL A,B,C,T READ

3、*,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 C=T ELSE END IF IF (A.GT.B) THEN T=A A=B B=T END IF PRINT *,A,B,C END4.運(yùn)用EISE IF語(yǔ)句。重做例子2PROGRAM 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.

4、 60.0)THEN PRINT *, GRADE,=C ELSE PRINT *, GARDE,=D END IF END5. 計(jì)算 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 END6.CONTINUE語(yǔ)句。為了使循環(huán)的起止范圍清晰，使用終端語(yǔ)句例：求5！program ex1 FACT=1 do 10, N=1,5,1 #DO 語(yǔ)句標(biāo)號(hào)，循環(huán)不變量=初值，終值，步長(zhǎng) FACT=FACT*N10 continue print*,FA

5、CT end求1！+2！+3！+N！READ *,N FACT=1 sum=0.0 do 10, I=1,N,1 FACT=FACT*I SUM=SUM+FACT10 continue print*,SUM end7.求EX=1+X+X2/2!+X3/3!+.+XN/N! READ *,X,N sum=1 TERM=1 do 10, I=1,N,1 TERM=TERM*X/I SUM=SUM+TERM10 continue print*,EXP(,X,),SUM end8.打印九九乘法表 do 10,I=1,9,1 do 20,J=1,I PRINT*,I,*,J,=,I*J20 CONTIN

6、UE10 CONTINUE end9.百錢買百雞公雞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 IF30 CONTINUE20 CONTINUE10 CONTINUE END10.DO循環(huán)用來(lái)處理循環(huán)次數(shù)確定的情況比較方便。在處理無(wú)法事先知道循環(huán)次數(shù)時(shí)，使用當(dāng)型循環(huán)。(1)DO WHILE DO 標(biāo)號(hào)，WHILE(邏輯表達(dá)式)循環(huán)體CONT

7、INUE 標(biāo)號(hào)例 統(tǒng)計(jì)若干同學(xué)的得分并計(jì)算平均值。INTEGER SCORE,COUNT READ*,SCORE SUM=0 COUNT=0 DO 10, WHILE(SCORE .GT. 0.0) SUM=SUM+SCORE COUNT=COUNT+1 READ *,SCORE10 CONTINUE AVER=SUM/COUNT PRINT *,AVER END11.輾轉(zhuǎn)相除法求最大公約數(shù)。如26和6的最大公約數(shù)。先用27除以6，余數(shù)是3，再用6除以3，余數(shù)為0，則3就是最大公約數(shù)。命令MOD()求余數(shù)。INTEGER M,N,R,T READ*,M,NIF (M .LT. N)THEN T

8、=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 END12.判斷一個(gè)數(shù)是否是素?cái)?shù)。SQRT（）開(kāi)根號(hào) READ *,N MAX=sqrt(real(N) I=2 DO 10, WHILE (I .LE. MAX .AND. MOD(N,I) .NE. 0) I=I+110 CONTINUE IF(I .GT. MAX)THEN PRINT*,N,IS A PRIME NUMBER ELSE PRINT *,N,IS NO

9、T A PRIME NUMVER END IF END13.求100和200之間能被3或7整除的自然數(shù)。 I=10010 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 END14. 直到型循環(huán)，UNTIL 語(yǔ)句實(shí)現(xiàn)直到型循環(huán)。 格式“ DO S1,UNTIL(條件) 循環(huán)體 標(biāo)號(hào)S1 CONTINUE例14用直到型循環(huán)求5！ FORTRAN不支持UNTIL語(yǔ)句。INTEGER I,PRO

10、DUCT I=1 PRODUCT=1 DO 10, UNTIL(I .GT. 5) PRODUCT=PRODUCT*I I=I+110 CONTINUE PRINT *,5!=,PRODUCT ENDFORTRAN不支持UNTIL語(yǔ)句。讓我們改成DO WHILE格式。INTEGER I,PRODUCT I=1 PRODUCT=1 DO 10, WHILE(I .LT. 5) PRODUCT=PRODUCT*I I=I+110 CONTINUE PRINT *,5!=,PRODUCT END15. INTEGER R,I,J,S,T,SUMDO 10, I=100,999,1 T=I SUM=0

11、 DO 20 ,J=1,3,1 R=MOD(T,10) SUM=SUM+R*3 T=T/1020 CONTINUE IF (SUM .EQ. I) PRINT *,I,is a daffodil NUMBER10 CONTINUE END15.求一元二次方程的根 。DOUBLE PRECISION A, B, CCOMPLEX X1,X2READ *,A,B,CIF(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

12、/(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 END16. 打印此圖像。CHARACTER*9 LINE DO 10, I=-4,4,1 K=1 DO 40 J=1,ABS(I),1 LINE(K:K)= K=K+140 CONTINUE DO 20, J=1, 2*(4-ABS(I)+1,1 LINE(K:K)=* K=K+120 CONTINUE PRINT*,LINE DO 30, J=1,9,1 LINE(J:J)= 30 CONTINUE10

13、 CONTINUE END 17.輸出數(shù)組中的數(shù)值。 INTEGER A(1:10)DO 10,I=1,10,1 READ *,A(I)10 CONTINUE PRINT *,A END18. 利用DO循環(huán)對(duì)數(shù)組進(jìn)行輸入和輸出。從鍵盤上輸入十個(gè)數(shù)，按順序逐個(gè)放入到A(1)到A(10)中。 INTEGER A(1:10) DO 10,I=1,10,1 READ *,A(I)10 CONTINUE DO 20,I=1,9,2 PRINT *,A(,I,),A(I)20 CONTINUE END19. 雙重循環(huán)將從終端讀入的數(shù)據(jù)賦值給W數(shù)組中的各元素。REAL W(1:3,1:2) #表示一個(gè)三行兩

14、列的數(shù)組，當(dāng)前面數(shù)一樣時(shí)，可以寫出 W(3,2) DO 10,I=1,3,1 DO 20, J=1,2,1 READ *,W(I,J)20 CONTINUE10 CONTINUE PRINT*,W END20。輸入10名學(xué)生的成績(jī)，統(tǒng)計(jì)各個(gè)分?jǐn)?shù)段的人數(shù)。 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)+110 CONTINUE DO 20,I=0,10,1 PRINT *,COUNTER(I),=,I*10,(I+1)*10-120 CONTINUE END21.順序查找。一組數(shù)據(jù)從第一個(gè)開(kāi)始找，折半查找，順序查找適用數(shù)