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 (R) code [https://arxiv.org/src/1907.05551/anc] 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.