347.sagews

y, x = var ('x,y')
plot3d(lambda x,y : (bessel_Y(y,x,"scipy",53)).real(), (-0.002,0.002),(-1,1))

implicit_plot3d(lambda x,y,z: z-(bessel_Y(y,x,"scipy",53)).real(), (-0.002,0.002),(-1,1), (-100,100))