| Differences between revisions 6 and 7 | Back to page |
|
Size: 922
Comment:
|
Size: 1029
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 5: | Line 5: |
| First, the Matlab code. | Here is the Matlab code for Example 2.3. This code is somewhat different from the code in the text book, but does the same thing. |
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.
Here is the Matlab code for Example 2.3. This code is somewhat different from the code in the text book, but does the same thing.
%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