1 files changed, 14 insertions, 0 deletions
diff --git a/1736/CH4/EX4.22/Ch04Ex22.sce b/1736/CH4/EX4.22/Ch04Ex22.sce
new file mode 100755
index 000000000..2afa025d8
--- /dev/null
+++ b/1736/CH4/EX4.22/Ch04Ex22.sce
@@ -0,0 +1,14 @@
+// Scilab Code Ex4.22: Page-137 (2006)
+clc; clear;
+rs_a0_ratio = 2.07;   // Solid radius to lattice parameter ratio for Al
+E_F = 50.1*(rs_a0_ratio)^(-2);  // Fermi level energy of Y, eV
+// According to Jellium model, h_cross*omega_P = E = 47.1 eV *(rs_a0_ratio)^(-3/2)
+E = 47.1*(rs_a0_ratio)^(-3/2);      // Plasmon energy of Al, eV
+printf("\nThe plasmon energy of Al = %4.2f eV", E);
+printf("\nThe experimental value is 15 eV");
+
+// Result 
+// The plasmon energy of Al = 15.81 eV
+// The experimental value is 15 eV 
+
+