公告:1)网站程序升级:Q2A升级到1.8.6,Wordpress升级到5.7.2
2)修复了头像加载慢与提交问题反应慢等问题
2021-06-16

欢迎来到 Mathematica 问答社区

提问时请贴上文本代码

语法高亮:在编辑器中点击

被禁止的话题:广告破解

请阅读:《提问的智慧》

备用域名:mma.ooo

支持LaTex数学公式
行内公式标识符:\$ 或“$\backslash ($”+“$\backslash )$”,
行间公式标识符:\$\$ 或 “$\backslash [$”+“$\backslash ]$”

社区建议QQ群:365716997

分类

–1 投票
356 浏览
M15[M_]=M/10^15;
mevap[M_]:=10^(-5)*M15[M]^(-2);
mk[M_,z]:=Which[z<=10^13,10^(-58)*(1+z)^4*M15[M]^2,z>10^13,0];
Tb[z_]:=Piecewise[{{2*(1+z),z>=170},{2*(1+z)^2/171},z<170}];
mb36[M_,z_]:=Piecewise[{{10^(-20)*(1+z)^1.5*M15[M]^2,M<=10^33*(1+z)^(-1.5)},{0.07*M15[M],M>=10^33*(1+z)^(-1.5)}}];
mb33[M_,z_]:=Piecewise[{{10^(-17)*M15[M]^2,M<=10^33},{0.07*M15[M],M>10^30}}];
mrb[M_,z_]:=Piecewise[{{10^7*10/86*(1+z)/171^4*M15[M],z>3000},{10^(-45)*(10^(-5)*9*10^(-30)*(1+z)^4/(10^(-5)*9*10^(-30)*(1+z)^4)+10^3*10^(-30)*(1+z)^3/4900)^1.5*(10^(-5)*10^(-30)*(1+z)^4+2206*10^(-30)*(1+z)^3/4900)*M^2,1000<z<3000},{10^(-48)*(1+z)^4*M^2/10^30+mb36[M,z],173<z<1000},{10^48*(1+z)^4*M^2/10^30+mb33[M,z],z<173}}]
mrb[10^33,2000]
LogLogPlot[mrb[10^33,z],{z,10^(-4),10^17}]



 

求一个比较罗嗦的函数mrb[M,z],然后求积分。中间过程用了Piecewise和Which写分段函数,分段函数画图的时候不显示0的值,但是单独求某个自变量的函数值没问题,但是到最后用分段函数的时候就显示结果中有True,而不是纯数值,我不知道怎么改,求教,谢谢。

部分源码,大概框架有了,省略了部分函数数值的定义。

分类:函数 | 用户: 大胖 (31 分)
修改于 用户:大胖
1. 问题要尽己所能,抽象成阅读者最容易理解的语言;
2. 问题中涉及代码时,粘贴上文本代码远比图片有用(不然,你让阅读者自己去一个字母一个字母地敲你的代码来调试?),况且社区是支持代码语法高亮的。
抱歉,没发过,刚开始贴了代码,发现全在一行,觉得看着费劲,就改成图片了
代码语法高亮,参见右侧方框中第3条:在编辑器中点击“<>”(最后一个按钮)

登录 或者 注册 后回答这个问题。

...