A Fortran program has been developed to optimize the design of coil configurations of axisymmetric magnets with an iron flux return yoke for stray field shielding such as those used in the areas of nuclear magnetic resonance, nuclear physics, and industrial processes. Some results and comments on the practical performance of the code are presented. The discussion covers the accuracy of the basic field routines, choice of points of interest, convergence of the optimization, influence of discretization of iron, successive-overrelaxation performance and numerical noise, and finding the final configuration.