Solving the Dirac equation is a formidable task due to the high frequency and the degrees of freedom involved. However, this high frequency allows one to obtain an approximation to the equation. Here, we directly solve the Dirac equation using an envelope method and derive analytical solutions of Dirac wave packets to first order for the small momentum spread. We apply the insight gained from this solution to the Zitterbewegung behavior in a Dirac-like system, where we show that it is crucial to include the first-order term in our solution to correctly describe a Dirac packet.