Please show all MATLAB script for this problem. Thank you.
10. Use the MATLAB command "taylor" to compute Taylor expansions of the Fermi-Dirac function:
where k is Boltzmann's constant, T is temperature, and EF is the Fermi energy. Take EF = 0 and kT to be 0.026 eV (its value at T=300 K). Taylor expand this function about E = 0 to different orders (from 1 to say, 10 terms). Plot the approximation of f(E) for different orders of the Taylor expansion. Overlap this plot on top of the Fermi-Dirac function. Show that as more terms are used in the Taylor expansion, the more accurate is the approximation to the Fermi function by Taylor expansion. (Include plot and MATLAB code with your homework.)