Numerical calculation for quantum discord in matlab