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)