Electric field with material present.
28.2 is the response in different directions and higher order dependencies.
When they introduce dependencies of arguments in the linear term (28.3), why is it done via convolution? Why not just write out E(r,r',..) etc.?
Are they assuming a different epsilon?