[vg] = groupVelocity(obj, fWavelength) [vgO, vgE] = groupVelocity(obj, fWavelength, fTheta) [vgS, vgF] = groupVelocity(obj, fWavelength, fTheta, fPhi)