# Orthocircles

bb = 3; ff = 0.075; i = -2; j = 2; ContourPlot3D[ ((x^2 + y^2 - 1)^2 + z^2)*((y^2 + z^2 - 1)^2 + x^2)*((z^2 + x^2 - 1)^2 + y^2) - ff^2*(1 + bb*(x^2 + y^2 + z^2)) == 0, {x, i, j}, {y, i, j}, {z, i, j}, PlotRange -> All, PlotPoints -> 100, Axes -> False, Boxed -> False, ContourStyle -> Directive[RandomColor[], Opacity[0.7], Specularity[White, 30]]]

The equation is derived by placing 3 circles perpendicular to each other.