diff options
Diffstat (limited to '413/CH4')
-rw-r--r-- | 413/CH4/EX4.1/Table_4_1.sce | 25 | ||||
-rw-r--r-- | 413/CH4/EX4.2/Table_4_2.sce | 16 | ||||
-rw-r--r-- | 413/CH4/EX4.4/Table_4_4.sce | 14 | ||||
-rw-r--r-- | 413/CH4/EX4.5/Table_4_5.sce | 14 |
4 files changed, 69 insertions, 0 deletions
diff --git a/413/CH4/EX4.1/Table_4_1.sce b/413/CH4/EX4.1/Table_4_1.sce new file mode 100644 index 000000000..709546706 --- /dev/null +++ b/413/CH4/EX4.1/Table_4_1.sce @@ -0,0 +1,25 @@ +clc
+clear
+
+function a=eco4(x)
+ a=1.000043+x+x.*x.*0.499219+x.*x.*x/6+x.*x.*x.*x.*0.043750
+endfunction
+function a=eco5(x)
+ a=1.000043+x+x.*x.*0.499219+x.*x.*x/6+x.*x.*x.*x.*0.043750+x.*x.*x.*x.*x/120
+endfunction
+function a=mac4(x)
+ a=1+x+x.*x/2+x.*x.*x/6+x.*x.*x.*x/24
+endfunction
+function a=mac5(x)
+ a=1+x+x.*x/2+x.*x.*x/6+x.*x.*x.*x/24+x.*x.*x.*x.*x/120
+endfunction
+function a=mac6(x)
+ a=1+x+x.*x/2+x.*x.*x/6+x.*x.*x.*x/24+x.*x.*x.*x.*x/120+x.*x.*x.*x.*x.*x/720
+endfunction
+function out=tab4(y)
+out=[y,exp(y),mac6(y),mac5(y),mac4(y),eco5(y),eco4(y)]
+
+endfunction
+for y=0:0.2:1
+disp(tab4(y))
+end
\ No newline at end of file diff --git a/413/CH4/EX4.2/Table_4_2.sce b/413/CH4/EX4.2/Table_4_2.sce new file mode 100644 index 000000000..4a76be593 --- /dev/null +++ b/413/CH4/EX4.2/Table_4_2.sce @@ -0,0 +1,16 @@ +clc
+clear
+
+function a=che(x)
+ a=0.9946+0.9973.*x+x.*x.*0.543+x.*x.*x.*0.1772
+ endfunction
+function a=mac(x)
+ a=1+x+x.*x/2+x.*x.*x/6
+endfunction
+function out=tab(y)
+out=[y,exp(y),che(y),exp(y)-che(y),mac(y),exp(y)-mac(y)]
+
+endfunction
+for y=-1:0.2:1
+disp(tab(y))
+end
\ No newline at end of file diff --git a/413/CH4/EX4.4/Table_4_4.sce b/413/CH4/EX4.4/Table_4_4.sce new file mode 100644 index 000000000..d261a1ed8 --- /dev/null +++ b/413/CH4/EX4.4/Table_4_4.sce @@ -0,0 +1,14 @@ +clc
+clear
+function a=pade(x)
+ a=(x+7*x.*x.*x/9 +64*x.*x.*x.*x.*x/945)/(1+10*x.*x/9 +5*x.*x.*x.*x/21)
+endfunction
+function a=mac(x)
+ a=x-x.*x.*x/3+x.*x.*x.*x.*x/5-x.*x.*x.*x.*x.*x.*x/7+x.*x.*x.*x.*x.*x.*x.*x.*x/9
+endfunction
+function out=tab4(y)
+out=[y,atan(y),pade(y),atan(y)-pade(y),mac(y),atan(y)-mac(y)]
+endfunction
+for y=0.2:0.2:1
+disp(tab4(y))
+end
\ No newline at end of file diff --git a/413/CH4/EX4.5/Table_4_5.sce b/413/CH4/EX4.5/Table_4_5.sce new file mode 100644 index 000000000..e38fbfbc0 --- /dev/null +++ b/413/CH4/EX4.5/Table_4_5.sce @@ -0,0 +1,14 @@ +clc
+clear
+function a=rat(x)
+ a=(1.0018+0.6727*x +0.1598*x.*x)/(1-0.3263*x)
+endfunction
+function a=che(x)
+ a=0.9946+0.9973.*x+x.*x.*0.543+x.*x.*x.*0.1772
+ endfunction
+function out=tab4(y)
+out=[y,exp(y),che(y),exp(y)-che(y), rat(y),exp(y)-rat(y)]
+endfunction
+for y=-1:0.2:1
+disp(tab4(y)')
+end
\ No newline at end of file |