%time
w = 1000
P = sum(point((j,j*(w-j)), color=(j/w,1-j/w,0)) for j in range(w))

%time P.save('a.png')

show(P)