Unlike the universal approximation theorems for functions mapping from a real-valued (RV) vector to a RV number or from a complex-valued (CV) vector to a CV number, in the field of electromagnetism, we need to approximate functions mapping from a RV vector to a CV number when we consider the electric field as a function of the spatial coordinate in the frequency domain. Typically, CV numbers contain phase information. When such phase information is handled properly, the performance of the neural networks (NNs) can be improved. This work derives a universal approximation theorem for functions mapping from a RV vector to a CV number. A deep NN, named as HV-DL, is designed accordingly, which consists of a RV input layer, a RV module containing two branches, a CV module and a CV output layer. The proposed universal approximation theorem is verified by numerical experiments on the HV-DL solution of the two-dimensional (2-D) electric field integral equation (EFIE). To integrate the underlying physics of electromagnetic scattering into the proposed HV-DL, the loss function is computed according to the EFIE.