關於fortran 77的題目 有誰會?? 20點

這個問題有誰會??? 需用fortran 77這個軟體(規定 ~"~)

可以的話 麻煩順便解釋各行的語法(老師會隨機問>""<

問題

以半間距法,計算苯在25C(298.15 K)下的飽和液態與飽和氣態的莫耳體積,該溫度下的壓力由安東尼方程式計算。苯的熱力物理性質:臨界溫度Tc為562.2 K、臨界壓力Pc為48.98 bar,其安東尼常數分別為A=13.8594、B=2773.78和C=220.07。

相關公式與方程式

1.van der Waals equation體積立方型方程式:

(1). 方程式型態一 :

(2). 方程式型態二:

(3). 方程式型態三:

(4). 或得自行推導適當型式之方程式

2.立方型狀態方程式之物質參數a、b的取得:

(1).整理,(2).用臨界壓力與(3).臨界溫度表示: 、

3.安東尼方程式

(1).方程式型態一

(2).方程式型態二

已更新項目:

我需要的只是"".for""黨 麻煩會寫的大大 寫成.for文字檔給我 >"<

拜託囉

裡面的語法需採用 fortran77 去寫

2 個已更新項目:

謝謝那位大大的提醒

小弟沒注意到方程式沒貼上去 >"<

我直接PO上題目

http://eagle90071.myweb.hinet.net/Q2.doc

1 個解答

評分
  • 1 0 年前
    最佳解答

    C MAIN CALCULATION,BY BISECTION METHOD

    DX = UPPER - LOW

    ISTEP = 0

    WARN=0

    X0=0.D0

    DO WHILE (ABS(DX).GT.ERROR)

    X0=(UPPER+LOW)/2.D0

    IF ((VDW(A,B,P,T1,LOW)*VDW(A,B,P,T1,X0)).LT.0.D0) THEN

    UPPER = X0

    DX = UPPER-LOW

    ELSE

    WARN=WARN+1

    LOW = X0

    DX = UPPER-LOW

    END IF

    ISTEP = ISTEP + 1

    END DO

    P.S.

    0.想知道為什麼你的發問後來都被移除?

    1.這裡只PO上最核心的計算。

    2.字數限制的關係,把註解都刪了。

    3.這個CORE只可計算BOUNDARY內恰有一根狀況。

    參考資料: ”Numerical Mathematics and Computing 3rd Edition”,Ward Cheney & David Kincaid,1994,Brooks/Cole Publishing Company
還有問題?馬上發問,尋求解答。