A numerical modeling of Earth's atmosphere is carried out by means of the Transmission Line Matrix (TLM) numerical method with the aim of calculating the Schumann resonance frequencies. The numerical results obtained are very close to the experimental ones and those obtained with the widely accepted two-scale-height ionospheric model, which allows us to affirm that this is a valid numerical tool for predicting the Schumann frequencies in the atmospheres of other planets and moons. The great flexibility of the TLM numerical method also allows the study of slight shifts in the Schumann resonance frequencies due to an increase in electrical conductivity at the Earth's poles, originating from high-energy particle precipitation emitted from the Sun in conjunction with solar flares. A slight increase in the first Schumann resonance frequency is observed during these events, which is associated with a reduction in the dimensions of the electromagnetic cavity.