作者Jesuscrash (被爱是一种幸福)
看板Fortran
标题[问题] 函数问题
时间Mon Feb 28 02:53:54 2011
请问各位先进，小弟写了此函数计算，发现函数式子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
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
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.125.202.117
1F：推 terryys:1. 这能compile吗？应该至少要把F也宣告为real 02/28 13:54
2F：→ terryys:2. F(x)里2969/20000等都是整数相除，所以都是0。 02/28 13:55
3F：→ terryys:要把这些都变成浮点数才行。 02/28 13:56
4F：→ Jesuscrash:谢大大~我试试看@@~ 02/28 18:16