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 ------------------------------------------------------------------------------ Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day trial. Simplify your report design, integration and deployment - and focus on what you do best, core application coding. Discover what's new with Crystal Reports now. http://p.sf.net/sfu/bobj-july _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
