我算了一个很简单的求积分的程序，得到的结果为什么是0呢？用imsl下的数据库函数QDAGS (F, a, b, RES, ERRA,ERRREL, ERREST)求的。高手帮忙指点下吧。不论函数f(x)取什么形式，得到结果总是0.
program jifenQDAGS
use imsl
implicit none
real a,b,erra,errest,f,res,errrel
external f
ERRA=1.0e-3
ERRREL=1.0e-2
a=0.0
b=3.0
CALL QDAGS (F, a, b, RES, ERRA,ERRREL, ERREST)
WRITE (*,*) RES, ERREST
END
REAL FUNCTION F (X)
REAL X
F = x*x
RETURN
END |
回复此楼
+1/83