This research studies the application of the nth order Runge-kutta method as a numerical solution to ordinary differential equations. This method was chosen because it is able to provide high accuracy and flexibility in various PDB problems. We implement the nth-order Runge-Kutta algorithm in MATLAB and compare with other numerical methods, such as Euler's method. The results show that the nth order Runge-Kutta method is able to produce more accurate solutions, especially for nonlinear systems. This research makes a significant contribution to the development of numerical solutions for PDB and shows the potential of MATLAB as an effective tool for numerical simulation. Sensitivity analyzes of parameters and time steps were also performed to understand the impact of variations on stability and convergence.