  General analytic arguments lead us to expect that in the modified dynamics
(MOND) self-gravitating disks are more stable than their like in Newtonian
dynamics. We study this question numerically, using a particle-mesh code based
on a multi-grid solver for the (nonlinear) MOND field equation. We start with
equilibrium distribution functions for MOND disk models having a smoothly
truncated, exponential surface-density profiles and a constant Toomre $Q$
parameter. We find that, indeed, disks of a given ``temperature'' are locally
more stable in MOND than in Newtonian dynamics. As regards global instability
to bar formation, we find that as the mean acceleration in the disk is lowered,
the stability of the disk is increased as we cross from the Newtonian to the
MOND regime. The degree of stability levels off deep in the MOND regime, as
expected from scaling laws in MOND. For the disk model we use, this maximum
degree of stability is similar to the one imparted to a Newtonian disk by a
halo three times as massive at five disk scale lengths.
