Students learn/review how to do integrals in a multivariable context, using the vector differential
\(d\vec{r}=dx\, \hat{x}+dy\, \hat{y}+dz\, \hat{z}\)
and its curvilinear coordinate analogues as a unifying strategy. This strategy is common among physicists, but is NOT typically taught in vector calculus courses and will be new to most students.
Students write python programs to compute the potential due to a square of surface charge, and then to visualize the result. This activity can be used to introduce students to the process of integrating numerically.