diff --git a/membrane_solution31.m b/membrane_solution.m
similarity index 68%
rename from membrane_solution31.m
rename to membrane_solution.m
index 1e1aca6..d0b1734 100644
--- a/membrane_solution31.m
+++ b/membrane_solution.m
@@ -1,4 +1,8 @@
-function [w] = membrane_solution31(T,P,n)% Set up initial matrix
+function [w] = membrane_solution(T,P,n)% Set up initial matrix
+    % T = given tension (microNewton/micrometer)
+    % P = given pressure (MPa)
+    % n = number of interior nodes
+    
     od = ones(n^2-1,1);
     od(n:n:end) = 0;
     k = -4*diag(ones(n^2,1))+diag(ones((n^2)-n,1),n)+diag(ones((n^2)-n,1),-n)+diag(od,1)+diag(od,-1);
diff --git a/membrane_solution3.m b/membrane_solution3.m
index 39211e2..ec4ca80 100644
--- a/membrane_solution3.m
+++ b/membrane_solution3.m
@@ -1,5 +1,7 @@
-function [w] = membrane_solution3(T,P)
-    % Set up initial matrix
+function [w] = membrane_solution3(T,P) % Set up initial matrix
+    % T = given tension (microNewton/micrometer)
+    % P = given pressure (MPa)
+    
     od = ones(8,1);
     od(3:3:end) = 0;
     k = -4*diag(ones(9,1))+diag(ones(9-3,1),3)+diag(ones(9-3,1),-3)+diag(od,1)+diag(od,-1);
diff --git a/membrane_solution3_script.m b/membrane_solution3_script.m
new file mode 100644
index 0000000..9a18ae1
--- /dev/null
+++ b/membrane_solution3_script.m
@@ -0,0 +1 @@
+[w] = membrane_solution3(0.006,0.001);
diff --git a/membrane_solution_script.m b/membrane_solution_script.m
new file mode 100644
index 0000000..15fc1b5
--- /dev/null
+++ b/membrane_solution_script.m
@@ -0,0 +1 @@
+[w] = membrane_solution(0.006,0.001,10)
\ No newline at end of file
diff --git a/script1.m b/script1.m
deleted file mode 100644
index d675de7..0000000
--- a/script1.m
+++ /dev/null
@@ -1,15 +0,0 @@
-% Set up initial matrix
-od = ones(8,1);
-od(3:3:end) = 0;
-k = -4*diag(ones(9,1))+diag(ones(9-3,1),3)+diag(ones(9-3,1),-3)+diag(od,1)+diag(od,-1);
-
-% Solve for unknown matrix, w
-P = 0.001; %MPa
-T = 0.006; %microN/microm
-y = -(10/4)^2*(P/T)*ones(9,1); %output vector 
-w = k\y; %solves for w in microm
-
-[x,y] = meshgrid(0:10/4:10,0:10/4:10);
-z = zeros(size(x));
-z(2:end-1,2:end-1) = reshape(w,[3 3]); 
-surf(x,y,z)
\ No newline at end of file
diff --git a/script2.asv b/script2.asv
deleted file mode 100644
index a1e15b6..0000000
--- a/script2.asv
+++ /dev/null
@@ -1,15 +0,0 @@
-% Set up initial matrix
-od = ones(n^2-1,1);
-od(n:n:end) = 0;
-k = -4*diag(ones(n^2,1))+diag(ones((n^2)-n,1),n)+diag(ones((n^2)-n,1),-n)+diag(od,1)+diag(od,-1);
-
-% Solve for unknown matrix, w
-P = 0.001; %MPa
-T = 0.006; %microN/microm
-y = -(10/(n+1))^2*(P/T)*ones(n^2,1); %output vector 
-w = k\y; %solves for w in microm
-
-[x,y] = meshgrid(0:10/4:10,0:10/4:10);
-z = zeros(size(x));
-z(2:end-1,2:end-1) = reshape(w,[3 3]); 
-surf(x,y,z)
\ No newline at end of file
diff --git a/script2.m b/script2.m
deleted file mode 100644
index fb841e8..0000000
--- a/script2.m
+++ /dev/null
@@ -1,15 +0,0 @@
-function [w] = membrane_solution31(T,P,n);% Set up initial matrix
-od = ones(n^2-1,1);
-od(n:n:end) = 0;
-k = -4*diag(ones(n^2,1))+diag(ones((n^2)-n,1),n)+diag(ones((n^2)-n,1),-n)+diag(od,1)+diag(od,-1);
-
-% Solve for unknown matrix, w
-P = 0.001; %MPa
-T = 0.006; %microN/microm
-y = -(10/(n+1))^2*(P/T)*ones(n^2,1); %output vector 
-w = k\y; %solves for w in microm
-
-[x,y] = meshgrid(0:10/(n+1):10,0:10/(n+1):10);
-z = zeros(size(x));
-z(2:end-1,2:end-1) = reshape(w,[n n]); 
-surf(x,y,z)
\ No newline at end of file