mathematical model for aura borialis