Hello,<br>         <br>         I am writing a 3D FEM code for learning purpose (experimenting with object oriented concepts in Fortran 2003). Can some one tell me (pseudo code) how to implement a non homogenous Neumann boundary condition. You can also point me to a book. I am using tetrahedral elements for solving  Poisson equation. I am using a 4 point Gauss quadrature, and linear basis functions for cell volume integrals. I am confused on how to do the integration on the facet. <br>
<br>Thanks<br>Reddy<br>        <br clear="all"><br>-- <br>-----------------------------------------------------<br>Dharmendar Reddy Palle<br>Graduate Student<br>Microelectronics Research center,<br>University of Texas at Austin,<br>
10100 Burnet Road, Bldg. 160<br>MER 2.608F, TX 78758-4445<br>e-mail: <a href="mailto:dharmareddy84@gmail.com">dharmareddy84@gmail.com</a><br>Phone: +1-512-350-9082<br>United States of America.<br><br>