We have simulated the magnetic relaxation and the formation of magnetic domains in very thin films with strong anisotropy perpendicular to the film plane. We have used Monte Carlo simulation based on a two-dimensional classical Ising model on a square lattice including the long-range dipole-dipole interaction and an external magnetic field. Starting with all magnetic moments aligned, the magnetic relaxation exhibits two distinct behaviors, depending upon the relation $\ensuremath{\alpha}$ between the dipole-dipole and the exchange interactions. For $\ensuremath{\alpha}l{\ensuremath{\alpha}}_{c}\ensuremath{\cong}0.37$, the magnetization follows an exponential decay, $\frac{M(t)}{M(0)}=\mathrm{exp}[\ensuremath{-}\frac{t}{\ensuremath{\tau}({T}^{\ensuremath{'}})}]$ ($\ensuremath{\tau}$ is the relaxation time and ${T}^{\ensuremath{'}}$ the reduced temperature) and the domain pattern at the beginning of the relaxation process is characterized by the nucleation of a few domains, followed by a rapid growth of the magnetic domain size. For $\ensuremath{\alpha}l{\ensuremath{\alpha}}_{c}$, the magnetization follows a power-law time decay, $\frac{M(t)}{M(0)}=a(\ensuremath{\alpha}, {T}^{\ensuremath{'}}){t}^{\ensuremath{-}\ensuremath{\gamma}({T}^{\ensuremath{'}})}$, with a demagnetization process associated with the nucleation of many domains at random positions in the film. In both cases, the system relaxes towards the ground state which, depending on the value of $\ensuremath{\alpha}$, exhibits a striped structure. With the present model we were able to obtain the energy of domain nucleation, ${E}_{N}$, for $\ensuremath{\alpha}l{\ensuremath{\alpha}}_{c}$, despite the complexity associated with long-range dipole-dipole interactions; we have also obtained the $\ensuremath{\alpha}$ dependence of the energy of domain nucleation, ${E}_{N}(\ensuremath{\alpha})$.