Fortran90幾個重要的實驗報告

1、Fortran90語言程序設(shè)計 上機實驗報告河南科技大學(xué)規(guī)劃與建筑工程學(xué)院班級： 工力112 姓名： 學(xué)號： 111409070210 成績： 2013年 5月 19 日實驗一（1） 熟悉FORTRAN90軟件開發(fā)環(huán)境一 實驗?zāi)康?#160;1 了解FORTRAN90與軟件開發(fā)環(huán)境的關(guān)系。2 掌握FORTRAN90上機實驗基本操作過程。3 掌握FORTRAN90軟件開發(fā)環(huán)境啟動方法。4 理解有關(guān)工作區(qū)、項目、文件的基本概念。5 掌握工作區(qū)、項目、文件創(chuàng)建方法。6 熟悉Compaq Visual Fortran軟件開發(fā)環(huán)境圖形用戶界面。7 掌握程序的輸入、編輯、編譯、構(gòu)建、運行方法。8 了解軟件

2、開發(fā)環(huán)境常用菜單、工具按鈕、環(huán)境窗口基本功能。9 掌握文件、文件夾的移動、復(fù)制、刪除等操作方法。 二 課堂實驗內(nèi)容 1 問題描述 現(xiàn)有長1500米，寬1000米的地塊要拍賣建造住宅小區(qū)，地塊中央有一半徑為100米的圓形公共綠地不屬于拍賣范圍，每平方米地價為5000元。編寫程序計算該地塊拍賣總價，并輸出。要求地塊長寬和綠地半徑數(shù)據(jù)從鍵盤輸入。 2 算法設(shè)計通過對該問題的分析研究，設(shè)計求解算法，并繪制流程圖如圖1-1所示輸入月銷售罐數(shù)：number輸入月銷售罐數(shù)：number計算油罐容積：cubage計算月銷售數(shù)量：volumes計算月銷售額：sale計算月

3、銷售額：sale輸出月銷售額：sale輸出月銷售額：sale輸出月銷售額：sale輸出月銷售額：sale開始結(jié)束圖1-1 計算地塊總價程序流程圖3. 程序編寫根據(jù)算法和流程圖，設(shè)計和編寫程序如下：Ø Ø  !班級：工力112Ø Ø  !姓名：Ø Ø  !時間：201304Ø Ø  !計算加油站月銷售額Ø Ø  PROGRAM mainØ Ø  PARAMETER (pi=3.1415)Ø Ø

4、60; INTEGER : numberØ Ø  REAL ： cubage,volumes,saleØ Ø  REAL ： radius=5.0,heigh=3.0,unitprice=500.0Ø Ø  PRINT*,請輸入月銷售油罐數(shù)量：Ø Ø  READ*，numberØ Ø  cubage=pi*radius*radius*heighØ Ø  volumes=cubage*numberØ 

5、6;  sale=volumes*unitpriceØ Ø  PRINT*, 加油站月銷售額：，sale，元。Ø Ø  END4. 運行結(jié)果 請輸入月銷售油罐數(shù)量:4加油站月銷售額： 471225.0 元。Press any key to continue 三課外練習(xí)1 問題描述. 1有一半徑為5，高為20的圓柱型物體，求圓柱體積，并輸出。編寫程序?qū)崿F(xiàn)之。 2有一半徑為5，高為20的圓錐型物體，求圓錐體積，并輸出。編寫程序?qū)崿F(xiàn)之。  2 算法設(shè)計輸入半徑和高：r,h計算圓柱體積：V輸出圓柱體積：V開始結(jié)束

6、輸入半徑和高：r,h計算圓錐體積：V輸出圓錐體積：V開始結(jié)束3. 程序編寫程序2：Ø Ø  !班級：工力112Ø Ø  !姓名：Ø Ø  !時間：2101304Ø Ø  !計算圓錐體積Ø Ø  PROGRAM mainØ Ø  PARAMETER (pi=3.1415)Ø Ø  REAL ： r,h,VØ Ø  PRINT*,請輸入圓錐底半徑和高：Ø

7、; Ø  READ*，r,hØ Ø  V=pi*r*r*h/3Ø Ø  PRINT*, 圓錐體積：，VØ Ø  END 根據(jù)流程圖，分別編寫求解上述兩問題的程序：程序1、程序2。程序1：Ø Ø  !班級：工力112Ø Ø  !姓名：Ø Ø  !時間：201304Ø Ø  !計算圓柱體積Ø Ø  PROGRAM mainØ

8、; Ø  PARAMETER (pi=3.1415)Ø Ø  REAL ： r,h,VØ Ø  PRINT*,請輸入圓柱底半徑和高：Ø Ø  READ*，r,hØ Ø  V=pi*r*r*hØ Ø  PRINT*, 圓柱體積：，VØ Ø  END4. 運行結(jié)果請輸入圓柱底半徑和高：5，20 圓柱體積： 1570.750Press any key to continue請輸入圓錐底半徑和高：5，20

9、 圓錐體積： 523.5833Press any key to continue實驗一（2） 掌握FORTRAN90基礎(chǔ)知識一 實驗?zāi)康?#160; 1 了解FORTRAN90字符集。2 掌握FORTRAN90實體名稱的語法描述規(guī)則和具體使用要求。3 了解FORTRAN90關(guān)鍵字意義，以及與名稱的區(qū)別。4 掌握FORTRAN90程序組織結(jié)構(gòu)和語句排列次序的基本原則。5 掌握FORTRAN90自由書寫格式及三個標(biāo)志符的主要作用。6 理解數(shù)據(jù)、數(shù)據(jù)表示、數(shù)據(jù)類型概念，熟悉FORTRAN90數(shù)據(jù)類型。7 進(jìn)一步熟 二 課堂實驗內(nèi)容 1 問題描述   

10、有一呈等腰梯形的魚塘，上底邊長為150米，下底邊長為210米，高為40米，在魚塘周圍建起了圍墻，如圖2-1所示。每平方米魚塘產(chǎn)魚4公斤，每公斤魚的價格為10元。編寫程序計算該魚塘的產(chǎn)量和產(chǎn)值，以及圍墻的長度。要求梯形的上底邊長、下底邊長、高從鍵盤輸入。  a:150b:210h:40c       圖2-1 魚塘示意圖 2 算法設(shè)計開始輸入梯形的上底和下底以及高：a,b,h計算梯形面積：area計算產(chǎn)量和產(chǎn)值：yield，sales計算梯形斜高：c計算魚塘的圍墻長度：perimeter結(jié)束3. 程序編寫Ø

11、 Ø  !班級：工力112Ø Ø  !姓名：Ø Ø  !時間：201304Ø Ø  !計算魚塘產(chǎn)量和產(chǎn)值Ø Ø  PROGRAM mainØ Ø  INTEGER ：p=4,q=10Ø Ø  REAL ： a,b,c,hØ Ø  REAL ： area,perimeter，yield，salesØ Ø  PRINT*，請輸入梯形魚塘的上底

12、、下底和高值：Ø Ø  READ*，a,b,hØ Ø  area=(a+b)*h/2Ø Ø  yield=area*pØ Ø  sales=yield*qØ Ø  c=SQRT(b-a)/2)*(b-a)/2)+h*h)Ø Ø  perimeter=a+b+2*cØ Ø  PRINT*， 魚塘產(chǎn)量：，yield，公斤。Ø Ø  PRINT*， 魚塘產(chǎn)值：，s

13、ales，元。Ø Ø  PRINT*， 圍墻長度：，perimeter，米。Ø Ø  ENDØ4. 運行結(jié)果請輸入梯形魚塘的上底、下底和高值：150,210,40 魚塘產(chǎn)量： 28800.00 公斤。 魚塘產(chǎn)值： 288000.0 元。 圍墻長度： 460.0000 米。Press any key to continue三. 課外練習(xí)1 問題描述：有一長5米、寬3米的長方形鐵皮，在四個角對稱挖了四個半徑為0.2米的圓洞， 根據(jù)需要對這塊鐵皮鍍一層銅，鍍銅的價格是每平方米50元。計算鍍銅費用。   2 算法

14、設(shè)計開始輸入矩形邊長及圓洞半徑：length、width、radius計算鐵皮的面積：area計算鍍銅費用：rate輸出3. 程序編寫程序3-1:Ø Ø  !班級：工力112Ø Ø  !姓名：Ø Ø  !時間：201304Ø Ø  !計算鍍銅費用Ø Ø  PROGRAM mainØ Ø  INTEGER ： p=50Ø Ø  REAL ： length,width,radiusØ

15、; Ø  REAL ： area,rateØ Ø  PRINT*,請輸入鐵皮長寬和圓洞半徑：Ø Ø  READ*，length,width,radiusØ Ø  area=length*width-4*3.1415*radius*2Ø Ø  rate=p*areaØ Ø  PRINT*, 鍍銅費用：，rateØ Ø  END4. 運行結(jié)果請輸入鐵皮長寬和圓洞半徑：5,3,0.2 鍍銅費用: 724

16、.8680Press any key to continue實驗二（1） 內(nèi)部數(shù)據(jù)類型與簡單順序程序設(shè)計1 實驗?zāi)康? 了解內(nèi)部數(shù)據(jù)類型及其KIND值參數(shù)概念。2 掌握常量、變量、表達(dá)式、賦值語句的語法規(guī)則和基本功能。3 掌握常用標(biāo)準(zhǔn)函數(shù)的基本功能和使用要求。4 掌握表控輸入輸出的基本功能和使用要求。5 掌握順序結(jié)構(gòu)概念和順序結(jié)構(gòu)程序設(shè)計方法。進(jìn)一步熟悉軟件開發(fā)環(huán)境和掌握程序的輸入、編輯、編譯、構(gòu)建和運行方法。  二 課堂實驗內(nèi)容 1 問題描述已知兩復(fù)數(shù)：3.5+4.8i、7.5-3.2i。計算其和、差、積、商，并輸出計算結(jié)果。編寫程序?qū)崿F(xiàn)之。要求兩復(fù)數(shù)從鍵盤輸入

17、。  2 算法設(shè)計開始輸入兩個復(fù)數(shù)：c1,c2計算兩數(shù)之和并賦值于a計算兩復(fù)數(shù)之差并賦值于b計算復(fù)數(shù)之積賦值給c計算復(fù)數(shù)之商賦值給d輸出a,b,c,d結(jié)束3. 程序編寫根據(jù)算法和程序流程圖，設(shè)計和編寫程序如下。Ø !班級：工力112Ø !姓名：Ø !時間：201304Ø !計算兩復(fù)數(shù)的和、差、積、商Ø PROGRAM mainØ COMPLEX ： c1,c2Ø complex： a,b,c,dØ PRINT*，請輸入兩個復(fù)數(shù)：Ø READ*，c1,c2Ø a=c1+c2；b

18、=c1-c2Ø c=c1*c2；d=c1/c2Ø PRINT*， 復(fù)數(shù)之和：，aØ PRINT*， 復(fù)數(shù)之差：，bØ PRINT*， 復(fù)數(shù)之積：，cØ PRINT*， 復(fù)數(shù)之商：，dØ ENDØ !復(fù)數(shù)變量占有8字節(jié)4. 運行結(jié)果 請輸入兩個復(fù)數(shù)：3.5004.8007.500-3.2 復(fù)數(shù)之和： (11.00000,1.600000) 復(fù)數(shù)之差： (-4.000000,8.000000) 復(fù)數(shù)之積： (41.61000,24.80000) 復(fù)數(shù)之商： (0.1637840,0.7098812)Press any key t

19、o continue三. 課外練習(xí)1 問題描述某炮兵陣地，某次使用大炮向敵人陣地猛烈發(fā)射炮彈。炮彈的運行軌跡滿足方程：f(x)=-0.0000125x2+0.07x-50,計算炮彈能打擊的最大高度和最遠(yuǎn)距離。編寫程序?qū)崿F(xiàn)之。問題分析：（提示）YXOx2x1設(shè)：y=f(x)=ax2+bx+c，軌跡線為一拋物線，如圖3-5所示。ymax圖3-5 一元二次方程軌跡曲線一元二次方程求根公式為：假定方程只有兩個不同實數(shù)根。  2 算法設(shè)計開始請輸入a,b,c計算兩根x1,x2計算射程length和最大高度輸出射程和最大高度：length、h結(jié)束!班級:工程力學(xué)112!姓名:!時間:2

20、0130426晚上!計算計算大炮的最大打擊高度和最遠(yuǎn)距離，并輸出PROGRAM mainREAL:a,b,c ,d,x1,x2,xREAL:length ,hPRINT*,'請輸入方程參數(shù)a,b,c'READ*,a,b,cd=b*b-4*a*cx1=(-b+SQRT(d)/(2*a)x2=(-b-sqrt(d)/(2*a)x=0.5*(x1+x2)h=a*x*x+b*x+clength=x2-x1PRINT*,'最大高度是:','h=',hPRINT*,'最遠(yuǎn)射程是:','length=',lengthEND3.

21、程序編寫4. 運行結(jié)果 請輸入方程參數(shù)a,b,c-0.0000127,0.07,-50 最大高度是:h= 46.45669 最遠(yuǎn)射程是:length= 3825.183Press any key to continue實驗二（2） 數(shù)據(jù)有格式輸入輸出一 實驗?zāi)康?  了解數(shù)據(jù)有格式輸入輸出的主要作用、基本方法和使用規(guī)則。2 掌握數(shù)據(jù)有格式輸入輸出語句(READ、PRINT、WRITE)的基本功能和語法規(guī)則。3 掌握格式說明語句(FORMAT)的基本功能和語法規(guī)則。4 掌握I、F、E、A、L、X、/、等若干常用格式編輯符的格式要求和基本功能。5 掌握數(shù)據(jù)有格式輸入輸出的數(shù)據(jù)組織和格式說

22、明方法。6 進(jìn)一步掌握順序結(jié)構(gòu)程序設(shè)計方法。7 進(jìn)一步熟悉軟件開發(fā)環(huán)境。 二 課堂實驗內(nèi)容 1 問題描述 通常需要計算(a+b)n的展開式系數(shù)，展開式系數(shù)滿足楊輝三角關(guān)系。計算并輸出n=1，2，3，4，5時的展開式系數(shù)及系數(shù)之和，即下面的楊輝三角。編寫程序?qū)崿F(xiàn)之。字母“a”、字母“b”、二項式階數(shù)n從鍵盤輸入。 1 1n=1 1 1 2n=2 1 2 1 4n=3 1 3 3 1 8n=4 1 4 6 4 1 16n=5 1 5 10 10 5 1 32要求輸入輸出數(shù)據(jù)以規(guī)范、整潔、直觀形式給出。程序中應(yīng)采用有格式輸入輸出。二項式階數(shù)n為2位整數(shù)，項符號為長度等于

23、1的字符，系數(shù)為3位整數(shù)，系數(shù)和4位整數(shù)。 2 算法設(shè)計開始根據(jù)提示輸入二項式第一、二項字母和階數(shù)：x、y、n計算n=1時二項式系數(shù)c11,c12及系數(shù)和s1計算n=2時二項式系數(shù)c21,c22,c23及系數(shù)和s2計算n=3時二項式系數(shù)c31,c32,c33,c34及系數(shù)和s3計算n=4時二項式系數(shù)c41,c42,c43,c44,c45及系數(shù)和s4。按輸出格式要求輸出數(shù)據(jù)。計算n=5時二項式系數(shù)c51,c52,c53,c54,c55,c56及系數(shù)和s5。結(jié)束3. 程序編寫!班級：工力112!姓名：!時間：20105!計算并輸出二項式系數(shù)（楊輝三角）PROGRAM yanghui_tr

24、iangleINTEGER(1) : nINTEGER: x,yINTEGER : c11=1,c12=1,s1INTEGER : c21=1,c22,c23=1,s2INTEGER : c31=1,c32,c33,c34=1,s3INTEGER: c41=1,c42,c43,c44,c45=1,s4INTEGER : c51=1,c52,c53,c54,c55,c56=1,s5WRITE(*,"(1X,'請輸入二項式第一項字母(1個字母):',)")READ(*,"(A)") xWRITE(*,"(1X,'請輸入二項式

25、第二項字母(2個字母):',)")READ(*,"(A)") yWRITE(*,"(1X,'請輸入二項式階數(shù)(2位整數(shù)):',)")READ(*,"(I2)") ns1=c11+c12c22=c11+c12; s2=c21+c22+c23c32=c21+c22; c33=c22+c23; s3=c31+c32+c33+c34c42=c31+c32;c43=c32+c33;c44=c33+c34; s4=c41+c42+c43+c44+c45c52=c41+c42; c53=c42+c43; c54=c

26、43+c44; c55=c44+c45; s5=c51+c52+c53+c54+c55+c56WRITE(*,"(1X,'下面給出(',A1,'+',A1,')1到',I2, '次方的展開式系數(shù)及系數(shù)之和。')") x,y,nWRITE(*,"(1X,'階數(shù)',9X,'楊輝三角',7X,'系數(shù)之和')") WRITE(*,"(1X,15X,I3)") 1WRITE(*,"(1X,'n=',I2,8X

27、,2(1X,I3),9X,I4)") 1,c11,c12,s1WRITE(*,"(1X,'n=',I2,6X,3(1X,I3),7X,I4)") 2,c21,c22,c23,s2WRITE(*,"(1X,'n=',I2,4X,4(1X,I3),5X,I4)") 3,c31,c32,c33,c34,s3WRITE(*,"(1X,'n=',I2,2X,5(1X,I3),3X,I4)") 4,c41,c42,c43,c44,c45,s4WRITE(*,"(1X,'n

28、=',I2, 6(1X,I3),1X,I4)") 5,c51,c52,c53,c54,c55,c56,s5END 請輸入二項式第一項字母(1個字母):1 請輸入二項式第二項字母(2個字母):1,1 請輸入二項式階數(shù)(2位整數(shù)):5 下面給出(1+1)1到 5次方的展開式系數(shù)及系數(shù)之和。 階數(shù) 楊輝三角 系數(shù)之和 1 n= 1 1 1 2 n= 2 1 2 1 4 n= 3 1 3 3 1 8 n= 4 1 4 6 4 1 16 n= 5 1 5 10 10 5 1 32Press any key to continue4. 運行結(jié)果三. 課外練習(xí)1 問題描述使用數(shù)據(jù)有格式輸入

29、輸出功能，輸出由4個字符abcd組成的圖形。編寫程序?qū)崿F(xiàn)之。 a b b c c ddddddddddddd c c b b a問題分析：（提示）a、b、c、d四個字符可以是任意可顯示字符，這四個字符從鍵盤輸入得到。若輸入的四個字符為*、#、&、$，則輸出結(jié)果為： * # # & & $ & & # # *由相同字符組成的字符串,可用標(biāo)準(zhǔn)函數(shù)REPEAT(str,I)得到。如REPEAT($,3)= $。  2 算法設(shè)計開始結(jié)束3. 程序編寫4. 運行結(jié)果實驗三（1） 選擇結(jié)構(gòu)序程序設(shè)計一 實驗?zāi)康?  熟悉關(guān)系表達(dá)式和邏輯

30、表達(dá)式。2 掌握邏輯IF、塊IF、多支塊IF、SELECT語句的語法規(guī)則和使用要求。3 掌握嵌套選擇結(jié)構(gòu)和嵌套選擇語句的使用要求和實現(xiàn)方法。4 掌握選擇結(jié)構(gòu)概念和選擇結(jié)構(gòu)程序設(shè)計方法。 二 課堂實驗內(nèi)容 1 問題描述計算下面分段函數(shù)的值，編寫程序?qū)崿F(xiàn)之。x<0x=0x>0  2 算法設(shè)計圖5-2 程序框結(jié)束輸出y的值用公式三計算y開始x < 0用公式一計算y假真x = 0假真用公式二計算y輸入x的值5 程序編寫Ø !班級：工力112Ø !姓名：Ø !時間：201305Ø !計算分段函數(shù)的值

31、6; PROGRAM exampleØ REAL(4) ： x,yØ PRINT*,'請輸入x的值：'Ø READ*，xØ IF (x<0) THENØ y=exp(2*abs(x)+sin(x*3.14159/180)Ø ELSEIF ( x=0) THENØ y=2Ø ELSEØ y=x/sqrt(1+x*x)Ø END IFØ PRINT *, 'y=', yØ END 請輸入x的值：-5.8 y= 109097.7Press an

32、y key to continue 請輸入x的值：0 y= 2.000000Press any key to continue 請輸入x的值：5.8 y= 0.9854601Press any key to continue4. 運行結(jié)果三. 課外練習(xí)1 問題描述  某商店對購貨額為1000元以上(含1000元,下同),八折優(yōu)惠;500元以上1000元以下,九折優(yōu)惠;200元以上500元以下,九五折優(yōu)惠;100元以上200元以下,九七折優(yōu)惠;100元以下,不優(yōu)惠。輸入購貨額,計算優(yōu)惠貨款,并輸出。2 算法設(shè)計開始輸入購物額x假真x>1000假真x>500y=0.

33、8*xx>200真y=0.97*x假x>0真y=xy=0.95*x輸出y!班級:工程力學(xué)112!姓名：!時間：201305!程序用途：計算超市購物優(yōu)惠金額PROGRAM mainreal(4): x,yPRINT*,'請輸入購貨金額'read*,xIF (x>=1000) THEN y=0.8*xELSE IF(X>=500) THEN y=0.9*x ELSE IF(X>=200) then y=0.95*x else IF(X>=100) THEN y=0.97*x else y=xEND IFPRINT*,'優(yōu)惠金額為：'

34、;,yEND結(jié)束3. 程序編寫請輸入購貨金額1200 優(yōu)惠金額為： 960.0000Press any key to continue 請輸入購貨金額600 優(yōu)惠金額為： 540.0000Press any key to continue 請輸入購貨金額250 優(yōu)惠金額為： 237.5000Press any key to continue 請輸入購貨金額150 優(yōu)惠金額為： 145.5000Press any key to continue 請輸入購貨金額50 優(yōu)惠金額為： 50.00000Press any key to continue4. 運行結(jié)果請輸入購貨金額1200 優(yōu)惠金額為：

35、960.0000Press any key to continue 請輸入購貨金額600 優(yōu)惠金額為： 540.0000Press any key to continue 請輸入購貨金額250 優(yōu)惠金額為： 237.5000Press any key to continue 請輸入購貨金額150 優(yōu)惠金額為： 145.5000Press any key to continue 請輸入購貨金額50 優(yōu)惠金額為： 50.00000Press any key to continue實驗三(2) 循環(huán)結(jié)構(gòu)序程序設(shè)計一 實驗?zāi)康?#160;1 掌握關(guān)系表達(dá)式和邏輯表達(dá)式。2 掌握循環(huán)控制條件的描述。3

36、掌握三種循環(huán)結(jié)構(gòu)的特征及相互之間的等價轉(zhuǎn)換。4 掌握DO、DO WHILE語句的語法規(guī)則和使用要求。5 掌握循環(huán)結(jié)構(gòu)概念和循環(huán)結(jié)構(gòu)程序設(shè)計方法。二 課堂實驗內(nèi)容 1 問題描述 驗證哥德巴赫猜想。哥德巴赫提出，一個不小于6的偶數(shù)必定能表示為兩個素數(shù)之和。如：8=3+5。求6至20之間的所有偶數(shù)的兩個素數(shù)之和。編寫程序?qū)崿F(xiàn)之。輸出格式為：6=3+38=3+510=3+710=5+5  開始2 算法設(shè)計請輸入偶數(shù)n從3開始驗證循環(huán)得出a,b輸出n,a,b結(jié)束3. 程序編寫 程序6-2：Ø !班級：工力112Ø !姓名：Ø !時間：20130

37、5Ø !驗證哥德巴赫猜想Ø PROGRAM mainØ INTEGER : a,b,i,m,nØ out: DO n=6,20,2Ø mid: DO a=3,n/2,2Ø m=INT(SQRT(REAL(a)Ø i=3Ø in1: DO WHILE ( (i<=m) .AND. (MOD(a,i)/=0)Ø i=i+2Ø ENDDO in1Ø IF (i>m) THENØ b=n-aØ m=0 i=3Ø in2: DO WHILE (i<

38、=m).AND.(MOD(b,i)/=0)Ø i=i+2Ø ENDDO in2Ø IF ( b>=0) THEN 6 = 3 + 3 8 = 3 + 5 10 = 3 + 7 10 = 5 + 5 12 = 3 + 9 12 = 5 + 7 14 = 3 + 11 14 = 5 + 9 14 = 7 + 7 16 = 3 + 13 16 = 5 + 11 16 = 7 + 9 18 = 3 + 15 18 = 5 + 13 18 = 7 + 11 20 = 3 + 17 20 = 5 + 15 20 = 7 + 13Press any key to cont

39、inueØ PRINT *,n,'=',a,'+',bØ ENDIFØ ENDIFØ END DO midØ END DO outØ END4. 運行結(jié)果三. 課外練習(xí)1 問題描述 輸入x值,按下列公式計算cos(x)。精度要求7位有效數(shù)字,最后一項<10-7。編寫程序?qū)崿F(xiàn)之。2 算法設(shè)計請輸入x開始n=n+1真T>eps假輸出cos x!班級:工力112!姓名：!時間：20130512!求cos xPROGRAM cos_xPARAMETER (PI=3.1415927,EPS=

40、1e-7) INTEGER :n=1,max=10REAL x,T,COSXWRITE(*,"(A)",ADVANCE='NO')'請輸入角度值（實數(shù)）:'READ*,xX=X*PI/180 !角度轉(zhuǎn)換公式T=1;COSX=TDO WHILE(ABS(T)>EPS).AND.(n<max) n=n+1 T=-T*x*2/(2*n-2)*(2*n-3) COSX=COSX+Tend doIF (ABS(T)>EPS)THENPRINT*,'不能用此方法計算出結(jié)果'else print*,'本程序計算的余

41、弦結(jié)果是：',COSXprint*,'標(biāo)準(zhǔn)函數(shù)計算所得公式:',cos(x)end ifend結(jié)束3. 程序編寫4. 運行結(jié)果請輸入角度值（實數(shù)）:30 本程序計算的余弦結(jié)果是： 0.8660254 標(biāo)準(zhǔn)函數(shù)計算所得公式: 0.8660254Press any key to continue實驗四（1） 數(shù)組一 實驗?zāi)康?#160;1. 進(jìn)一步掌握數(shù)組聲明、數(shù)組元素引用、數(shù)組輸入輸出、數(shù)組賦初值、數(shù)組構(gòu)造器、隱含DO循環(huán)的使用規(guī)則。2. 掌握有關(guān)數(shù)組運算的常用標(biāo)準(zhǔn)函數(shù)。3. 掌握動態(tài)數(shù)組概念和動態(tài)數(shù)組的聲明、存儲分配的使用規(guī)則。4. 掌握設(shè)計和編寫多維數(shù)組和動態(tài)數(shù)組應(yīng)用

42、程序。 二 課堂實驗內(nèi)容 1 問題描述生成“魔幻方陣。魔幻方陣是指方陣的每一行、每一列、對角線元素之和均相等,其值為n(n2+1)/2,n代表方陣的行數(shù)和列數(shù)。一個n階奇數(shù)方陣由1到n2個自然數(shù)構(gòu)成,其組成規(guī)律是：(1)將1放在第1行的中間一列;(2)由m=1開始依次從左下到右上方向放連續(xù)數(shù),如m放在A(i,j),則m+1放在A(i-1,j+1);(3)如果i-1<1,則i-1為n;(4)如果j+1>n,則j+1為1;(5)如果A(i-1,j+1)已放數(shù),則將數(shù)放在A(i+1,j)中。以下是五階魔幻方陣：17 24 1 8 1523 5 7 14 16 4 6

43、13 20 2210 12 19 21 311 18 25 2 9  2 算法設(shè)計真假輸出：魔幻方陣A結(jié)束開始初始化：1Þi,3=jmn2i-1Þi,j+iÞj初始化：1ÞA(1,3), 2Þm真假 i<1nÞi真假 j>n1Þj真假 A(i,j)=0mÞA(i,j),m=m+1i+3Þi,j-2Þj真假 i>ni-nÞi真假 j<1j+nÞj圖7-1 程序框圖3. 程序編寫Ø !班級：工力112Ø !姓名：&#

44、216; !時間：201305Ø !生成魔幻方陣Ø PROGRAM mainØ PARAMETER (n=5)Ø INTEGER : A（：）=0,i=1,j=3,mØ A(i,j)=1Ø m=2Ø DO WHILE (m<=n*n)Ø i=i-1; j=j+1Ø IF (i<1) i=nØ IF (j>n) j=nØ IF (A(i,j)=0) THENØ A(i,j)=mØ m=m+1Ø ELSEØ i=i+3;j=j-2&

45、#216; IF (i>n) i=i-nØ IF (j<1) j=j+nØ End ifØ ENDDOØ DO i=1,nØ PRINT (<n>(I3,1X),(A(i,j),j=1,n)Ø ENDDOØ END4. 運行結(jié)果17 24 1 8 1523 5 7 14 16 4 6 13 20 2210 12 19 21 311 18 25 2 9   Press any key to continue三. 課外練習(xí)1 問題描述   已知兩組有序整數(shù)(個數(shù)不定

46、),將兩組整數(shù)合并為一組有序整數(shù)。合并后的該組整數(shù)沒有相同整數(shù)。從鍵盤輸入兩組整數(shù)，輸出合并前后整數(shù)序列。輸入使用靜態(tài)數(shù)組和動態(tài)數(shù)組兩種方式，編寫程序?qū)崿F(xiàn)之。提示：使用三個數(shù)組，兩個保存合并前兩組整數(shù)，一個保存合并后整數(shù)。2 算法設(shè)計3. 程序編寫!姓名：!班級：工力112!時間：201305!合并兩個有序數(shù)組program mainreal a(1:50),b(1:50),c(1:50)read(*,*)m,nread(*,*)(a(i),i=1,m),(b(j),j=1,n)i=1j=1k=1do while(（i.le.m）.and.（j.le.n)）if (a(i).gt.b(j)th

47、enc(k)=b(j)j=j+1else c(k)=a(i)i=i+1endifk=k+enddowrite(*,*)(c(k),k=1,m+n)end1，2，3，4，55，6，7，8，91，2，3，4，5，6，7，8，9Press any key to continue4. 運行結(jié)果實驗四(2) 語句函數(shù)一 實驗?zāi)康?#160;1. 了解函數(shù)的有關(guān)基本概念，如：函數(shù)名稱、函數(shù)自變量、函數(shù)調(diào)用、函數(shù)表達(dá)式等。2. 掌握語句函數(shù)定義、調(diào)用、參數(shù)的基本內(nèi)容和使用規(guī)則。3. 掌握使用語句函數(shù)設(shè)計和編寫程序的方法。二 課堂實驗內(nèi)容1 問題描述 求解一元二次方程2x2+3x-2=0的兩個實數(shù)根

48、，并輸出。使用語句函數(shù)編寫程序?qū)崿F(xiàn)之。 開始2 算法設(shè)計輸入方程系數(shù)A,B,C計算判別式D!姓名：!班級：工力112!時間：201305!語句函數(shù)計算方程兩根PROGRAM mainDISC(a,b,c)=SQRT(b*b-4*a*c)F(a,b,d)=(-b+d/(2*a)real A,B,C,D,X1,X2print*,'請輸入方程系數(shù)A,B,C'READ*A,B,CD=DISC(A,B,C)x1=F(A,B,D)x2=F(A,B,-D)print*,'x1=',x1print*,'x2=',x2語句函數(shù)計算根x1,x2輸出x1，x

49、2結(jié)束3. 程序編寫4. 運行結(jié)果 請輸入方程系數(shù)A,B,C2,3,-2 x1= 0.5000000 x2= -2.000000Press any key to continue三. 課外練習(xí)1 問題描述   求解以下問題：使用語句函數(shù)編寫程序，求解牛頓迭代法求方程x3-2x2+x-1=0在x=1.5附近的實數(shù)根。牛頓迭代公式有： x0=1.5 xi+1=xi-f(xi)/f¢(xi) I=1,2,3,其中：f¢(x)是f(x)的導(dǎo)數(shù)，當(dāng)|xi+1-xi|<10-7時，xi+1為方程的近似實數(shù)根。開始!姓名：!班級：工力112!時間：201305PROGRAM main IMPLICIT NONE REAL:ROOT,A,B,C,D,X0,X write( *,* ) "PLEASE INPUT A,B,C,D,X0" READ (*,*) A,B,C,D,X0 X=ROOT(A,B,C,D,X0) write(*,*) X return END PROGRAM MAIN FUNCTION ROOT(A,B,C,D,X0) RESULT(ROOT_RESULT) IMPLICIT NONE REAL,INTENT(IN):A,B,C,D,X0 REAL:ROOT_RESULT real:X1,