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/CH7/EX7.3/Ex7_3.sce | 52 +++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 52 insertions(+)
 create mode 100755 845/CH7/EX7.3/Ex7_3.sce

(limited to '845/CH7/EX7.3/Ex7_3.sce')

diff --git a/845/CH7/EX7.3/Ex7_3.sce b/845/CH7/EX7.3/Ex7_3.sce
new file mode 100755
index 000000000..1d02509df
--- /dev/null
+++ b/845/CH7/EX7.3/Ex7_3.sce
@@ -0,0 +1,52 @@
+//Example 7.3
+
+clc
+clear
+
+x = 0:4;
+y = [6.9897 7.4036 7.7815 8.1281 8.451];
+
+n = length(x);
+del = %nan*ones(n,5);
+del(:,1) = y';
+for j = 2:6
+    for i = 1:n-j+1
+        del(i,j) = del(i+1,j-1) - del(i,j-1);
+    end
+end
+del(:,1) = [];
+n0 = length(del(1,:));
+
+X = 2;
+i = find(x==X);
+dowy = 0;
+
+for j = 1:n0
+    if j==2*int(j/2) then
+        add = del(i,j);
+    else
+        add = (del(i-1,j) + del(i,j))/2;
+        i = i-1;
+        if i==0 then
+            break
+        end
+    end
+
+    if add == %nan then
+        break
+    else
+        dowy(j) = add;
+    end
+end
+mprintf("%5s %6s %10s %9s %9s %9s",'x','y','dy','d2y','d3y','d4y')
+disp([x' y' del])
+
+mu = 1;
+h = x(2) - x(1);
+dy2 = mu/h*(dowy(1) - 1/6*dowy(3));
+d2y2 = mu/h^2*(dowy(2)-1/12*dowy(4));
+dy2 = round(dy2*10^4)/10^4;
+d2y2 = round(d2y2*10^4)/10^4;
+
+disp(dy2,"y''(2) = ")
+disp(d2y2,"y''''(2) = ")
-- 
cgit