Finding and Plotting Lorenz Solution using MATLAB

I wrote a function, LorenzRK4IVP(), that takes the system of three differential equations as input and solves the system using the Runge-Kutta method with step size . I plot the strange attractor as well as use MATLAB to produce a GIF of the solution.