1

cpl11001 · Dec 13, 2017 · c1713ed · c1713ed
1 parent 0ee309e
commit c1713ed
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 49 deletions.
diff --git a/SE_diff.asv b/SE_diff.asv
diff --git a/SE_diff.m b/SE_diff.m
@@ -8,7 +8,7 @@ function [pw_se,w] = SE_diff(T,P,n)
     m(2:n+1,2:n+1) = c; % insert matrix c into middle of m matrix
     % m = symmetric matrix 
 
-     h = 10/(n+1);
+     h = 10/(n+1); % h = delta_x = delta_y (micron)
 
     % Compute dw_dx
     dx = zeros(n+2,n+1);
@@ -38,7 +38,7 @@ function [pw_se,w] = SE_diff(T,P,n)
 
     pw_se = zeros(n+1,n+1);
 
-    for i = 1:16
+    for i = 1:(n+1)^2
         pw_se(i) = (1000*10^3*0.3*10^-3*h^2)/(2*(1-0.31^2))*((dw_dx(i))^4/4+(dw_dy(i))^4/4+0.5*(dw_dx(i))^2+0.5*(dw_dy(i))^2);
     end 
 
diff --git a/membrane_solution.m b/membrane_solution.m
@@ -9,7 +9,7 @@ function [w] = membrane_solution(T,P,n)% Set up initial matrix
 
     % 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
+    w = k\y; %solves for w in micron
 
     [x,y] = meshgrid(0:10/(n+1):10,0:10/(n+1):10);
     z = zeros(size(x));