For many, mathematics is just the language of numbers. But, what people forget is that the universe is made of numbers, thus mathematics is its language. That is why we used a mathematical approach to study the long-term behavior of a moving ball influenced by the Earth’s gravitational field. For our purpose, we experimentally defined gravitational billiards to be vertical billiards in which the earth gravitational field affects the billiards ball’s motion. We also assumed that the billiards boundary could be any mathematical function of our choice, and the experimental system was conservative. Once in motion, the billiards ball followed a two-dimensional projectile trajectory until collision with the boundary occurred. After that, the ball engaged in a new projectile trajectory. As a matter of fact, each trajectory between two points of ball-boundary collision could be considered as an isolated two-dimensional gravitational motion with its own set of conditions. We took advantage of that aspect to build different simulations in parabolic, circular, paraboloidal and spherical gravitational billiards. The long-term behavior of the ball under gravity in our different gravitational billiards presented some aspects that so far corroborated the physics behind our study.