initial commit / add all books

1 files changed, 36 insertions, 0 deletions

diff --git a/1670/CH7/EX7.5/7_5.sce b/1670/CH7/EX7.5/7_5.sce
new file mode 100755
index 000000000..e879350d9
--- /dev/null
+++ b/1670/CH7/EX7.5/7_5.sce
@@ -0,0 +1,36 @@
+//Example 7.5

+//Stirlings Central Difference Derivatives

+//Page no. 238

+clc;close;clear;

+printf('   x\t\t    y\t\t    d\t\t    d2\t\t    d3\t\t    d4\n')

+printf('-----------------------------------------------------------------------------------------')

+h=0.1;s=1;

+e=[1,6,30]

+deff('y=f1(x,s)','y=((z(x,3)+z(x-1,3))/2+s*z(x-1,4)+(z(x-1,5)+z(x-2,5))*(3*s^2-1)/12)/h')

+deff('y=f2(x,s)','y=(z(x-1,4))/h^2')

+deff('y=f3(x,s)','y=(z(x-1,5)+z(x-2,5))/(2*h^3)')

+z=[0.7,0.644218;0.8,0.717356;0.9,0.783327;1,0.841471;1.1,0.891207;1.2,0.932039;1.3,0.963558];

+for i=3:6

+    for j=1:9-i

+        z(j,i)=z(j+1,i-1)-z(j,i-1)

+    end

+end

+printf('\n')

+for i=1:7

+    for j=1:6

+        if z(i,j)==0 then

+            printf(' \t')

+        elseif j==1

+            printf('   %.1f\t\t',z(i,j))

+        else

+            printf('%.6f\t',z(i,j))

+        end

+    end

+    printf('\n')

+end

+fp=0;i=5;

+for j=2:2:6

+    fp=fp+((-1)^(j/2+1))*(z(i,j)-z(i-2,j))/(2*h*e(j/2))

+    i=i-1;

+end

+printf('\n\nf`p (sin`(x))= %g',fp)
\ No newline at end of file