Representation Formulas and Stability Analysis for Hilfer–Hadamard Proportional Fractional Differential Equations

This paper introduces a novel version of the Gronwall inequality specifically related to the Hilfer–Hadamard proportional fractional derivative. By utilizing Picard’s method of successive approximations along with the definition of Mittag–Leffler functions, we derive a representation formula for the solution of the Hilfer–Hadamard proportional fractional differential equation featuring constant coefficients, expressed in the form of the Mittag–Leffler kernel. We establish the uniqueness of the solution through the application of Banach’s fixed-point theorem, leveraging several properties of the Mittag–Leffler kernel. The current study outlines optimal stability, a new Ulam-type concept based on classical special functions. It aims to improve approximation accuracy by optimizing perturbation stability, offering flexible solutions to various fractional systems. While existing Ulam stability concepts have gained interest, extending and optimizing them for control and stability analysis in science and engineering remains a new challenge. The proposed approach not only encompasses previous ideas but also emphasizes the enhancement and optimization of model stability. The numerical results, presented in tables and charts, are provided in the application section to facilitate a better understanding.