Proton und Deuteron

mp = 1.6726231*10^(-27); me = 9.1093897*10^(-31); 
    h = 6.626076*10^(-34); ϵ0 = 8.854188*10^(-12); 
    e = 1.602177*10^(-19); 
    md = 2*mp; μp = me*mp/(me + mp); μd = me*md/(me + md); 
    Epd = 3*10^6*e; V0 = 10*10^6*e; a = 10^(-14);

Proton:

Tp = (1 - 
            Epd/V0)/((1 - Epd/V0) + V0/4/Epd*
              Sinh[(8*π^2*mp (V0 - Epd)/h^2)^(1/2)*
                    a]^2)

Tp2 = 
      16*Epd/V0^2*(V0 - Epd)*
        Exp[-2*(8*π^2*mp (V0 - Epd)/h^2)^(1/2)*
            a]

Deuteron:

Td = (1 - 
            Epd/V0)/((1 - Epd/V0) + V0/4/Epd*
              Sinh[(8*π^2*md (V0 - Epd)/h^2)^(1/2)*
                    a]^2)

Td2 = 
      16*Epd/V0^2*(V0 - Epd)*
        Exp[-2*(8*π^2*md (V0 - Epd)/h^2)^(1/2)*
            a]