From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 845/CH6/EX6.22/Ex6_22.sce | 42 ++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100755 845/CH6/EX6.22/Ex6_22.sce

(limited to '845/CH6/EX6.22/Ex6_22.sce')

diff --git a/845/CH6/EX6.22/Ex6_22.sce b/845/CH6/EX6.22/Ex6_22.sce
new file mode 100755
index 000000000..e1bf6e2bb
--- /dev/null
+++ b/845/CH6/EX6.22/Ex6_22.sce
@@ -0,0 +1,42 @@
+//Example 6.22
+
+clc
+clear
+
+x = 3:8;
+y = [0.205 0.24 0.259 0.262 0.25 0.224];
+
+n = length(x);
+del = %nan*ones(n,5);
+del(:,1) = y';
+for j = 2:5
+    for i = 1:n-j+1
+        del(i,j) = del(i+1,j-1) - del(i,j-1);
+    end
+end
+
+X = poly(0, "X");
+x0 = x(1);
+y0 = y(1);
+h = x(2) - x(1);
+p = (X-x0)/h;
+del0 = del(1,:);
+del0 = round(del0*10^4)/10^4;
+del0 = del0(1:find(del0==0)-1);
+
+Y = 0;
+for i = 0:length(del0)-1
+    t = 1;
+    for j = 1:i
+        t = t * (p-j+1);
+    end
+    Y = Y + t*del0(i+1)/factorial(i);
+end
+disp(Y,"y = ")
+
+dydx = derivat(Y);
+minx = roots(dydx);
+miny = round(horner(Y,minx)*10^5)/10^5;
+disp(minx,"min_x = ")
+disp(miny,"min_y = ")
+//min_y value is incorrectly displayed in textbook as 0.25425 instead of 0.26278
-- 
cgit