This paper presents the numerical analysis for a variational formulation of rate-independent phase transformations in elastic solids due to Mielke et al. The new model itself suggests an implicit time-discretization which is combined with thefinite element method in space. A priori error estimates are established for the quasioptimal spatial approximation of the stress fieldwithin one time-step. A posteriorierror estimates motivate an adaptive mesh-refining algorithm for efficient discretization. The proposed scheme enables numerical simulations which show that the model allows for hysteresis.