In this work, the quasistatic thermoviscoelastic thermistor problem isconsidered. The thermistor model describes the combination of the effects due tothe heat, electrical current conduction and Joule's heat generation. The variationalformulation leads to a coupled system of nonlinear variational equations for whichthe existence of a weak solution is recalled.Then, a fully discrete algorithm is introduced based on the finite elementmethod to approximate the spatial variable and an Euler scheme to discretizethe time derivatives. Error estimates are derived and, under suitableregularity assumptions, the linear convergence of the scheme is deduced.Finally, some numerical simulations are performed in order to show the behaviourof the algorithm.