Lessons · circuits

The field of a charged rod: where algebra runs out

Regime 2 — calculus does what algebra cannot. The integrand isn't constant, so no single algebra product gives the answer; the accumulated quantity is the area under the curve — and SymPy proves that area is exactly the closed-form result.

Regime 2 · calculus does more Math machine-derived & checked by SymPy 1 modeling assumptions (author-asserted)

A rod of length L=0.3L = 0.3 m with linear charge density λ=100\lambda = 100 nC/m lies along the axis, its near end a=0.1a = 0.1 m from the point where we want the field. Algebra gives the field of a point charge, kq/r2kq/r^2 — but a rod is not a point, and there is no algebra formula for charge spread along a length. Calculus is the only way in: slice the rod into point charges dq=λdxdq = \lambda\,dx, each contributing kdq/x2k\,dq/x^2, and integrate. The shaded area under those contributions is the field E=kλ(1/a1/(a+L))E = k\lambda(1/a - 1/(a+L)), which collapses to the point-charge kQ/a2kQ/a^2 when the rod is short.

Total field E=kQa(a+L)E = \dfrac{kQ}{a(a+L)}
6742 N/C
Point-charge approx kQ/a2kQ/a^2 (overestimate)
26970 N/C
Field from the near half of the rod
5394 N/C
Slice the rod into point charges and add their fields
dE=kdqx2=kλdxx2dE = \frac{k\,dq}{x^2} = \frac{k\lambda\,dx}{x^2}

Cut the rod into slices of length dxdx, each carrying charge dq=λdxdq = \lambda\,dx and acting like a point charge at distance xx. Its field is kdq/x2k\,dq/x^2. Summing over the rod is an integral — the only way to handle a continuous distribution.

check dE/dx=kλ/x2dE/dx = k\lambda/x^2; check E=axkλ/x2dxE = \int_a^x k\lambda/x'^2\,dx'; total kQ/(a(a+L))kQ/(a(a+L)); short-rod limit kQ/a2kQ/a^2

  • The accumulated field's slope is the contribution density: dEdx=kλx2\dfrac{dE}{dx} = \dfrac{k\lambda}{x^2}. [structural]
  • The field is the area under the contributions: E(x)=axkλx2dxE(x) = \int_a^x \dfrac{k\lambda}{x'^2}\,dx'. [structural]
  • The whole rod gives E=kλ(1a1a+L)=kQa(a+L)E = k\lambda\left(\dfrac1a - \dfrac1{a+L}\right) = \dfrac{kQ}{a(a+L)} with Q=λLQ = \lambda L. [structural]
  • For a short rod the field reduces to a point charge: E/Qk/a2E/Q \to k/a^2 as L0L \to 0. [structural]

Dimensional homogeneity: checked by SymPy (holds).

Common misconception: “To get the rod's field, just treat all its charge as sitting at the rod's centre (or near end) and use kQ/r2kQ/r^2.

Lumping the charge at one point gives the wrong answer, because the field is non-linear in distance: the near slices (xx small) contribute far more than the far ones (k/x2k/x^2 falls off steeply), so you cannot replace the spread-out charge with a single point at the average position. Treating the rod as a point at its near end overestimates the field; the true field is the integral E=kλ(1/a1/(a+L))=kQ/(a(a+L))E = k\lambda(1/a - 1/(a+L)) = kQ/(a(a+L)), which is genuinely smaller. Only when the rod is much shorter than its distance (LaL \ll a) do the slices all sit at nearly the same distance and the integral collapses back to the point charge kQ/a2kQ/a^2. The cursor sweeps how much of the rod is included; the area under the contribution curve is the field.

Modeling assumptions — author-asserted, disclosed not discharged
  • A thin rod with uniform linear charge density λ\lambda, and the field evaluated at a point on the rod's axis beyond its near end (so every slice lies along the same line and the contributions add as scalars).

The dE/dx–x graph, fully annotated

A static rendering (Matplotlib): the shaded area under dE/dx is the accumulated integral E, and the slope of E is dE/dx. The interactive version with a draggable cursor is in the Graph tab above.

2026-06-27T05:45:53.007705 image/svg+xml Matplotlib v3.11.0, https://matplotlib.org/