This code uses Newton's Method for implicit problems.
I recommend checking solutions with ODE45 MATLAB function (which aproximates the solutions with Runge-Kutta4). I also have a RK4 solver in my other directory that can be used for parameter checking.