  In this paper, we construct anisotropic spherical solutions from known
isotropic solutions through extended gravitational decoupling method in the
background of self-interacting Brans-Dicke theory. The field equations are
decoupled into two sets by applying geometric deformations on radial as well as
temporal metric components. The first array corresponds to isotropic fluid
while influence of the anisotropic source is confined to the second set. The
isotropic sector is determined through metric functions of isotropic solutions
(Tolman IV/Krori-Barua) whereas two constraints on the anisotropic source are
required to close the second system. The impact of the massive scalar field as
well as the decoupling parameter on the physical characteristics of the
anisotropic solutions is analyzed graphically. We also check the viability,
compactness, surface redshift and stability of the obtained solutions. It is
found that the resulting solutions follow accepted physical trend for some
values of the decoupling parameter.
