Despite the COVID-19 pandemic and wearing masks in many countries, women are keen on elegance, beauty and the use of face foundations. Assessment of health risks associated with the regular use of face foundation by females is dynamic due to the emerging products. The most common international 14 brands of face foundation powders were collected and the concentrations of different elements (Ag, Al, As, B, Ba, Be, Ca, Cd, Co, Cr, Cu, Fe, Hg, K, Li, Mg, Mn, Mo, Na, P, Pb, Sb, Se, Sn, V and Zn) in each sample were determined. A combined approach merging the conventional and computational tools was used for investigating the risk of exposure to toxic elements. Monte Carlo simulations were applied to calculate risks associated with twenty elements. We attempted different probability distribution functions for concentrations because the actual distribution functions are not known, and the only data available are the mean value and standard deviation of concentrations obtained from experiment. Our results indicate that the total non-carcinogenic health risk through exposure to different elements (Hazardous Index, HI) does not strongly depend on the choice of the probability distribution function for the concentrations. We also show that taking into account probability distributions of other variables and parameters such as body weight, exposed skin area, skin adhesion, etc. does not significantly change the main result rather just slightly broadening the final Hazardous Index distribution function. We found that calculated HI is well below unity for all considered samples, i.e., the dermal exposure to toxic elements in the considered facial powders is negligible and the considered face foundation powders are quite safe to use.