Gravitational decoupling via the Minimal Geometric Deformation (MGD) approach has been used extensively in General Relativity (GR), mainly as a simple method for generating exact anisotropic solutions from perfect fluid seed solutions. Recently this method has also been used to generate exact spherically symmetric solutions of the Einstein-scalar system from the Schwarzschild vacuum metric. This was then used to investigate the effect of scalar fields on the Schwarzschild black hole solution. We show that this method can be extended to higher order theories. In particular, we consider fourth order Einstein–Weyl gravity, and in this case by using the Schwarzschild metric as a seed solution to the associated vacuum field equations, we apply the MGD method to generate a solution to the Einstein–Weyl scalar theory representing a hairy black hole solution. This solution is expressed in terms of a series using the Homotopy Analysis Method (HAM).