应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Python » scipy.integrate.quad 积分错误原因？
51/1返回列表
查看: 1846  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

swuli

木虫 (小有名气)

[求助] scipy.integrate.quad 积分错误原因？ 已有1人参与

有以下积分，n2和n3的值非常接近，但积分结果却大相径庭。
n2=22144  积分结果=0.4256854383924181
n3=22145 积分结果=11774635265351.027

用matlab算，后面是准确的。是什么原因导致n2所出现的错误？


from math import sqrt, log, exp
import numpy as np
from scipy.integrate import quad

def integrand(x, r1, r2):
    a = r1 + r2
    b = r1 * r2
    R = a * x / 2
    R1 = R ** 2 - a ** 2
    R2 = R ** 2 - (r1 - r2) ** 2
    vdw = -7e-21 * (2 * b / R1 + 2 * b / R2 + log(R1 / R2))
    edl = 7.8e-12 * b / a * log(1 + exp(-328774227 * (R - a)))
    force = vdw + edl
    return exp(force / 4.11447e-21) / x ** 2

n1 = 5010
n2 = 22144
n3 = 22145
radius1 = 1.05 * 10 ** -8 * n1 ** (1 / 1.8)
radius2 = 1.05 * 10 ** -8 * n2 ** (1 / 1.8)
radius3 = 1.05 * 10 ** -8 * n3 ** (1 / 1.8)
w2 = quad(integrand, 2, np.inf, args=(radius1, radius2))
w3 = quad(integrand, 2, np.inf, args=(radius1, radius3))
print('w2=', w2[0])
print('w3=', w3[0])
回复此楼

» 猜你喜欢
1楼 2021-05-21 15:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Certificate


发自小木虫Android客户端
2楼2021-05-21 16:59:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bcsnow

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
swuli: 金币+10, ★★★★★最佳答案 2021-05-23 08:45:33
没有错误啊：w2= 11772040939067.406  w3= 11774635265351.027
赞一下 回复此楼
3楼2021-05-22 10:36:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swuli

木虫 (小有名气)
引用回帖:
3楼: Originally posted by bcsnow at 2021-05-22 10:36:12
没有错误啊：w2= 11772040939067.406  w3= 11774635265351.027

我这里得到的结果却是w2= 0.4256854383924181；w3= 11774635265351.027。奇怪了。
和版本有关吗？我的python 3.6.6 scipy 1.6.0。您有没有做什么别的处理呢？
赞一下 回复此楼
4楼2021-05-23 08:48:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bcsnow

铁杆木虫 (著名写手)
引用回帖:
4楼: Originally posted by swuli at 2021-05-23 08:48:12
我这里得到的结果却是w2= 0.4256854383924181；w3= 11774635265351.027。奇怪了。
和版本有关吗？我的python 3.6.6 scipy 1.6.0。您有没有做什么别的处理呢？...

没有处理，直接拷贝运行的，我的是3.7.7和1.5.2，按说版本间不应该有这么严重错误。你把n2 n3的值调换一下，看看结果如何。
赞一下 回复此楼
5楼2021-05-24 00:23:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 swuli 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程调剂 +3 咪咪空空 2026-03-12 3/150 2026-03-12 16:36 by Equinoxhua
[基金申请] 面上和青基一样限30页不合理 +4 wowsunflower 2026-03-10 6/300 2026-03-12 15:59 by Equinoxhua
[考研] 一志愿中科大材料与化工。353分能调剂到哪？ +5 MayUxw1 2026-03-08 10/500 2026-03-12 15:46 by EBSD
[考研] 302求调剂 +4 负心者当诛 2026-03-11 4/200 2026-03-12 15:03 by 每天只摆一小会
[考研] 311求调剂 +7 zchqwer 2026-03-10 7/350 2026-03-12 14:48 by ruiyingmiao
[考研] 274求调剂0856材料化工 +9 z2839474511 2026-03-11 10/500 2026-03-11 23:46 by BruceLiu320
[考研] 289求调剂 +4 步川酷紫123 2026-03-11 4/200 2026-03-11 19:41 by peike
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 7x4iz20zqa 2026-03-11 8/400 2026-03-11 19:32 by s9xj0iuolc
[考研] 材料工程330分求调剂，一志愿985 +5 小材化本科 2026-03-07 5/250 2026-03-11 16:51 by 535743368
[考研] 材料专硕288分求调剂 一志愿211 +3 在家想你 2026-03-11 3/150 2026-03-11 08:13 by ms629
[考研] 265求调剂 +8 小木虫085600 2026-03-09 11/550 2026-03-10 19:44 by 绿色低碳化学
[考研] 328，0703考生求调剂，一志愿为东北师范大学 +3 观素律 2026-03-09 4/200 2026-03-10 15:58 by houyaoxu
[考研] 08工科 +5 li李乐成 2026-03-06 5/250 2026-03-10 14:51 by 在风落中
[考研] 求调剂 一志愿苏州大学，0856化工323分 | 本科应化 | 有专利/竞赛/科研助手经历 | +7 橙子cyx 2026-03-06 9/450 2026-03-10 12:03 by kiokin
[考研] 320求调剂 +4 魏zy 2026-03-08 4/200 2026-03-09 16:14 by ruiyingmiao
[考研] 求调剂，数一英一274分 +4 小菲会努力 2026-03-08 4/200 2026-03-09 12:40 by 一定上岸哟_
[考研] 理学287分求调剂 +4 看看我. 2026-03-05 4/200 2026-03-08 22:11 by qingfeng258
[考研] 085701环境工程317分求调剂 +9 6汆尼9 2026-03-07 9/450 2026-03-08 06:41 by 刘兵
[考研] 276求调剂 本科天津工业大学 一志愿北京化工大学0855 +4 路lyh123 2026-03-06 6/300 2026-03-07 11:37 by shdliugang
[考研] 材料工程调剂 +6 zbcm_zbcm 2026-03-05 6/300 2026-03-06 15:36 by @飒飒飒飒
信息提示
关闭
请填处理意见
关闭