Update README.md

README.md

@@ -11,25 +11,32 @@ function [w] = membrane_solution3(T,P)
 % P = Pressure (MPa)
 % output:
 % w = vector of displacement of interior nodes
-    
-    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);
-   
-    
-    y = -(10/4)^2*(P/T)*ones(9,1);
-    w = k\y; 
-    % Solves for displacement (micrometers)
-    % Output w is a vector
-    % Solution represents a 2D data set w(x,y)
+
+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);
+
+
+y = -(10/4)^2*(P/T)*ones(9,1);
+w = k\y;
+
+% Solves for displacement (micrometers)
+% Solution represents a 2D data set w(x,y)
+
 
     [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)
-    title('Membrane Displacement')
-    zlabel('Displacement (micrometer)')
+    title('Membrane Displacement (3 x 3)')
+    zlabel('Z Position (micron)')
+    ylabel('Y Position (micron)')
+    xlabel('X Position (micron)')
     % Membrane displacement is shown on chart
+
+% Membrane displacement is shown on chart
+
+end
     
 end 
 ```