properly implement getting the Lagrange nodes.  Works without a priori information for all orders and dimensions

[[Imported from SVN: r8869]]