找回密码
 注册
Simdroid-非首页
查看: 138|回复: 4

[数值计算] 用mathematica计算一个广义积分的困惑

[复制链接]
发表于 2013-3-29 22:45:20 | 显示全部楼层 |阅读模式 来自 广东深圳
本帖最后由 szldh2005 于 2013-3-30 11:40 编辑

现有一个广义积分。并用变量代换的方法计算其等价积分及其截断积分(使用NIntegrate函数及其缺省参数)。发现,当u<10e-14时,有的计算结果相差太大,不知原因何在?非常困惑。计算结果见附件。红色结果为异常计算结果。(maple计算的结果非常好)

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2013-3-30 12:30:23 | 显示全部楼层 来自 陕西西安
Simdroid开发平台
NIntegrate[ 1/Log[y], {y, 0, Exp[-10^-34]}, WorkingPrecision -> 100]

Integrate[ 1/Log[y], {y, 0, Exp[-10^-34]}] // N[#, 20] &

maple的确很强,有空帮忙看看maple版的这个帖子吧:

http://forum.simwe.com/thread-1072065-1-1.html
回复 不支持

使用道具 举报

 楼主| 发表于 2013-3-30 18:34:36 | 显示全部楼层 来自 广东深圳
版主:您好!非常感谢!&的作用是什么?不加&,计算结果是∞。
回复 不支持

使用道具 举报

 楼主| 发表于 2013-3-30 22:35:38 | 显示全部楼层 来自 广东深圳
本帖最后由 szldh2005 于 2013-3-30 23:51 编辑

谢谢版主!已搞明白&的作用,但N[Integrate[]与 Integrate[]//N[] & 的区别在哪里?前者计算第2个积分为∞。
回复 不支持

使用道具 举报

 楼主| 发表于 2013-4-2 11:28:36 | 显示全部楼层 来自 广东深圳
已清楚,区别在于N[expr]及N[expr,n]上。不指定精度时,结果为∞,指定精度时可得到数值结果。但为什么?看帮助文件也未弄明白。
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-24 16:22 , Processed in 0.031071 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表