8 node quadrilateral basis functions and calculation of the Jacobian.
Tests using sympy.