Use c code in getcovzpd.m. Make it faster under octave.