Just finished the implementation of RK4 method to my vacuum projectile motion. Since acceleration is constant, the RK4 method is under utilized, and so produces a exact replica of the heun graph. This is to be expected since constant acceleration leads to RK2 and RK3 being the same quantity, and therefore the formula simply becomes heun’s method. I plotted two graphs, one was the trajectory of the 3 methods, and you can clearly see euler’s method starts to drift away slowly from the heun and RK4 curve. This is even more clear in the error graph with log axis. The large curve of the euler’s error shows how quickly the error accumulates, as compared to the rough curve of the heun and RK4 method. This roughness is noise caused by the program’s mathematical rounding limits. The error is so small that the program is not able to properly round it. (floating point arithmetic noise). AFter this, I will most likely be implementing drag.