1 files changed, 26 insertions, 0 deletions
diff --git a/3689/CH12/EX12.9/12_9.sce b/3689/CH12/EX12.9/12_9.sce
new file mode 100644
index 000000000..4678e367c
--- /dev/null
+++ b/3689/CH12/EX12.9/12_9.sce
@@ -0,0 +1,26 @@
+//////Variable Declaration
+N = [10,50,100]       //Valures for N
+
+//Calculations
+printf("\n   N        ln(N!)       ln(N!)sterling      Error')
+for i =10
+
+ lnN = log(factorial(i))
+    lnNs = i*log(i)-i
+    err = abs(lnN-lnNs)
+    printf('\n%3d       %5.2f        %5.2f              %4.2f',i,lnN,lnNs, err)
+end
+for i =50
+
+ lnN = log(factorial(i))
+    lnNs = i*log(i)-i
+    err = abs(lnN-lnNs)
+    printf('\n%3d       %5.2f        %5.2f              %4.2f',i,lnN,lnNs, err)
+end
+for i =100
+
+ lnN = log(factorial(i))
+    lnNs = i*log(i)-i
+    err = abs(lnN-lnNs)
+    printf('\n%3d       %5.2f        %5.2f              %4.2f',i,lnN,lnNs, err)
+end