2332.sagews

x,y,z=var('x,y,z')
s1=implicit_plot3d(x^2+y^2+(z/2)^2==1, (x,-2,2), (y,-2,2), (z,-2,2))
s2=point3d((2,0,0),size=5,color='red')
s3=line3d([(2,0,0),(1,0,0)],color='green')
show(s1+s2+s3)