Neumann_polynomial loading ...