红移距离转换也是常用的一个操作。现把代码放在下面

function y=zD(z) % z-distance relation, distance in units of Mpc
h=0.72;
DH=3000./h; % in Mpc
OmegaM=0.274;
OmegaLam=0.726;
%fun=strcat('',num2str(OmegaM),'*(1+x)^3+',num2str(OmegaLam),'')
DC=DH.*comdis(OmegaM,OmegaLam,z);
DM=DC; % for zero curvature
y=DM.*(1+z);


%=================================================

function y=comdis(OmegaM,OmegaLam,z0)
function f=E(z)
f=1.0./sqrt(OmegaM.*(1+z).^3+OmegaLam);
end
y=quadl(@E,0,z0);
end

转载本文请联系原作者获取授权，同时请注明本文来自钱磊科学网博客。
链接地址：https://m.sciencenet.cn/blog-117333-301651.html

上一篇：笔记（二）几个matlab脚本
下一篇：『读文献』（二）长时标射电瞬变源与中子星