; see power_solid_angle

(eq E
    (/ I (exp r 2)))