I have also found some difference in calculating force in axisymetric model: when calculating repulsion forces between two cylindrical blocks, one, a superconductor in mixed state, and the other, a NdFeB magnet, I've found different results using FEMM and student's Quickfield. I don't know how to explain this difference. Agnaldo.