Create a matlab program to determine and plot the energy bands, Ec, Ey and EF in units of eV as a function of position in an n-p abrupt junction diode. Calculate the Fermi level, EF, and plot it on the graph. Assume the parameters indicated: i) Semiconductor: silicon ii) n-side doping: No = 2 x 1015 (cm-3), NA = 0 (cm-3) iii) p-type doping: NA = 1 x 1015 (cm-3), N, = 0 (cm-3) iv) Temperature: T = 300 (K) Permittivity: Esi = 11.9