2 files changed, 17 insertions, 0 deletions
diff --git a/3159/CH16/EX16.1/Ex16_1.sce b/3159/CH16/EX16.1/Ex16_1.sce
new file mode 100755
index 000000000..4f35fcf3d
--- /dev/null
+++ b/3159/CH16/EX16.1/Ex16_1.sce
@@ -0,0 +1,13 @@
+// Calculate the net magnetic moment per iron atom in crystal
+clc
+a = 2.87 // lattice parameter in angstrom
+n = 2 // number of atoms per unit cell
+m = 1750 // Saturation magnetization in kAm^-1
+mu = 9.273e-24 // bohr magneton 
+printf("\n Example 16.1")
+m_atom = m*1e3*(a*1e-10)^3 /n
+mu_b = m_atom/mu
+
+printf("\n Net magnetic moment per iron atom in crystal is %.3e Am^2",m_atom)
+printf("\n In unit of mu_b, Net magnetic moment per iron atom in crystal is %.1f ",mu_b)
+
diff --git a/3159/CH16/EX16.1/Ex16_1.txt b/3159/CH16/EX16.1/Ex16_1.txt
new file mode 100755
index 000000000..3f3b86f8f
--- /dev/null
+++ b/3159/CH16/EX16.1/Ex16_1.txt
@@ -0,0 +1,4 @@
+
+ Example 16.1
+ Net magnetic moment per iron atom in crystal is 2.068e-23 Am^2
+ In unit of mu_b, Net magnetic moment per iron atom in crystal is 2.2
+\ No newline at end of file