Age-specific dose coefficients are required to assess internal exposure to the general public. This study utilizes reference age-specific biokinetic models of iodine to estimate total numbers of nuclear disintegrations ã(rS,τ) occurring in source regions (rS) during the commitment time (τ). Age-specific S values are estimated for 35 target regions due to 131I present in 22 rS using 10 paediatric reference computational (PRC) phantoms (representing 5 ages for both sexes) data recently published by International Commission of Radiation Protection (ICRP). Monte Carlo transport simulations are performed in FLUKA code. The estimated ã(rS,τ) and S values are then used to compute committed tissue equivalent dose HT(τ) for 27 radiosensitive tissues and dose coefficients e(τ) for all 5 ages due to inhalation and ingestion of 131I. The derived ã(rS,τ) values in the thyroid source are observed to be increasing with age due to increased retention of iodine in the thyroid. S values are found to be decreasing with age mainly due to increase in target masses. Generally, HT(τ) values are observed to decrease with age, indicating the predominant behaviour of S values over ã(rS,τ) . On an average, ingestion dose coefficients are 63% higher than inhalation for all ages. The maximum contribution to dose coefficients is from thyroid, accounting for 96% in the case of new-born and 98-99% for all other ages. Furthermore, the estimated e(τ) values for the reference public are observed to be lower than previously published values by the ICRP. The estimated S, HT(τ), and e(τ) values can be used to improve the estimation of internal doses to organs/whole body for members of the public in cases of 131I exposure. The estimated dose coefficients can also be interpolated for other ages to accurately evaluate the doses received by the general public during 131I therapy or during a radiological emergency.