calc-rho   A few short examples for calling rho()
