the electric field in a homogeneous dispersive nonisotropic medium. Assuming the pulse is propagating in the direction of the z-axis, it can be shown that the envelope for one of the most general of cases, namely a biaxial crystal, is governed by the PDE: where the coefficients contains diffraction and dispersion effects which have been determined analytically with computer algebra and verified numerically to within third order for both isotropic and non-isotropic media, valid in the near-field and far-field. is the inverse of the group velocity projection. The term in is the group velocity dispersion (GVD) or second-order dispersion; it increases the pulse duration and chirps the pulse as it propagates through the medium. The term in is a third-order dispersion term that can further increase the pulse duration, even if vanishes. The terms in and describe the walk-off of the pulse; the coefficient is the ratio of the component of the group velocity and the unit vector in the direction of propagation of the pulse (z-axis). The terms in and describe diffraction of the optical wave packet in the directions perpendicular to the axis of propagation. The terms in and containing mixed derivatives in time and space rotate the wave packet about the and axes, respectively, increase the temporal width of the wave packet (in addition to the increase due to the GVD), increase the dispersion in the and directions, respectively, and increase the chirp (in addition to that due to ) when the latter and/or and are nonvanishing. The term rotates the wave packet in the plane. Oddly enough, because of previously incomplete expansions, this rotation of the pulse was not realized until the late 1990s but it has been experimentally confirmed. To third order, the RHS of the above equation is found to have these additional terms for the uniaxial crystal case: The first and second terms are responsible for the curvature of the propagating front of the pulse. These terms, including the term in are present in an isotropic medium and account for the spherical surface of a propagating front originating from a point source. The term can be expressed in terms of the index of refraction, the frequency and derivatives thereof and the term also distorts the pulse but in a fashion that reverses the roles of and (see reference of Trippenbach, Scott and Band for details). So far, the treatment herein is linear, but