Diff for "Homework 4"

Differences between revisions 1 and 5 (spanning 4 versions)

Back to page

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 }}}

-  ⇤ ← Revision 1 as of 3:06AM, Apr 25, 2011 → 
  Size: 25
  Editor: Gweon
  Comment:
+   ← Revision 5 as of 3:47AM, Apr 25, 2011 → ⇥
  Size: 921
  Editor: Gweon
  Comment:
-Deletions are marked like this.
+Additions are marked like this.
 Line 1:
-[[ Matlab and Python ]]
+[[ Matlab and Python | 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
}}}