Long time dissipative MHD simulations have been done of high beta, large aspect ratio tokamaks with bean shaped cross-section, resembling the PBX experiment. A disruption free path from the first to the second regimes is followed explicitly. The amplitudes of long wavelength resistive MHD modes increase with beta up to a critical value. As beta increases above this value, the mode amplitudes decrease. If the pressure and current profiles are not broad enough, long wavelength modes produce a disruption at moderate beta, which rapidly flattens the pressure profile. In the computational model, time dependent equilibria are produced by a balance of MHD quasilinear effects, transport, a pressure source and the toroidal electric field. Dissipation is provided by Spitzer resistivity and perpendicular thermal conduction. A rigid conducting wall boundary condition is used, but there is a highly resistive edge layer, which simulates a thin vacuum region. The numerical method makes use of conformal coordinates, which are generated by prescribing a set of external coil currents.