NUMERICAL BGK-SOLUTIONS OF LARGE SCALE ELECTROSTATIC POTENTIAL IN AURORAL PLASMAS