summaryrefslogtreecommitdiff
path: root/1736/CH8/EX8.9/Ch08Ex9.sce
blob: 9a268bc8f8e1bd6bf89c3d45cdafb1d6f1e05bc1 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
// Scilab code Ex8.9 Page:250 (2006)
clc; clear;
e = 1.6e-019;   // Energy equivalent of 1 eV, J/eV
mu = 9.29e-024;    // Bohr magneton, J/T
mu_0 = 1.26e-006;    // Permeability of free space, Sq. tesla cubic meter per joule
E_F= 11.63*e;    // Fermi energy, J
N = 6.02e+028;    // Atomic concentration, atoms per cubic meter 
chi_Total = 2.2e-005;   // Paramagnetic susceptibility of Mg, S.I. units
chi_Pauli = 3*N*mu^2*mu_0/(2*E_F);    // Pauli diamagnetism, S.I. units
chi_dia = chi_Total - chi_Pauli;    // Diamagnetic contribution to magnetic susceptibility

printf("\nThe Pauli spin susceptibility of Al = %5.3e S.I. units", chi_Pauli);
printf("\nThe diamagnetic contribution to magnetic susceptibility of Al = %5.3e S.I. units", chi_dia);

// Result 
// The Pauli spin susceptibility of Al = 5.277e-06 S.I. units
// The diamagnetic contribution to magnetic susceptibility of Al = 1.672e-05 S.I. units