First, here is some base-line discussion for Matlab and Python.

Now, for this homework, you basically need to calculate some functions and plot up the results. Example 2.3 of the textbook is a pretty good "template" for this task, and so I will discuss that example.

First, the Matlab code.

{{{!highlight matlab %E2.3: Fermi Function Calculation, f(E-EF,T)

%Initialization clear close

%Constant %25.85 meV for 300 K is an equivalent way to remember it. k=8.617e-5;

%Google for "Matlab linspace", to find out what linspace does! dE=linspace(-0.2,0.2); for ii=1:4;

• T=100*ii; kT=k*T; f(ii,:)=1./(1+exp(dE./kT)); end

%Plotting result close plot(dE,f); grid; xlabel('E - E_F (eV)'); ylabel('f (E)'); text(.05,.22,'T=400K'); text(-.03,.12,'T=100K');

%Octave-specific -- uncomment this line for using with Octave. %print -djpg E2.3.oct.jpg }}}