var('x,y,z')
plott1 = implicit_plot3d(x^2+y^2+z^2==4,(x,-2,2),(y,-2,2),(z,-2,2),opacity=0.4)
plott2 = plot3d(1-x-y,(x,-2,2),(y,-2,2),opacity=0.6,color='red')
plott2 += point3d((1/3,1/3,1/3),size=20)
plott2 += point3d((0,0,0),size=10)
show(plott1+plott2,aspect_ratio=1)