作者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