Listed below is the source code for compiling the two torus link.
The command (package): <<Graphics`ParametricPlot3D` loads in the function: ParametricPlot3D.


<<Graphics`ParametricPlot3D`;
htorus[R_, r_, psi_, phi_] :=

{(R+r Cos[psi]) Cos[phi], (R+r Cos[psi]) Sin[phi], r Sin[psi]};
horiz = ParametricPlot3D[
Evaluate[htorus[6,1,psi,phi]],
{phi, 0, 2Pi, 2Pi/18}, {psi, 0, 2Pi, 2Pi/8},
Axes -> False,
Boxed -> False,
DisplayFunction -> Identity];
vtorus[R_, r_, psi_, phi_] :=
{r Sin[psi], (R+r Cos[psi]) Sin[phi]+3, (R+r Cos[psi]) Cos[phi]-1};
vertic = ParametricPlot3D[
Evaluate[vtorus[6,1,psi,phi]],
{phi, 0, 2Pi, 2Pi/18}, {psi, 0, 2Pi, 2Pi/8},
Axes -> False,
Boxed -> False,
DisplayFunction -> Identity];
Show[horiz, vertic, DisplayFunction -> $DisplayFunction];