diff options
Diffstat (limited to '1332/CH5/EX5.34')
-rwxr-xr-x | 1332/CH5/EX5.34/5_34.pdf | bin | 0 -> 6418 bytes | |||
-rwxr-xr-x | 1332/CH5/EX5.34/5_34.sce | 35 |
2 files changed, 35 insertions, 0 deletions
diff --git a/1332/CH5/EX5.34/5_34.pdf b/1332/CH5/EX5.34/5_34.pdf Binary files differnew file mode 100755 index 000000000..84d494e09 --- /dev/null +++ b/1332/CH5/EX5.34/5_34.pdf 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 |