PTT數位生活區 / Fortran

Re: [問題] 函數問題

看板Fortran作者 (被愛是一種幸福)時間13年前 (2011/02/28 18:22), 編輯推噓1(1推 0噓 3→)
留言4則, 3人參與, 最新討論串3/4 (看更多)
※ 引述《Jesuscrash (被愛是一種幸福)》之銘言： : 請問各位先進，小弟寫了此函數計算，發現函數式子debug似乎要帶入 : 的數字並沒有帶入函數式子中，希望各位先進可以能給個幫忙~謝謝 : program INTEGRAL : Real A,B,W1,W2,W3,W4,W5,W6,X1,X2,X3,X4,X5,X6,I,X,PI : PI=3.141593 : A=0.0 : B=1.0 : W1=0.0337652429 : W2=0.1693953068 : W3=0.380690407 : W4=0.619309593 : W5=0.8306046932 : W6=0.9662347571 : X1=W1*B+W6*A : X2=W6*B+W1*A : X3=W2*B+W5*A : X4=W5*B+W2*A : X5=W3*B+W4*A : X6=W4*B+W3*A : I=(B-A)*(0.0856622462*(F(X1)+F(X2))+ : < 0.1803807865*(F(X3)+F(X4))+ : < 0.2339569673*(F(X5)+F(X6))) : WRITE(*,*) I : END PROGRAM : function F(X) : implicit none : real X,DYDX,T,R,F : T=12 : DYDX=5*T*(2969/20000/(x**0.5)-63/500-879/1250*x+8529/10000*(x**2)- : > 203/500*(x**3)) : F=SQRT(1+DYDX**2) : RETURN : END 我將F浮點數補上後compile所得到的答案仍然是1.0，還是怪怪的@@~ 懇請各位幫個忙> <~謝謝!! -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.125.202.117
→
02/28 19:08, , 1F
主程式要宣告F為函式，你要學會開implicit none
02/28 19:08, 1F
推
02/28 21:20, , 2F
ijklim開頭的變數沒特別宣告 會當作整數用
02/28 21:20, 2F
→
02/28 21:20, , 3F
mn
02/28 21:20, 3F
→
03/01 17:24, , 4F
感謝~已經解決了@@"
03/01 17:24, 4F
‣ 更多 Jesuscrash 的文章
文章代碼(AID): #1DQtTMq9 (Fortran)
討論串 (同標題文章)
以下文章回應了本文：
完整討論串 (本文為第 3 之 4 篇)：
排序： 最舊先 | 最新先 | 留言數
在新視窗開啟完整討論串 (共4篇)
Fortran 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 Jesuscrash 的文章
文章代碼(AID): #1DQtTMq9 (Fortran)