Using variable-order fractional derivatives in differential equations is essential. It enables more precise modeling of complex phenomena with varying memory and long-range dependencies, improving our ability to describe real-world processes reliably. This study investigates the properties of solutions for a two-point boundary value problem associated with φ-Caputo fractional derivatives of variable order. The primary objectives are to establish the existence and uniqueness of solutions, as well as explore their stability through the Ulam-Hyers concept. To achieve these goals, Banach’s and Krasnoselskii’s fixed point theorems are employed as powerful mathematical tools. Additionally, we provide numerical examples to illustrate results and enhance comprehension of theoretical findings. This comprehensive analysis significantly advances our understanding of variable-order fractional differential equations, providing a strong foundation for future research. Future directions include exploring more complex boundary value problems, studying the effects of varying fractional differentiation orders, extending the analysis to systems of equations, and applying these findings to real-world scenarios, all of which promise to deepen our understanding of Caputo fractional differential equations with variable order, driving progress in both theoretical and applied mathematics.