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

---
 1332/CH5/EX5.34/5_34.pdf | Bin 0 -> 6418 bytes
 1332/CH5/EX5.34/5_34.sce |  35 +++++++++++++++++++++++++++++++++++
 2 files changed, 35 insertions(+)
 create mode 100755 1332/CH5/EX5.34/5_34.pdf
 create mode 100755 1332/CH5/EX5.34/5_34.sce

(limited to '1332/CH5/EX5.34')

diff --git a/1332/CH5/EX5.34/5_34.pdf b/1332/CH5/EX5.34/5_34.pdf
new file mode 100755
index 000000000..84d494e09
Binary files /dev/null and b/1332/CH5/EX5.34/5_34.pdf differ
diff --git a/1332/CH5/EX5.34/5_34.sce b/1332/CH5/EX5.34/5_34.sce
new file mode 100755
index 000000000..e1ca08e04
--- /dev/null
+++ b/1332/CH5/EX5.34/5_34.sce
@@ -0,0 +1,35 @@
+//Example 5.34
+//QD Method
+//Page no. 202
+clc;clear;close;
+a=[1,2,10,-20]
+for i=1:5
+    e(i,1)=0;
+    e(i,4)=0;
+end
+q(1,1)=-a(2)/a(1);
+q(1,2)=0;q(1,3)=0;
+e(1,2)=a(3)/a(2);
+e(1,3)=a(4)/a(3);
+for i=2:7
+    for j=1:3
+        q(i,j)=e(i-1,j+1)+q(i-1,j)-e(i-1,j)
+    end
+    for j=1:2
+        e(i,j+1)=e(i-1,j+1)*q(i,j+1)/q(i,j)
+    end
+end
+printf('e0\t\tq1\t\te1\t\tq2\t\te2\t\tq3\t\te3\n')
+printf('------------------------------------------------------------------------------------------------------------\n')
+for i=1:7
+    for j=1:3
+        printf('\t\t%.10f\t',q(i,j))
+    end
+    printf('\n')
+    for j=1:4
+        printf('%.10f\t\t\t',e(i,j))
+    end
+    printf('\n')
+end
+printf('\t\t%.10f\t\t\t%.10f\t\t\t%.10f\n',q(7,1),q(7,2),q(7,3))
+printf('\nThe exact roots are \t%.10f    and    %.10f',q(7,1),q(7,3))
\ No newline at end of file
-- 
cgit