Create calculateKEff.m

fcn to calculate the local stiffness matrix