MAXX = 10 # we specify what ranges we wanna look at for x and y
MINX = -4
MAXY = 20
MINY = 0
x = var('x')
p = plot( sqrt(x^3 + 16) , (x,MINX,MAXX)) # plot
# we add some artificial "integer coordinate points" by drawing little circles
for xv in range(MINX,MAXX+1):
for yv in range(MINY,MAXY+1):
p += circle((xv,yv),0.04,color='red') # draws a red circle centered at (xv,yv) of given radius
p.show(xmin=MINX, ymin=MINY, xmax=MAXX, ymax=MAXY, figsize=6)