Re: [Libmesh-users] Non-homogeneous Neumann Boundary Condition

[Ted K. <ted...@go...>]

Thanks guys but the alterations in code that you gave are not working
either.

This is the actual equation I'm trying to solve:

d/dx (x * du/dx) = 2/x^2                   on  the domain   1 <  x < 2

subject to these boundary conditions:

u(1) = 2
(-x * du/dx) = 0.5   @ x = 2

Could you possibly show me how to solve this in libmesh with exemplar code.
This would really help me get through the boook I'm using to learn the FEM
using libmesh.

Thank you.

Regards

Ted