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

---
 1673/CH1/EX1.13/1_13.sce | 25 +++++++++++++++++++++++++
 1 file changed, 25 insertions(+)
 create mode 100755 1673/CH1/EX1.13/1_13.sce

(limited to '1673/CH1/EX1.13')

diff --git a/1673/CH1/EX1.13/1_13.sce b/1673/CH1/EX1.13/1_13.sce
new file mode 100755
index 000000000..61aea03e1
--- /dev/null
+++ b/1673/CH1/EX1.13/1_13.sce
@@ -0,0 +1,25 @@
+//taylor series
+//example 1.13
+//page 16
+clc;clear;close;
+deff('y=f(x)','y=sin(x)');
+deff('y=f1(x)','y=cos(x)');
+deff('y=f2(x)','y=-sin(x)');
+deff('y=f3(x)','y=-cos(x)');
+deff('y=f4(x)','y=sin(x)');
+deff('y=f5(x)','y=cos(x)');
+deff('y=f6(x)','y=-sin(x)');
+deff('y=f7(x)','y=-cos(x)');
+D=[f(%pi/6) f1(%pi/6) f2(%pi/6) f3(%pi/6) f4(%pi/6) f5(%pi/6) f6(%pi/6) f7(%pi/6)];
+S1=0;
+h=%pi/6;
+printf('order of approximation  computed value of sin(pi/3)   absolute eror\n\n');
+for j=1:8
+for i=1:j
+    S1=S1+h^(i-1)*D(i)/factorial(i-1);
+end
+printf('%d                       %0.9f                        %0.9f\n',j,S1,abs(sin(%pi/3)-S1));
+S1=0;
+end
+
+
-- 
cgit