Modelling a soft iron disk

[M4a1x]

Hi there!

Great library! It has been a pleasure to try out different things and build different magnet assemblies (or Collections).

I have recently stumbled across this paper here A simple and low-cost permanent magnet system for NMR, where the authors use two simple disc magnets and homogenize the field by inserting two "Pure Iron" (SUY) disks between the two magnets made out of Y30BH.

I was wondering: Would it be possible to simulate this using magpylib and magpylib-material-response? Or would I need to buy or find a FEM tool for this?
My idea would be to vary the polarization angle and strength of the disk magnets to see what kind of effect I could expect from the iron.

Looking at your soft magnet tutorial it looks like I could simply make a cylinder with no polarity? Is this correct? what would I use for the susceptibility?

iron_disk = magpy.magnet.Cylinder(polarization=(0, 0, 0), dimension=(0.12, 0.02))
iron_disk.susceptibility = 3999  # µr=4000

I hope this makes sense. Thanks for reading!

Cheers,
Max