--- /dev/null
+import ElectroMag;
+
+real u = 1cm;
+
+Charge a = pCharge(center=(0.5u,2u));
+Charge b = pCharge(center=(0,0));
+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$");
+Distance dbc = Distance(c.center, b.center, scale=u, L="$r_2$");
+Distance ds[] = {dab};
+Angle abc = Angle(a.center, b.center, c.center, radius=.5u, L="$\theta_T$");
+Angle bac = Angle(b.center, a.center, c.center, radius=-0.4u, L="$\phi_x$");
+Angle as[] = {abc, bac};
+
+//write(a.center.x);
+//write(a.center.y);
+for (int i=0; i<cs.length; i+=1)
+ cs[i].draw();
+for (int i=0; i<ds.length; i+=1)
+ ds[i].draw();
+for (int i=0; i<as.length; i+=1)
+ as[i].draw();
+
+//Vector Fba = CoulombForce(b,a, scale=2u, unit=u);
+//Fba.draw();
+CoulombForces(c=cs, scale=2u, unit=u);
+
+Charge a = pCharge(center=(-2u, 2u));
+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();