MoM-based MATLAB program for a charged cube. Write a computer program in MATLAB for...

MoM-based MATLAB program for a charged cube. Write a computer program in MATLAB for the method-of-moments analysis of a charged metallic cube, Fig.1.22, with edge length a = 1 m, and compute the total charge of the cube for V₀ = 1 V and ten, or as many as possible (given available computational resources), subdivisions per cube edge (N = 600 if ten subdivisions per edge are adopted). (ME1- 41.m on IR)

HINT: See the tutorial to MATLAB Exercise 1.39. The following lines of MATLAB code perform meshing of the surface of the cube. We mesh each side of the cube in the same way, with Nside standing for the number of patches on a side (Nside = n²). The local coordinates of centers of patches are calculated only once, and then properly used to compute the global x,y,z coordinates for all sides (on each side, one of the global coordinates is constant), which are merged (for all sides) together.

Once the cube surface is meshed, most lines of the program written in MATLAB Exercise 1.39 can be directly applied. Note that, upon the MoM matrix equation is solved, parts of the column matrix rhos corresponding to individual cube sides should be extracted from it and repacked into a square matrix in order to represent them graphically (as in Fig.1.23). Note also that the obtained charge distributions on all the sides should come out to be the same, due to the symmetry of the problem.

REference: Fig.1.22