  Parker Solar Probe observations reveal that the near-Sun space is almost
filled with magnetic switchbacks (``switchbacks'' hereinafter), which may be a
major contributor to the heating and acceleration of solar wind. Here, for the
first time, we develop an analytic model of an axisymmetric switchback with
uniform magnetic field strength. In this model, three parameters control the
geometry of the switchback: height (length along the background magnetic
field), width (thickness along radial direction perpendicular to the background
field), and the radial distance from the center of switchback to the central
axis, which is a proxy of the size of the switchback along the third dimension.
We carry out three-dimensional magnetohydrodynamic simulations to investigate
the dynamic evolution of the switchback. Comparing simulations conducted with
compressible and incompressible codes, we verify that compressibility, i.e.
parametric decay instability, is necessary for destabilizing the switchback.
Our simulations also reveal that the geometry of the switchback significantly
affects how fast the switchback destabilizes. The most stable switchbacks are
2D-like (planar) structures with large aspect ratios (length to width),
consistent with the observations. We show that when plasma beta ($\beta$) is
smaller than one, the switchback is more stable as $\beta$ increases. However,
when $\beta$ is greater than one, the switchback becomes very unstable as the
pattern of the growing compressive fluctuations changes. Our results may
explain some of the observational features of switchbacks, including the large
aspect ratios and nearly constant occurrence rates in the inner heliosphere.
