# Problem with Neumann boundary conditions for discontinuous Galerkin method

**URL:** <https://community.freefem.org/t/problem-with-neumann-boundary-conditions-for-discontinuous-galerkin-method/2325>\
**Category:** General Discussion\
**Created:** [March 12, 2023, 1:43pm UTC](https://community.freefem.org/t/problem-with-neumann-boundary-conditions-for-discontinuous-galerkin-method/2325 "2023-03-12T13:43:20Z")\
**Posts on this page:** 1\
**Page:** 1

<div class="post-metadata">

**Author:** ![prabir.barman](https://avatars.discourse-cdn.com/v4/letter/p/e0b2c6/32.png) [@prabir.barman](https://community.freefem.org/u/prabir.barman)\
**Post date:** [March 12, 2023, 1:43pm UTC](https://community.freefem.org/t/problem-with-neumann-boundary-conditions-for-discontinuous-galerkin-method/2325/1 "2023-03-12T13:43:20Z")

</div>

I am trying to impose flux boundary conditions (Neumann type) in the 2D heat conduction equation. The int1d(Th,3)(…) works fine for the Galerkin method, but this integration does not affect boundary values in discontinuous Galerkin.  
Please help me to impose the flux boundary condition for Discontinuous Galerkin in FreeFem++.