To illustrate discord and its use, Zurek in his paper on discord (NB pdf) gives example of a partially decohered bell state i.e. $$\rho_{AB}=\frac{1}{2}(|00\rangle\langle 00|+|11\rangle\langle 11|) + \frac{z}{2}(|00\rangle \langle 11|+|11\rangle \langle 00|)$$ using the measurement basis for Alice's side. $$\{\cos\theta |0\rangle + e^{i\phi} \sin\theta |1\rangle,\ e^{-i\phi} \sin\theta |0\rangle + \cos\theta |1\rangle\}$$ He plotted it in the Fig.1.
Now my problem is that I have not been able to reproduce this result. Here is what I got. I can not seem to locate what could be the source problem. The Mathematica code which I used is here. Can anyone see what might be going wrong here.