add the definition of FemNodalFunctionMatrix
for now disable this part.
define the matrix to store the derivative and initialize them with zeros.
this is where the derivative should be calculated.