Nonunitary equal-spin pairing magnetization demonstration
========================================================

Pairing:
  Delta_upup(k,T)     = Delta(T) [sin(kx)+i sin(ky)]
  Delta_downdown(k,T) = r Delta(T) [sin(kx)+i sin(ky)]

Target density = 1.1008995627
Representative T/Tc = 0.050000

r = 1.00
  M_z = 0.0000000000 mu_B/site
  <q_Delta,z> = 0.0000000000
  normalized pair spin = 0.0000000000

r = 0.60
  M_z = 0.0066836100 mu_B/site
  <q_Delta,z> = 0.1151998810
  normalized pair spin = 0.4705882353

r = 0.25
  M_z = 0.0123211965 mu_B/site
  <q_Delta,z> = 0.1687498256
  normalized pair spin = 0.8823529412

r = 0.00
  M_z = 0.0144335081 mu_B/site
  <q_Delta,z> = 0.1799998140
  normalized pair spin = 1.0000000000
