We numerically simulate the dynamical relaxation to equilibrium of a sunspot. Treating the sunspot as a thick axisymmetric flux tube in cylindrical geometry, we solve the time dependent MHD equation to examine the evolution ...