I am trying to solve an evolution equation (Allen - Cahn type) with ZERO Neumann boundary conditions over all boundaries. Function SOLVE, admits three arguments including variational problem, desired solution and bcs. Zero Neumann bcs are handled automatically using FEM, therefore I don't need to put the third arguments. Without the third arguments, Solve does not work. How can I resolve this problem?