Students calculate electrostatic fields (\(V\), \(\vec{E}\)) and magnetostatic fields (\(\vec{A}\), \(\vec{B}\)) from charge and current sources with a common geometry. The sequence of activities is arranged so that the mathematical complexity of the formulas students encounter increases with each activity. Several auxiliary activities allow students to focus on the geometric/physical meaning of the distance formula, charge densities, and steady currents. A meta goal of the entire sequence is that students gain confidence in their ability to parse and manipulate complicated equations.

A short improvisational role-playing skit based on the *Star Trek* series in which students explore the definition and notation for position vectors, the importance of choosing an origin, and the geometric nature of the distance formula.
\[\vert\vec{r}-\vec{r}^\prime\vert=\sqrt{(x-x^\prime)^2+(y-y^\prime)^2-(z-z^\prime)^2}\]

Students work in small groups to use the superposition principle
\[V(\vec{r}) = \frac{1}{4\pi\epsilon_0}\sum_i \frac{q_i}{\vert\vec{r}-\vec{r}_i\vert}\]
to find the electrostatic potential \(V\) everywhere in space due to a pair of charges (either identical charges or a dipole).
Different groups are assigned different arrangements of charges and different regions of space to consider: either on the axis of the charges or in the plane equidistant from the two charges, for either small or large values of the relevant geometric variable.
Each group is asked to find a power series expansion for the electrostatic potential, valid in their group's assigned region of space.
The whole class wrap-up discussion then compares and contrasts the results and discuss the symmetries of the two cases.

The distance \(\left\vert\vec r -\vec r\,{}'\right\vert\) between the point \(\vec r\) and the point \(\vec r'\) is a coordinate-independent, physical and geometric quantity. But, in practice, you will need to know how to express this quantity in different coordinate systems.

- Find the distance \(\left\vert\vec r -\vec r\,{}'\right\vert\) between the point \(\vec r\) and the point \(\vec r'\) in rectangular coordinates.
- Show that this same distance written in cylindrical coordinates is: \begin{equation} \left|\vec r -\vec r\,{}'\right| =\sqrt{s^2+s\,{}'^2-2ss\,{}'\cos(\phi-\phi\,{}') +(z-z\,{}')^2} \end{equation}
- Show that this same distance written in spherical coordinates is: \begin{equation} \left\vert\vec r -\vec r\,{}'\right\vert =\sqrt{r'^2+r\,{}^2-2rr\,{}' \left[\sin\theta\sin\theta\,{}'\cos(\phi-\phi\,{}') +\cos\theta\cos\theta\,{}'\right]} \end{equation}
- Now assume that \(\vec r\,{}'\) and \(\vec r\) are in the \(x\)-\(y\) plane. Simplify the previous two formulas.

Students, pretending they are point charges, move around the room acting out various prompts from the instructor regarding charge densities, including linear \(\lambda\), surface \(\sigma\), and volume \(\rho\) charge densities, both uniform and non-uniform. The instructor demonstrates what it means to measure these quantities. In a remote setting, we have students manipulate 10 coins to model the prompts in this activity and the we demonstrate the answers with coins under a doc cam.

Students work in small groups to use the superposition principle \[V(\vec{r}) =\frac{1}{4\pi\epsilon_0}\int\frac{\rho(\vec{r}^{\,\prime})}{\vert \vec{r}-\vec{r}^{\,\prime}\vert} \, d\tau^{\prime}\] to find an integral expression for the electrostatic potential, \(V(\vec{r})\), everywhere in space, due to a ring of charge.

In an optional extension, students find a series expansion for \(V(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.

Students work in small groups to use Coulomb's Law \[\vec{E}(\vec{r}) =\frac{1}{4\pi\epsilon_0}\int\frac{\rho(\vec{r}^{\,\prime})\left(\vec{r}-\vec{r}^{\,\prime}\right)}{\vert \vec{r}-\vec{r}^{\,\prime}\vert^3} \, d\tau^{\prime}\] to find an integral expression for the electric field, \(\vec{E}(\vec{r})\), everywhere in space, due to a ring of charge.

In an optional extension, students find a series expansion for \(\vec{E}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.

Students, pretending they are point charges, move around the room so as to make an imaginary magnetic field meter register a constant magnetic field, introducing the concept of *steady current*. Students act out linear \(\vec{I}\), surface \(\vec{K}\), and volume \(\vec{J}\) current densities. The instructor demonstrates what it means to measure these quantities by counting how many students pass through a gate.

Students work in small groups to use the superposition principle \[\vec{A}(\vec{r}) =\frac{\mu_0}{4\pi}\int\frac{\vec{J}(\vec{r}^{\,\prime})}{\vert \vec{r}-\vec{r}^{\,\prime}\vert}\, d\tau^{\prime}\] to find an integral expression for the magnetic vector potential, \(\vec{A}(\vec{r})\), due to a spinning ring of charge.

In an optional extension, students find a series expansion for \(\vec{A}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.

Students work in small groups to use the Biot-Savart law \[\vec{B}(\vec{r}) =\frac{\mu_0}{4\pi}\int\frac{\vec{J}(\vec{r}^{\,\prime})\times \left(\vec{r}-\vec{r}^{\,\prime}\right)}{\vert \vec{r}-\vec{r}^{\,\prime}\vert^3} \, d\tau^{\prime}\] to find an integral expression for the magnetic field, \(\vec{B}(\vec{r})\), due to a spinning ring of charge.

In an optional extension, students find a series expansion for \(\vec{B}(\vec{r})\) either on the axis or in the plane of the ring, for either small or large values of the relevant geometric variable. Add an extra half hour or more to the time estimate for the optional extension.