# Notes of interest for the course

Official version of Course
description .

A
note on math for thermodynamics (from M. Ediger)

Golden Sheet Autumn 2011
version 0.1

#### Old Exam questions used in the course

9 ects course was on the masters level, 7.5 ects is a bachelor
course.
Difference between "English" and "Englisk (translated)" is that the
former is real English sets from exams given in English, where the
later is my (sometimes rather quick) translation.

The collection of old exam questions (from 1976 - 2004) in Danish: Text from
IMFUFA 429 contains a few more old problems

#### In Danish

Et sæt noter skrevet af en studerende på et tidligere
Termo-kursus kan findes her

En matlab-note findes her

Videnskabens verden (DR P1 program) har et sjovt inslag om entropi 22
minutter inde i deres
podcast fra den 26-07-2011.

## Other simulation stuff

This page has
lots of fun to play with. Try e.g. this.

## Some matlab stuff to play with

### Two state paramagnet

Matlab/octave program for
calculating U,S,T,M,C for Two state paramagnet (problem 3.17 and 3.18).
### Multiplicity of two weekly coupled "Einstein solids"

A matlab/octave program for
calculating the multiplicity of two weekly coupled "Einstein solids"
(problem 2.9). It can be used as::

Save the program as TwoEinsteinSolids.m

> data = TwoEinsteinSolids(200,100,100);

> plot(data(:,1),data(:,6))

Alternatively one can use the function nchoosek
which directly calculates the Binomial coefficient for large
numbers.

### Random Walk for a ensemble of particles

A matlab/octave program for simulating Random Walk
for a ensemble of particles.

Example of use::

> x0 = zeros(1000000,1); % Start a million
particles at x=0

> x = RW(x0,100);
% Let each particle do 100 random jumps

> [yy,xx] = hist(x,-100:2:100); % Plot a histogram

> plot(xx', yy'/1e6, 'o', xx', sqrt(2/pi/100)*exp(-xx.^2/2/100)) %
Plot result and analytical solution

If you want to see how the particles are moving :

> x = cumsum( sign(rand(1000,100)-0.5)); % 100 particles, 1000
steps eachr

> plot(x(:,1:10)) % Plot 10 first particles.