這問題錯大了,由此基本可以認(rèn)為你沒有系統(tǒng)學(xué)習(xí)過 Fortran 的任何語法.
建議你買一本教材先看看.
隨便選幾個錯誤：
第一：Function 不能寫在主程序內(nèi)部,程序單元必須獨(dú)立.
也就是：
主程序
結(jié)束主程序
函數(shù) 1
結(jié)束函數(shù)1
函數(shù)2
結(jié)束函數(shù)2
這樣的結(jié)構(gòu).不能把函數(shù)寫在其他函數(shù)的 “肚子” 里.
第二：M1,W1,M2,W2 這 4 個函數(shù).函數(shù)倒是給出了各自的虛參.但是調(diào)用時,沒有給出任何對應(yīng)的實(shí)參.
你確定你的程序需要使用 function 老師上課水的要命，還要求用子程序計(jì)算，你可不可以教教我怎么修改？我不是學(xué)力學(xué)的，專業(yè)上的東西，我?guī)筒涣四?。我只能給你從格式和結(jié)構(gòu)上幫助你：以下的代碼，我不能保證算出來的是正確的。因?yàn)槟愕拇a問題太多，很多運(yùn)算方法，我不確定。比如你的代碼里有 W2=(F8L**3)/(48*EI)我根本不知道 F8L 是什么東西。Q 在這個函數(shù)里沒有出現(xiàn)過，卻在虛參表里？以下代碼，需要你根據(jù)自己的算法進(jìn)行校對。Program MainImplicit None !// 這個十分關(guān)鍵，每個程序單元都要加Real Q,L,F,EI,M,WReal M1 , W1 , M2 , W2 !// 所有變量都要定義Read( * , * ) Q,F,L,EIM=M1(Q,L)+M2(Q,L)!// 調(diào)用函數(shù)，必須給出參數(shù)W=W1(Q,L,EI)+W2(Q,L,EI) !// 調(diào)用函數(shù)，必須給出參數(shù)PRINT *,"最大彎矩值M=",MPRINT *,"最大撓度值W=",WEnd Program MainReal Function M1(Q,L)Implicit None !// 這個十分關(guān)鍵，每個程序單元都要加Real Q , L !// 所有變量都要定義M1=(1./8.)*Q*L**2!//// 八分之一要寫成 1.0/8.0 或者 1./8.，如果寫成 1/8，表示是整數(shù)，整數(shù)的 1/8 = 0!!!!End Function M1Real Function W1(Q,L,EI)Implicit None !// 這個十分關(guān)鍵，每個程序單元都要加Real Q , L , EI !// 所有變量都要定義W1=(5.*Q*L**4)/(384.*EI)End Function W1Real Function M2(Q,L)Implicit None !// 這個十分關(guān)鍵，每個程序單元都要加Real Q , L !// 所有變量都要定義M2=(1./4.)*Q*L!// 四分之一也是End Function M2Real Function W2(Q,L,EI)Implicit None !// 這個十分關(guān)鍵，每個程序單元都要加Real Q , L , EI !// 所有變量都要定義W2=(Q*8*L**3)/(48*EI)End Function W2