\[\rho(s,t)=(s-t)(st)^{\gamma-1}\left((1-s)(1-t)\right)^{\delta-1}\*\left((1-(s/%
a))(1-(t/a))\right)^{\epsilon-1},\]