stuff = [line(L,rgbcolor=(1,1/4,1/2)) for L in [[f(a,b*pi/20) for a in range(-5,4)] for b in range(-21,20)]]
new_plot=stuff[0]
for j in stuff[1:]:
new_plot+=j
stuff = [line(L,rgbcolor=(0,1/4,1/2)) for L in [[f(a,b*pi/20) for b in range(-21,20)] for a in range(-5,4)]]
for j in stuff:
new_plot+=j
new_plot.show()