code: def newton_refine1(s, new_s, curve1, t, new_t, curve2):
      
      if NO_IMAGES:
          return
  
      points = np.hstack([curve1.evaluate(s), curve2.evaluate(t)])
      points_new = np.hstack([curve1.evaluate(new_s), curve2.evaluate(new_t)])
      ax = curve1.plot(256)
      curve2.plot(256, ax=ax)
      ax.plot(
          points[0, :],
          points[1, :],
          color="black",
          linestyle="None",
          marker="o",
          markeredgewidth=1,
          markerfacecolor="None",
      )
      ax.plot(
          points_new[0, :],
          points_new[1, :],
          color="black",
          linestyle="None",
          marker="o",
      )
      ax.axis("scaled")
      save_image(ax.figure, "newton_refine1.png")
docs: Image for :func:`.newton_refine` docstring.