大俗人
2007-05-04, 02:53
小弟最近遇到一个数值积分问题，请教一下：
mm=(sqrt(s^2+795648849807027/737869762948382064640)/s-1)*sin(4258168138844955/1125899906842624*s)*cos(2036772065443477/562949953421312*s)
这个表达式要从0积到无穷。
问题是：
1。这个式子在0时得值是无穷。
2。我选择用命令：
maple('int',mm,'s=0.01..5'),
或者
vpa(int(qa(i),'s',0.01,5))；
来积分，花很长得时间，结果报错，而且matlab重启之前不能用。
3，据说可以用quad来积分，但是那个函数里面的相应参量是一个函数而不是一个表达式，不知该怎么处理，因为我的这个式子是用一个循环算出来的，有120多个，所以每次手工定义一个相关的函数不是很现实。
忘高手指教。
还有，我最后是用maple 软件算出来的，两个软件之间不知道能不能互相调用？
mm=(sqrt(s^2+795648849807027/737869762948382064640)/s-1)*sin(4258168138844955/1125899906842624*s)*cos(2036772065443477/562949953421312*s)
这个表达式要从0积到无穷。
问题是：
1。这个式子在0时得值是无穷。
2。我选择用命令：
maple('int',mm,'s=0.01..5'),
或者
vpa(int(qa(i),'s',0.01,5))；
来积分，花很长得时间，结果报错，而且matlab重启之前不能用。
3，据说可以用quad来积分，但是那个函数里面的相应参量是一个函数而不是一个表达式，不知该怎么处理，因为我的这个式子是用一个循环算出来的，有120多个，所以每次手工定义一个相关的函数不是很现实。
忘高手指教。
还有，我最后是用maple 软件算出来的，两个软件之间不知道能不能互相调用？