First, consider the approximation of the potential field 
by a
two dimensional Taylor's series expansion about a point (x, y):
where h and k are the maximum x and y distances within an element respectively. Using the first two terms (up to first order terms) in the above Taylor's expansion, we can obtain the standard linear interpolation function for a triangle:
where A is the area of the triangle.
