| 查看: 575 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者呆萌工科男将赠送您 30 个金币 | ||
呆萌工科男新虫 (正式写手)
|
[求助]
麻烦帮忙计算一个公式
|
|
|
我用mathematica计算一个积分的公式，可是算了很多次，很久没有结果，还报错，麻烦大神帮我码一下这个公式该怎么计算，谢谢了。公式如下： 我写的公式如下，可以直接帮我修改也行。  Untitled.png |
回复此楼» 猜你喜欢
普通院校药学硕士，做合成的，感觉找不到工作
已经有11人回复
-
求博导
已经有10人回复
-
Nature一直在编辑手里，考虑好几天了，是悬了吗
已经有10人回复
-
博士白读了
已经有41人回复
-
一直找不到审稿人
已经有6人回复
-
聊天
已经有12人回复
-
with editor 两个月了，什么原因？
已经有7人回复
-
研二光催化6月底四篇二区什么水平
已经有11人回复
-
某种做法不行。说过几遍了。同学还那样做。再那样做就给低分
已经有3人回复
-
两类问题算是白选了~
已经有11人回复
呆萌工科男
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 4644.9
- 散金: 19
- 红花: 2
- 沙发: 1
- 帖子: 522
- 在线: 137.8小时
- 虫号: 3956940
- 注册: 2015-07-06
- 性别: GG
- 专业: 声学
|
我写的mathematca如下： SystemOpen@ Export["a.dat", First@Cases[ Normal@Plot[-NIntegrate[-w Cos[ w t] Exp[-w^2/4]/2/(Log[2]/0.2^2) (BesselJ[0, 2 w] + Re[HankelH1[0, 2 w]/Pi/w/(-HankelH1[1, w/2])]* NIntegrate[(0.5 - 4 Cos[y]) BesselJ[1, Sqrt[0.25 + 16 - 4 Cos[y]] w]/ Sqrt[0.25 + 16 - 4 Cos[y]], {y, 0, Pi}] + Sum[Re[2 HankelH1[k, 2 w]/Pi/ w/(HankelH1[k, w/2] 2 k/w - HankelH1[k + 1, w/2])] NIntegrate[(0.5 - 4 Cos[y]) BesselJ[1, Sqrt[0.25 + 16 - 4 Cos[y]] w] Cos[k y]/ Sqrt[0.25 + 16 - 4 Cos[y]], {y, 0, Pi}], {k, 1, \[Infinity]}]), {w, 0, \[Infinity]}], {t, 0, 10}, PlotRange -> {-0.05, 0.12}] Line[t_] -> t, Infinity]] |
2楼2019-06-19 23:31:45
5