diff --git a/README.md b/README.md
index 991d45b..bf372d9 100644
--- a/README.md
+++ b/README.md
@@ -71,3 +71,34 @@ end
 [w] = membrane_solution(0.006,0.001,10)
 ```
 ![](https://github.uconn.edu/ltd13002/ME3255_Final_Project/blob/master/Part%20D/PartDFigure.png)
+
+
+# Part E
+```matlab
+function [pw_se,w]=SE_diff(T,P,n)
+E = 1; %TPa Units may need to be changed
+v = .31; %Poissons ratio
+t = .3; %nm
+h = 10/(n+1); %nm
+w = membrane_solution(T,P,n);
+z = zeros(n+2);
+z(2:end-1,2:end-1) = reshape(w,[n n]);
+num = n + 1;
+wbar = zeros(num);
+for i = 1:num
+    for j = 1:num
+        wbar(i,j) = mean([z(i,j),z(i+1,j),z(i,j+1),z(i+1,j+1)]);
+    end
+end
+pw = sum(sum(wbar.*h^2.*P));
+dwdx = zeros(num);
+dwdy = zeros(num);
+for i = 1:num
+    for j = 1:num
+        dwdx(i,j) = mean([z(i+1,j)-z(i,j),z(i+1,j+1)-z(i,j+1)]);
+        dwdy(i,j) = mean([z(i,j+1)-z(i,j),z(i+1,j+1)-z(i+1,j)]);
+    end
+end
+se = E*t*h^2/(2*(1-v^2))*sum(sum(0.25.*dwdx.^4+.25.*dwdy.^4+0.5.*(dwdx.*dwdy).^2));
+pw_se = pw-se;
+```