The similarity between the 2D Helmholtz equation in elliptical coordinates and the Schrödinger equation for the simple mechanical pendulum inspires us to use light to mimic this quantum system. When optical beams are prepared in Mathieu modes, their intensity in the Fourier plane is proportional to the quantum mechanical probability for the pendulum. Previous works have produced a two-dimensional pendulum beam that oscillates as a function of time through the superpositions of Mathieu modes with phases proportional to pendulum energies. Here we create a three-dimensional pendulum wavepacket made of a superposition of Helical Mathieu-Gaussian modes, prepared in such a way that the components of the wave-vectors along the propagation direction are proportional to the pendulum energies. The resulting pattern oscillates or rotates as it propagates, in 3D, with the propagation coordinate playing the role of time. We obtained several different propagating beam patterns for the unbound-rotor and the bound-swinging pendulum cases. We measured the beam intensity as a function of the propagation distance. The integrated beam intensity along elliptical angles plays the role of quantum pendulum probabilities. Our measurements are in excellent agreement with numerical simulations.