diff --git a/README.md b/README.md
index 7e7d315b439a0f10b805909d20581c4420bae44c..0417c981d39848790fe1a7d26a35e99d0991cf58 100644
--- a/README.md
+++ b/README.md
@@ -72,5 +72,10 @@ To characterize the use of the AS5013 as a three-dimensional tracking device, we
 <img src='as5013-test/testing2.jpg' width=350px>
 <img src='as5013-test/nrf52-as5013-board.jpg' width=400px>
 
+To maximize resolution (and minimize stroke), we are also considering alternative magnet arrangements.  I wrote a very simple (and slow!) 2D B field finite difference simulation, which is in <a href='as5013-test/sim/'>sim</a>.  It uses successive overrelaxation to compute the spatial distribution of the z component of the magnetic field.  Below we should output of the simulation for a single magnet, as well as for a pair of opposing magnets.  The question is whether we can create a region of high field gradient over each of the four sensors.  I think a central rod magnet in one orientation surrounded by four with the opposite orientation would work nicely for this.
+
+<img src='as5013-test/sim/run_single_magnet.png' width=350px>
+<img src='as5013-test/sim/run_double_magnet.png' width=350px>
+