a1(x, y) = x - y - ((x^(2) + y^(2)) * x)
a2(x, y) = x + y - ((x^(2) + y^(2)) * y)
a = CurveCartesian(cos(t), sin(t), t, 0, (2 * π))
t1 = Sequence(((2 * π) * k / 10), k, 0, 9)
r: {0.5, 1, 1.5}
p1 = Sequence((r(1) * a(t1(i))), i, 1, Length[t1])
p3 = Sequence((r(3) * a(t1(i))), i, 1, Length[t1])
p2 = Sequence((r(2) * a(t1(i))), i, 1, Length[t1])
p = Join({p1, p2, p3})
vec1 = Sequence(Vector[p1(i), p1(i) + (0.5 * (a1(p1(i)), a2(p1(i))))], i, 1, Length[p1])
vec3 = Sequence(Vector[p3(i), p3(i) + (0.2 * (a1(p3(i)), a2(p3(i))))], i, 1, Length[p3])
vec2 = Sequence(Vector[p2(i), p2(i) + (0.35 * (a1(p2(i)), a2(p2(i))))], i, 1, Length[p2])
vec = Join({vec1, vec2, vec3})