Question: How do I define a unit vector normal to a LineSegments objects

December 27 2013 Adriaan 10


I'm writing a simple Maple program to test the Generalized Finite Element Method:

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!


Please Wait...