l*********i
发帖数: 483 | 1 用GSL的adaptive Runge-Kutta routine数值解出了f(x)在x=0到1的数值，
如果想要求解g(f(x),x)在区间x=0到1之间的积分，有什么现成的C/C++ routine可以
用么？因为最初f(x)是用adaptive的routine求解的，所以x的步长并不是一个
固定的值，所以一般的Newton-Cotes的库没法直接用。哪位大侠指点一下？ |
s**b
发帖数: 169 | 2 can you get arbitrary f(x) value?
【在 l*********i 的大作中提到】
: 用GSL的adaptive Runge-Kutta routine数值解出了f(x)在x=0到1的数值，
: 如果想要求解g(f(x),x)在区间x=0到1之间的积分，有什么现成的C/C++ routine可以
: 用么？因为最初f(x)是用adaptive的routine求解的，所以x的步长并不是一个
: 固定的值，所以一般的Newton-Cotes的库没法直接用。哪位大侠指点一下？
|
l*********i
发帖数: 483 | 3 不行，因为解f(x)的code是别人的，我这里只有一个data file.
【在 s**b 的大作中提到】
: can you get arbitrary f(x) value?
|
s**b
发帖数: 169 | 4 then you can use interpolation or do a fit.
assuming the curve is smooth.
【在 l*********i 的大作中提到】
: 不行，因为解f(x)的code是别人的，我这里只有一个data file.
|
l*********i
发帖数: 483 | 5 多谢。
【在 s**b 的大作中提到】
: then you can use interpolation or do a fit.
: assuming the curve is smooth.
|
