分段函数应用于计算时结果出现“True”

Question

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，而不是纯数值，我不知道怎么改，求教，谢谢。

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

Comments

1. 问题要尽己所能，抽象成阅读者最容易理解的语言；
2. 问题中涉及代码时，粘贴上文本代码远比图片有用（不然，你让阅读者自己去一个字母一个字母地敲你的代码来调试？），况且社区是支持代码语法高亮的。

---

抱歉，没发过，刚开始贴了代码，发现全在一行，觉得看着费劲，就改成图片了

---

代码语法高亮，参见右侧方框中第3条：在编辑器中点击“<>”（最后一个按钮）