Fixed zernike.m vectorization. Was using matrix multiplication between radial and azimuthal parts - this has been changed to pointwise multiplication.