Convergence to equilibrium for the Modified Phase Field Crystal equation (see Ref. 20) 
βu_{tt}+u_{t}=Δ[∆^{2}u+ 2∆u+u^{3}+(1ε)u]

on the square (0,6π)x(0,6π) with periodic boundary conditions

β=0.1 / ε=2 / time step dt=0.25 / 200 iterations

initial value u_{0}(x,y)=0.2+0.2cos(x)cos(y) and u_{t}(0)=0

P1P1P1 Finite element discretization in space

2nd order discretization in time (scheme of Gomez and Hughes)

Software: FreeFem++
