Possible error in exponentiation when computing force for stiffness

[gregouar]

According to the paper, when computing force in CalcStiffness, we should take
force<-(dat$Chant*dat$Cwant)^(3/2)
instead of
force<-(dat$Chant*dat$Cwant)^3/2