  Here we construct approximate analytical forms for the metric coefficients
and fields representing the scalarized Einstein-Maxwell black holes with
various couplings of the scalar field, once the parameters of the system are
fixed. By increasing approximation order, one can obtain the analytic
representation with any desired accuracy, what was tested via calculations of
shadows for these black holes by using approximate analytical and accurate
numerical metric functions. We share the Mathematica code which allows one to
find an appropriate analytical form of the metric for any couplings and values
of parameters. Scalarization increases the radius of the black-hole shadow for
all the considered coupling functions.
