红移和角距离的关系

def z_to_pc(z,mas):
    """
    :param z: 天体的红移，单位为km/(s*Mpc)
    :param mas: 望远镜的角分辨率
    :return: 望远镜在天体上所张的角距离，以pc为单位
    """
    H_0=74
    c = 3e5
    ang=mas/1000/60/60/180*3.1415926
    d_pc=((z+1)**2-1)/((z+1)**2+1)*c/H_0
    a_pc=d_pc*ang*1e6
    print(a_pc)
z_to_pc(0.059,mas=0.55)