real u = 1cm;
Charge a = pCharge(center=(0.5u,2u));
-Charge b = pCharge(center=(0,0));
+Charge b = aCharge(center=(0,0), q=1);
Charge c = nCharge(center=(-3u,0), q=-3);
Charge cs[] = {a,b,c};
Distance dab = Distance(b.center, a.center, scale=u, L="$r_1$");
//Fba.draw();
CoulombForces(c=cs, scale=2u, unit=u);
-Charge a = pCharge(center=(-2u, 2u));
+Charge a = aCharge(center=(-2u, 2u), q=0);
a.draw();
Vector v = EField(a.center, mag=u/2, dir=0,L="E"); v.draw();
Vector v = BField(a.center, mag=u/2, dir=90,L="B"); v.draw();
Vector v = Velocity(a.center, mag=u/2, dir=180,L="$v$"); v.draw();
+
+a = aCharge(center=(-3u, 2u), q=-1);
+a.draw();