Renamed MJ to CJ for Cauchy (and also TJ for tau)
Added derivation using sympy