Index out of bounds in Runge-Kutta algorithm for parabolic shot

Question

Hi I am trying to get this program to run. The z structure is assumed to be size 4, however it gives the error. I also notice that if I print z in 1 it gives me a value, so I do