Main program to solve deformation of DE balloon. from Delayed electromechanical instability of a viscoelastic dielectric elastomer balloon

This main program together with two subprograms are used to solve the evolution of deformation field of a dielectric elastomer balloon subject to fixed amount of air and fixed voltage.