In 1909 Theodor Wulf made the first detections of cosmic rays. Over a century later, the origin of cosmic rays is still an open problem for physicists. This problem is due to the difficulty of determining the trajectories these particles made on their way to Earth. This is why the study of cosmic ray trajectories is such an interesting and important field. In the early 1900’s Adriaan Fokker and Max Planck developed a partial differential equation to describe the probability distribution function of particles experiencing both Brownian motion as well as outside forces. This equation can be applied to cosmic rays trajectories through turbulent magnetic fields. When applied to this situation the equation does not yet have a known simple analytical solution. Usually to solve this problem numerical methods are used, however, numerical methods are not practical solutions for everyone who needs to use these solutions in their work. We have therefore used Legendre polynomials to approximate a solution to the Fokker-Plank equation. The first, second and third order solutions were found as well as some of the solutions’ expectation values. These expectation values are useful both in future analytic work as well as in experimental work. These analytic solutions were compared to a numerical solution to show that for very early times the third order approximation is significantly better than the second order approximation, but at late times only a second order Legendre expansion is needed to achieve good results. This solution is so close to the numerical solution will allow people to more easily track the trajectories of cosmic rays.