summaryrefslogtreecommitdiff
path: root/413/CH5/EX5.8/Table_5_3.sce
diff options
context:
space:
mode:
Diffstat (limited to '413/CH5/EX5.8/Table_5_3.sce')
-rw-r--r--413/CH5/EX5.8/Table_5_3.sce42
1 files changed, 42 insertions, 0 deletions
diff --git a/413/CH5/EX5.8/Table_5_3.sce b/413/CH5/EX5.8/Table_5_3.sce
new file mode 100644
index 000000000..5964235bf
--- /dev/null
+++ b/413/CH5/EX5.8/Table_5_3.sce
@@ -0,0 +1,42 @@
+clc
+clear
+clc
+clear
+x=[1.7 1.8 2 2.35 2.50];
+function a=expsin(x)
+ a=exp(x).*sin(x)
+endfunction
+for i=1:1:4
+function b=firstdivdiff(x)
+b=(expsin(x(1,i+1))-expsin(x(1,i)))/(x(1,i+1)-x(1,i))
+endfunction
+A(1,i)=firstdivdiff(x)
+end
+for i=1:1:3
+function b=Secdivdiff(x)
+b=(A(1,i+1)-A(1,i))/(x(1,i+2)-x(1,i))
+endfunction
+B(1,i)=Secdivdiff(x)
+end
+for i=1:1:2
+function b=thdivdiff(x)
+b=(B(1,i+1)-B(1,i))/(x(1,i+3)-x(1,i))
+endfunction
+C(1,i)=thdivdiff(x)
+end
+for i=1:1:1
+function b=fodivdiff(x)
+b=(C(1,i+1)-C(1,i))/(x(1,i+4)-x(1,i))
+endfunction
+D(1,i)=fodivdiff(x)
+end
+out=[0,x(1,1),expsin(x(1,1)) ]
+disp(out)
+out1=[1,x(1,2),expsin(x(1,2)), A(1,1) ]
+disp(out1)
+out2=[2,x(1,3),expsin(x(1,3)), A(1,2), B(1,1) ]
+disp(out2)
+out3=[3,x(1,4),expsin(x(1,4)), A(1,3), B(1,2),C(1,1),D(1,1) ]
+disp(out3)
+out4=[4,x(1,5),expsin(x(1,5)), A(1,4), B(1,3),C(1,2) ]
+disp(out4) \ No newline at end of file