I'm writing a simple Maple program to test the Generalized Finite Element Method: main_screened_Poisso.mw
When trying to define the Neumann boundary conditions, I have to define a directional derivative dudn=dudx*n. However, I can't seem to define a unit vector normal to Gamma, which is defined by a LineSegments objects.
Other than that, the row reduction is very slow, even though I'm using floating point arithmatic and not exact arithmatic, I believe.
How can I solve these problems? Thanks in advance!