Need Help on Matlab code PLEASE

This code is using the 4th Order Runge-Kutta Method to solve for the orbit of Halley's comet over time. I don't know why it won't work. I'm given the initial x0 vector and h is the step size. Any help would be appreciated.