ON A NUMERICAL METHOD FOR THE ZEEMAN-STARK MODEL OF LIGAND BINDING UNDER EM EXPOSURE