var('y z')

(y, z)

R=1
r=0.25

t1=((x-1.5)^2+y^2+z^2+R^2-r^2)^2-4*R^2*((x-1.5)^2+y^2)
t2=((x-1.5*cos(2*pi/3.0))^2+(y-1.5*sin(2*pi/3.0))^2+z^2+R^2-r^2)^2-4*(R^2)*((x-1.5*cos(2*pi/3.0))^2+y^2)
t3=((x-1.5*cos(4*pi/3.0))^2+(y-1.5*sin(4*pi/3.0))^2+z^2+R^2-r^2)^2-4*(R^2)*((x-1.5*cos(4*pi/3.0))^2+y^2)
t=t1*t2*t3-10

t

(((x - 1.50000000000000*cos(1.33333333333333*pi))^2 + (y - 1.50000000000000*sin(1.33333333333333*pi))^2 + z^2 + 0.937500000000000)^2 - 4*(x - 1.50000000000000*cos(1.33333333333333*pi))^2 - 4*y^2)*(((x - 1.50000000000000*cos(0.666666666666667*pi))^2 + (y - 1.50000000000000*sin(0.666666666666667*pi))^2 + z^2 + 0.937500000000000)^2 - 4*(x - 1.50000000000000*cos(0.666666666666667*pi))^2 - 4*y^2)*(((x - 1.50000000000000)^2 + y^2 + z^2 + 0.937500000000000)^2 - 4*(x - 1.50000000000000)^2 - 4*y^2) - 10

implicit_plot3d(t,(x,-6,6),(y,-6,6),(z,-6,6),plot_points=120,color='lemonchiffon',viewer='threejs',online=True,frame=False)