triangle ([ 0 0 ], [ 1 0 ], [ 0 1 ]) Jacobian = [ 1 0 ] [ 0 1 ] Jacobi determinant = 1 inverse Jacobian = [ 1 -0 ] [ -0 1 ] Laplace elMat = [ 1 -0.5 -0.5 ] [ -0.5 0.5 0 ] [ -0.5 0 0.5 ] mass elMat = [ 0.0833333 0.0416667 0.0416667 ] [ 0.0416667 0.0833333 0.0416667 ] [ 0.0416667 0.0416667 0.0833333 ] element load vector = [ 0.166667 0.166667 0.166667 ] triangle ([ 1 1 ], [ 1.5 1 ], [ 1.25 1.5 ]) Jacobian = [ 0.5 0.25 ] [ 0 0.5 ] Jacobi determinant = 0.25 inverse Jacobian = [ 2 -1 ] [ -0 2 ] Laplace elMat = [ 0.625 -0.375 -0.25 ] [ -0.375 0.625 -0.25 ] [ -0.25 -0.25 0.5 ] mass elMat = [ 0.0208333 0.0104167 0.0104167 ] [ 0.0104167 0.0208333 0.0104167 ] [ 0.0104167 0.0104167 0.0208333 ] element load vector = [ 0.0416667 0.0416667 0.0416667 ]