# paste this code at the end of VectorFieldPlot 1.0
doc = FieldplotDocument('VFPt_cylindrical_magnets_repelling', commons=True)
field = Field({'coils':[[-1.5, 0, 0, 0.25, 0.75 ,1], [1.5, 0, 0., 0.25, 0.75 ,-1]]})
doc.draw_magnets(field)
n = 15
for i in range(n):
    for x0 in (-1.5, 1.5):
        a = -0.25 + 0.5 * (0.5 + i) / n
        if i in [5, 6, 8, 9]: d = 'backward'
        else: d = 'both'
        if i == n/2: of = [1.0, 0.5, 0.5, 0.2]
        else: of = [1.0, 0.5, 0.5, 1.0]
        line = FieldLine(field, [x0, a], directions=d, maxr=10.0)
        doc.draw_line(line, arrows_style={'dist':1.7, 'offsets':of})
for x0, y0 in ((-0.11, -1), (0.11, -1), (-0.11, 1), (0.11, 1)):
    line = FieldLine(field, [x0, y0], directions='both',
        maxr=10.0, bounds_func=lambda p: abs(p[0])-1.5)
    doc.draw_line(line, arrows_style={'dist':1.7, 'offsets':[0.5, 0.5, 0.8, 1.0]})
doc.write()