[問題] 函數問題
請問各位先進，小弟寫了此函數計算，發現函數式子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
推
02/28 13:54, , 1F
02/28 13:54, 1F
→
02/28 13:55, , 2F
02/28 13:55, 2F
→
02/28 13:56, , 3F
02/28 13:56, 3F
→
02/28 18:16, , 4F
02/28 18:16, 4F
討論串 (同標題文章)
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
