answered 11 years ago
The following seems to be what you are describing.
var('y') f=sqrt(25-y^2) g=-sqrt(25-y^2) plot([f,g],(y,-5,5),aspect_ratio=1)
Note that for numerical reasons, this will often not display well at y=±5.