add some mathematica code for analytical derivations