Field computation in postprocessing of the MoM solution. Write a function fieldE() i...

Field computation in postprocessing of the MoM solution. Write a function fieldE() in MATLAB that evaluates, in postprocessing, the electric field intensity vector at an arbitrary point in space due to a charged body (e.g., the cube in Fig.1.22), whose charge distribution, approximately described by Eq.(1.55), is determined by the MoM analysis. The input to the program contains EPS0, x0, y0, z0, rhos,S,x,y,z, where EPS0 is the permittivity of air (or another dielectric surrounding the body) and (x0,y0,z0) are Cartesian coordinates of the field point, while rhos,S,x,y,z are arrays containing the surface charge densities and surface areas of MoM patches, and coordinates of their centers. (fieldE.m on IR)

HINT: We approximate the field integral for a surface charge distribution [Fig.1.6 (b)], analogous to the one in Eq.(1.14), in the same way the potential integral in Eq.(1.54) is reduced to its approximate form in Eq.(1.56), based on the charge density approximation in Eq.(1.55), as follows:

and then approximate the integrand by its value at the center of the patch ΔS_i (i = 1, 2, . . . ,N), which results in

MATLAB code based on this equation is similar to that written in MATLAB Exercise 1.36.

Reference: Fig.1.22, Eq. (1.55)

MATLAB