  We consider Einstein-Weyl gravity with a minimally coupled scalar field in
four dimensional spacetime. By using the Minimal Geometric Deformation (MGD)
approach, we split the highly nonlinear coupled field equations into two
subsystems that describing the background geometry and scalar field source,
respectively. Regarding the Schwarzschild-AdS metric as a background geometry,
we derive analytical approximate solutions of scalar field and deformation
metric functions with Homotopy Analysis Method (HAM), providing their
analytical approximations to fourth order. Moreover, we discuss the accuracy of
the analytical approximations, showing they are sufficiently accurate
throughout the exterior spacetime.
