initial commit / add all books

3 files changed, 40 insertions, 0 deletions

diff --git a/964/CH26/EX26.1/26_1.sce b/964/CH26/EX26.1/26_1.sce
new file mode 100755
index 000000000..0af6c6d95
--- /dev/null
+++ b/964/CH26/EX26.1/26_1.sce
@@ -0,0 +1,40 @@
+clc;

+clear;

+function yp=f(t,y)

+    yp=-1000*y+3000-2000*exp(-t)

+endfunction

+y0=0;

+//explicit euler

+h1=0.0005;

+y1(1)=y0;

+count=2;

+t=0:0.0001:0.006

+for i=0:0.0001:0.0059

+    y1(count)=y1(count-1)+f(i,y1(count-1))*h1

+    count=count+1;

+end

+plot(t,y1)

+h2=0.0015;

+y2(1)=y0;

+count=2;

+t=0:0.0001:0.006

+for i=0:0.0001:0.0059

+    y2(count)=y2(count-1)+f(i,y2(count-1))*h2

+    count=count+1;

+end

+plot(t,y2,'.-')

+xtitle("y vs t","t","y")

+h=legend(["h-0.0005","h=0.0015"])

+clf();

+//implicit order

+h3=0.05;

+y3(1)=y0;

+count=2;

+t=0:0.01:0.4

+for j=0:0.01:0.39

+    y3(count)=(y3(count-1)+3000*h3-2000*h3*exp(-(j+0.01)))/(1+1000*h3)

+    count=count+1;

+end

+

+plot(t,y3)

+xtitle("y vs t","t","y")

diff --git a/964/CH26/EX26.1/26_1_1.jpeg b/964/CH26/EX26.1/26_1_1.jpeg
new file mode 100755
index 000000000..29a2401f8
--- /dev/null
+++ b/964/CH26/EX26.1/26_1_1.jpeg
diff --git a/964/CH26/EX26.1/26_1_2.jpeg b/964/CH26/EX26.1/26_1_2.jpeg
new file mode 100755
index 000000000..33d1479ac
--- /dev/null
+++ b/964/CH26/EX26.1/26_1_2.jpeg