Differential equations are powerful tools for understanding and predicting the behavior of dynamical systems in nature, engineering, and society. A dynamical system is a system with a state that evolves over time, such as a pendulum, disease spread, or weather. These systems are governed by mathematical rules, often expressed as differential equations. While algebraic equations solve for numbers, differential equations solve for functions and their derivatives. For example, $ f'(x) = f(x) $ has the solution $ e^x $. This chapter introduces solving first-order ordinary differential equations and systems using the Forward Euler method, demonstrated through population growth and disease spread. It then explores second-order differential equations for mechanical oscillations, where the Forward Euler method is less effective, leading to the development of more advanced methods like Euler-Cromer, Runge-Kutta, and finite difference schemes. The chapter begins with population growth, where the differential equation $ N'(t) = rN(t) $ models exponential growth. The solution is $ N(t) = N_0 e^{rt} $, derived from the difference equation $ N(t+\Delta t) = N(t) + \Delta t r N(t) $. The Forward Euler method approximates this as $ N^{n+1} = N^n + \Delta t r N^n $, which is accurate for small $ \Delta t $. The chapter also discusses more realistic models, such as the logistic model where growth rate depends on population size, leading to a carrying capacity. The SIR model for disease spread involves three categories: susceptible (S), infected (I), and recovered (R), with differential equations describing transitions between these states. The Forward Euler method is applied to these systems, with parameters like $ \beta $ and $ \gamma $ determining disease spread and recovery rates. The chapter emphasizes the importance of numerical methods for solving differential equations, especially when analytical solutions are not feasible. It includes programming examples in Python, demonstrating how to implement the Forward Euler method for both scalar and system differential equations. Verification techniques, such as consistency checks, ensure the accuracy of numerical solutions. The chapter concludes by highlighting the power of mathematical modeling, where specific examples lead to generalizable tools applicable to various problems. The SIR model, for instance, can be adapted to different diseases and populations, illustrating the broad applicability of differential equations in understanding complex systems.Differential equations are powerful tools for understanding and predicting the behavior of dynamical systems in nature, engineering, and society. A dynamical system is a system with a state that evolves over time, such as a pendulum, disease spread, or weather. These systems are governed by mathematical rules, often expressed as differential equations. While algebraic equations solve for numbers, differential equations solve for functions and their derivatives. For example, $ f'(x) = f(x) $ has the solution $ e^x $. This chapter introduces solving first-order ordinary differential equations and systems using the Forward Euler method, demonstrated through population growth and disease spread. It then explores second-order differential equations for mechanical oscillations, where the Forward Euler method is less effective, leading to the development of more advanced methods like Euler-Cromer, Runge-Kutta, and finite difference schemes. The chapter begins with population growth, where the differential equation $ N'(t) = rN(t) $ models exponential growth. The solution is $ N(t) = N_0 e^{rt} $, derived from the difference equation $ N(t+\Delta t) = N(t) + \Delta t r N(t) $. The Forward Euler method approximates this as $ N^{n+1} = N^n + \Delta t r N^n $, which is accurate for small $ \Delta t $. The chapter also discusses more realistic models, such as the logistic model where growth rate depends on population size, leading to a carrying capacity. The SIR model for disease spread involves three categories: susceptible (S), infected (I), and recovered (R), with differential equations describing transitions between these states. The Forward Euler method is applied to these systems, with parameters like $ \beta $ and $ \gamma $ determining disease spread and recovery rates. The chapter emphasizes the importance of numerical methods for solving differential equations, especially when analytical solutions are not feasible. It includes programming examples in Python, demonstrating how to implement the Forward Euler method for both scalar and system differential equations. Verification techniques, such as consistency checks, ensure the accuracy of numerical solutions. The chapter concludes by highlighting the power of mathematical modeling, where specific examples lead to generalizable tools applicable to various problems. The SIR model, for instance, can be adapted to different diseases and populations, illustrating the broad applicability of differential equations in understanding complex systems.