Adding plots

17 files changed, 10 insertions, 7 deletions

diff --git a/Part 3/figures/abs_sin_150_5000_2.png b/Part 3/figures/abs_sin_150_5000_2.png
new file mode 100644
index 0000000..9688ed0
--- /dev/null
+++ b/Part 3/figures/abs_sin_150_5000_2.png
diff --git a/Part 3/figures/arcsinh3D_150_5000.png b/Part 3/figures/arcsinh3D_150_5000.png
new file mode 100644
index 0000000..8ec71b8
--- /dev/null
+++ b/Part 3/figures/arcsinh3D_150_5000.png
diff --git a/Part 3/figures/arcsinh3D_150_5000_bc12_bc0.png b/Part 3/figures/arcsinh3D_150_5000_bc12_bc0.png
new file mode 100644
index 0000000..8d752be
--- /dev/null
+++ b/Part 3/figures/arcsinh3D_150_5000_bc12_bc0.png
diff --git a/Part 3/figures/arcsinh3D_150_5000_bcfunc.png b/Part 3/figures/arcsinh3D_150_5000_bcfunc.png
new file mode 100644
index 0000000..a2fdc8b
--- /dev/null
+++ b/Part 3/figures/arcsinh3D_150_5000_bcfunc.png
diff --git a/Part 3/figures/arcsinh_150_5000.png b/Part 3/figures/arcsinh_150_5000.png
new file mode 100644
index 0000000..dd5c390
--- /dev/null
+++ b/Part 3/figures/arcsinh_150_5000.png
diff --git a/Part 3/figures/arcsinh_150_5000_bc12_bc0.png b/Part 3/figures/arcsinh_150_5000_bc12_bc0.png
new file mode 100644
index 0000000..9ea35a0
--- /dev/null
+++ b/Part 3/figures/arcsinh_150_5000_bc12_bc0.png
diff --git a/Part 3/figures/arcsinh_150_5000_bcfunc.png b/Part 3/figures/arcsinh_150_5000_bcfunc.png
new file mode 100644
index 0000000..b7f6260
--- /dev/null
+++ b/Part 3/figures/arcsinh_150_5000_bcfunc.png
diff --git a/Part 3/figures/cos.png b/Part 3/figures/cos.png
new file mode 100644
index 0000000..5d44379
--- /dev/null
+++ b/Part 3/figures/cos.png
diff --git a/Part 3/figures/cos3D.png b/Part 3/figures/cos3D.png
new file mode 100644
index 0000000..f1810bd
--- /dev/null
+++ b/Part 3/figures/cos3D.png
diff --git a/Part 3/figures/cos3D_150_5000.png b/Part 3/figures/cos3D_150_5000.png
new file mode 100644
index 0000000..1d6a55e
--- /dev/null
+++ b/Part 3/figures/cos3D_150_5000.png
diff --git a/Part 3/figures/cos_150_5000.png b/Part 3/figures/cos_150_5000.png
new file mode 100644
index 0000000..34905bc
--- /dev/null
+++ b/Part 3/figures/cos_150_5000.png
diff --git a/Part 3/figures/exp3D_150_5000.png b/Part 3/figures/exp3D_150_5000.png
new file mode 100644
index 0000000..d522df0
--- /dev/null
+++ b/Part 3/figures/exp3D_150_5000.png
diff --git a/Part 3/figures/exp_150_5000.png b/Part 3/figures/exp_150_5000.png
new file mode 100644
index 0000000..f224153
--- /dev/null
+++ b/Part 3/figures/exp_150_5000.png
diff --git a/Part 3/figures/sin_150_5000_2.png b/Part 3/figures/sin_150_5000_2.png
new file mode 100644
index 0000000..a08b44c
--- /dev/null
+++ b/Part 3/figures/sin_150_5000_2.png
diff --git a/Part 3/figures/tirang_150_7500.png b/Part 3/figures/tirang_150_7500.png
new file mode 100644
index 0000000..5712fd6
--- /dev/null
+++ b/Part 3/figures/tirang_150_7500.png
diff --git a/Part 3/scripts/finite_script.m b/Part 3/scripts/finite_script.m
index a6bdb5d..96ea4ce 100644
--- a/Part 3/scripts/finite_script.m
+++ b/Part 3/scripts/finite_script.m
@@ -3,13 +3,13 @@
 
 clear;
 
-steps = 100;
+steps = 300;
 
 % Set the number of samples to take
 N = 150;
 
 % Set the total time to run
-m = 5000;
+m = 3000;
 
 % Declare final size of matrix for speed
 res = zeros(N+1, m+1);
@@ -22,11 +22,14 @@ h = 1/N;
 k = h^2 * v;
 
 % Set Initial condition
-res(:, 1) = get_function(N, h, 3)';
+res(:, 1) = get_function(N, h, 2)';
 
 % Set boundary conditions
-res(1, :) = zeros(1, m+1);
-res(N+1, :) = zeros(1, m+1);
+%res(1, :) = zeros(1, m+1);
+%res(N+1, :) = zeros(1, m+1);
+
+res(1, :) = sin(2*pi*(0:1/m:1));
+res(N+1, :) = sin(2*pi*(0:1/m:1));
 
 % Calculate M+1 and plot it continuously
 for c = 1:m
@@ -52,7 +55,7 @@ hold off;
 xlabel('x');
 ylabel('y');
 title('Plots of 1D Heat equation over time, bc = 1');
-legend('m = 0', 'm = 100', 'm = 200', 'm = ...', 'm = 5000');
+legend('m = 0', 'm = 500', 'm = 1000', 'm = ...', 'm = 5000');
 
 % Plot 3D
 figure;
diff --git a/Part 3/scripts/get_function.m b/Part 3/scripts/get_function.m
index dc6c823..9f357da 100644
--- a/Part 3/scripts/get_function.m
+++ b/Part 3/scripts/get_function.m
@@ -13,7 +13,7 @@ if f == 1
         count = count + 1;
     end
 elseif f == 2
-    y = sin(2*pi*(0:h:1));
+    y = cos(2*pi*(0:h:1));
 elseif f == 3
     y = abs(sin(2*pi*(0:h:1)));
 elseif f == 4