1 files changed, 35 insertions, 0 deletions
diff --git a/845/CH6/EX6.13/Ex6_13.sce b/845/CH6/EX6.13/Ex6_13.sce
new file mode 100755
index 000000000..b7c5b84e5
--- /dev/null
+++ b/845/CH6/EX6.13/Ex6_13.sce
@@ -0,0 +1,35 @@
+//Example 6.13
+
+clc
+clear
+
+x = 1974:2:1982;
+y = [40 43 48 52 57];
+
+n = length(x);
+del = %nan*ones(n,5);
+del(:,1) = y';
+for j = 2:5
+    for i = 1:n-j+1
+        del(i+j-1,j) = del(i+j-1,j-1) - del(i+j-2,j-1);
+    end
+end
+
+X = 1979;
+h = x(2) - x(1);
+p = (X-x(n)) / h;
+xn = x(n);
+yn = y(n);
+delyn = del(n,:);
+
+Y = 0;
+
+for i = 0:length(delyn)-1
+    t = 1;
+    for j = 1:i
+        t = t * (p+j-1);
+    end
+    Y = Y + t*delyn(i+1)/factorial(i);
+end
+Y = round(Y*10^4)/10^4;
+disp(Y,"Estimated sales for the year 1979: ")