Gravitational waves are produced by orbiting massive binary objects, such as black holes and neutron stars, and propagate as ripples in the very fabric of space-time. As the waves carry off orbital energy, the two bodies spiral into each other and eventually merge. They are described by Einstein's equations of general relativity. For the early phase of the orbit, called the inspiral, Einstein equations can be linearised and solved through analytical approximations, while for the late phase, near the merger, we need to solve the fully nonlinear Einstein's equations on supercomputers. In order to recover the gravitational wave for the entire evolution of the binary, a match is required between the inspiral and the merger waveforms. Our objectives are to establish an educational oriented toy model for a streamlined matching method, that will allow an analytical calculation of the complete gravitational waveform, while developing a gravitational wave modelling tutorial for undergraduate physics students. We use post-Newtonian (PN) theory for the inspiral phase, which offers an excellent training ground for students, and rely on Mathematica for our calculations, a tool easily accessible to undergraduates. For the merger phase we bypass Einstein's equations by using a simple analytic toy model named the implicit rotating source. After building the inspiral and merger waveforms, we construct our matching method and validate it by comparing our results with the waveforms for the first detection, GW150914, available as open source. Several future projects can be developed based from this project: building complete waveforms for all the detected signals, extending the PN model to take into account non-zero eccentricity, employing and testing a more realistic analytic model for the merger, building a separate model for the ringdown, and optimising the matching technique.