コード例 #1
0
    res = 300
    chi = (0, py.pi)
    phi = (0, 2 * py.pi)
    v = 1.000001 / n  #
    w = 0.9 * f.MaxW(n, m, v)  # Example values

n, m, v, w, chi, phi, a, proj, axes = p.VarProd(res,
                                                variables=(n, m, v, w, chi,
                                                           phi, a))

ch = f.ChAngle(n, m, v, w)
conc = f.Concurrence(n, m, v, w, chi, phi, a)
conc2 = f.Concurrence(n, m, v, w, chi, py.pi, a)

if sphere is True:
    k = 1.1 * p.Delta(ch - axes[0], 0.02) * p.Delta(axes[1], 0.02) * py.amax(z)
    if type(z) is not tuple:
        z = (z, )
    z = z + (k, )

labels = ('$C(\chi)$', '$\\theta_C$', '$|\mathcal{M}_{1,2-}|^2$',
          '$|\mathcal{M}_{1,2+}|^2$')
fig, ax = p.MultPlot((chi, ),
                     conc,
                     title='',
                     labels=labels,
                     xlabel='$\chi$',
                     ylabel='',
                     labelsize=25,
                     projection='polar',
                     sphere=sphere)