The present research correlates with a fuzzy hybrid approach merged with a homotopy perturbation transform method known as the fuzzy Shehu homotopy perturbation transform method (SHPTM). With the aid of Caputo and Atangana–Baleanu under generalized Hukuhara differentiability, we illustrate the reliability of this scheme by obtaining fuzzy fractional Cauchy reaction–diffusion equations (CRDEs) with fuzzy initial conditions (ICs). Fractional CRDEs play a vital role in diffusion and instabilities may develop spatial phenomena such as pattern formation. By considering the fuzzy set theory, the proposed method enables the solution of the fuzzy linear CRDEs to be evaluated as a series of expressions in which the components can be efficiently identified and generating a pair of approximate solutions with the uncertainty parameter λ∈[0,1]. To demonstrate the usefulness and capabilities of the suggested methodology, several numerical examples are examined to validate convergence outcomes for the supplied problem. The simulation results reveal that the fuzzy SHPTM is a viable strategy for precisely and accurately analyzing the behavior of a proposed model
