From 86466678b35deb496f6fd5ad08e527412290b3d3 Mon Sep 17 00:00:00 2001
From: Vasil Zlatanov <v@skozl.com>
Date: Sun, 5 Mar 2017 20:25:30 +0000
Subject: add more initial conditions

---
 coursework17/ex4.m | 45 +++++++++++++++++++++++++++++++++++----------
 1 file changed, 35 insertions(+), 10 deletions(-)

diff --git a/coursework17/ex4.m b/coursework17/ex4.m
index b82b60b..6f44271 100644
--- a/coursework17/ex4.m
+++ b/coursework17/ex4.m
@@ -23,20 +23,45 @@ u = zeros(length(t),length(x));
 u(1,:) = 0;
 u(1,length(x)) = 0;
 
-% Initial condition
-for i = 1:length(x)
-    if x(i) <= 0.5
-        u(1,i) = 2*x(i);;
-    else
-        u(1,i) = 2*(1-x(i)); 
-    end
+% % Initial condition 1
+% for i = 1:length(x)
+%     if x(i) <= 0.5
+%         u(1,i) = 2*x(i);;
+%     else
+%         u(1,i) = 2*(1-x(i)); 
+%     end
+% end
+
+% % Initial condition 2
+% for i = 1:length(x)
+% 	u(1,i) = abs(sin(2*pi*x(i)));
+% %	u(1,i) = sin(2*pi*x(i));
+% end
+
+% % Initial Condition 3 
+%  for i = 0.25/dx:0.75/dx
+%  	u(1,i) = 1;
+%  end
+
+
+% Initial Condition 4
+for i = 1:0.5/dx
+	u(1,i) = -1;
 end
-        
+for i = 0.5/dx+1:length(x)
+	u(1,i) = 1;
+end
+
+
+
+       
 
 for m = 1:length(t)
     % Set boundaries
-    u(m+1,1) = 0.5 * m /length(t); 
-    u(m+1,length(x)) = 0.5 * m /length(t); 
+%    u(m+1,1) = 0.5 * m /length(t); 
+%    u(m+1,length(x)) = 0.5 * m /length(t); 
+    u(m+1,1) = 0;
+    u(m+1,length(x)) = 0;
     for j = 2:(length(x)-1)
         % multiply out (1-2v) and factorise out v
         u(m+1,j) = u(m,j) + ((dt/(dx^2))*(u(m,j+1) - 2*u(m,j) + u(m,j-1))); 
-- 
cgit v1.2.3