From fe2d8b507a4fa4a45984093049df4112b6e47f23 Mon Sep 17 00:00:00 2001
From: Leahy <cpl11001@ad.engr.uconn.edu>
Date: Thu, 7 Dec 2017 12:19:12 -0500
Subject: [PATCH] 1

---
 membrane_solution3.m  | 15 +++++++++++++++
 membrane_solution31.m | 14 ++++++++++++++
 script1.m             | 15 +++++++++++++++
 script2.asv           | 15 +++++++++++++++
 script2.m             | 15 +++++++++++++++
 5 files changed, 74 insertions(+)
 create mode 100644 membrane_solution3.m
 create mode 100644 membrane_solution31.m
 create mode 100644 script1.m
 create mode 100644 script2.asv
 create mode 100644 script2.m

diff --git a/membrane_solution3.m b/membrane_solution3.m
new file mode 100644
index 0000000..39211e2
--- /dev/null
+++ b/membrane_solution3.m
@@ -0,0 +1,15 @@
+function [w] = membrane_solution3(T,P)
+    % 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
+    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)
+end
\ No newline at end of file
diff --git a/membrane_solution31.m b/membrane_solution31.m
new file mode 100644
index 0000000..1e1aca6
--- /dev/null
+++ b/membrane_solution31.m
@@ -0,0 +1,14 @@
+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
+    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)
+end 
\ No newline at end of file
diff --git a/script1.m b/script1.m
new file mode 100644
index 0000000..d675de7
--- /dev/null
+++ b/script1.m
@@ -0,0 +1,15 @@
+% 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
new file mode 100644
index 0000000..a1e15b6
--- /dev/null
+++ b/script2.asv
@@ -0,0 +1,15 @@
+% 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
new file mode 100644
index 0000000..fb841e8
--- /dev/null
+++ b/script2.m
@@ -0,0 +1,15 @@
+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