nii John Smith 6 years ago I need only the MATLAB codes to produce the figures identically that in the published paper. You can use Runge _Kutta (RK) with shooting method or Finite difference method.