Added Magnetic Polarization code

Added Error (difference) plot