作者Jesuscrash (被愛是一種幸福)
看板Fortran
標題Re: [問題] 函數問題
時間Mon Feb 28 18:22:13 2011
※ 引述《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
1F：→ Semisphere:主程式要宣告F為函式，你要學會開implicit none 02/28 19:08
2F：推 youlieme:ijklim開頭的變數沒特別宣告 會當作整數用 02/28 21:20
3F：→ youlieme: mn 02/28 21:20
4F：→ Jesuscrash:感謝~已經解決了@@" 03/01 17:24