2419.sagews

v=cremona_optimal_curves([1..100])

w = [E.plot(thickness=1,
   rgbcolor=(random(),random(),random())) for E in v]

graphics_array(w, 10, 10).show(axes=False)