blob: 9c1dc0b2ec1dc42c8a564f1dc71082a3b4ed061f (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
// Scilab Code Ex6.1: Page-345 (2008)
clc; clear;
n = 14; // Total number of particles
C = 2; // Total number of compartments
N_micro = C^n; // Total number of microstates
n1 = [10 7 14]; // Set of number of particles in first compartment
n2 = [4 7 0]; // Set of number of particles in second compartment
for i = 1:1:3
W = factorial(n1(i) + n2(i))/(factorial(n1(i))*factorial(n2(i)));
P = W/N_micro;
printf("\nThe probability of microstate (%d, %d) = %8.6f", n1(i), n2(i), P);
end
// Result
// The probability of microstate (10, 4) = 0.061096
// The probability of microstate (7, 7) = 0.209473
// The probability of microstate (14, 0) = 0.000061
|
