Fixed some typos that fucked up the basis functions. When the Jacobian was evaluated on the coordinates of the natural element it was different to the identity matrix.

Also added a code snipett that checks if the basis functions work correctly.