From 08aa3968710f97075c5b7b3df867a0d97e563b49 Mon Sep 17 00:00:00 2001
From: mattmaliniak <31718098+mattmaliniak@users.noreply.github.com>
Date: Thu, 14 Dec 2017 15:28:51 -0500
Subject: [PATCH] Part F

---
 Part F/F_sol.m       | 2 +-
 Part F/SE_diff.m     | 6 +++---
 Part F/tension_sol.m | 2 +-
 3 files changed, 5 insertions(+), 5 deletions(-)

diff --git a/Part F/F_sol.m b/Part F/F_sol.m
index e7e8a13..1a749ca 100644
--- a/Part F/F_sol.m	
+++ b/Part F/F_sol.m	
@@ -1,4 +1,4 @@
- n=20:5:40;
+ n=[3,20:5:40];
  P=0.001; %MPa
  T = zeros(1,length(n));
  ea = zeros(1,length(n));
diff --git a/Part F/SE_diff.m b/Part F/SE_diff.m
index 9338a84..bf9e129 100644
--- a/Part F/SE_diff.m	
+++ b/Part F/SE_diff.m	
@@ -18,9 +18,9 @@ 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)]);
+        dwdx(i,j) = mean([z(i+1,j)-z(i,j),z(i+1,j+1)-z(i,j+1)])./h;
+        dwdy(i,j) = mean([z(i,j+1)-z(i,j),z(i+1,j+1)-z(i+1,j)])./h;
     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 = abs(pw-se);
\ No newline at end of file
+pw_se = pw-se;
\ No newline at end of file
diff --git a/Part F/tension_sol.m b/Part F/tension_sol.m
index dd534c9..bc38e29 100644
--- a/Part F/tension_sol.m	
+++ b/Part F/tension_sol.m	
@@ -1,3 +1,3 @@
 function [T,ea] = tension_sol(P,n)
 y =@(T) SE_diff(T,P,n);
-[T,fx,ea,iter]=bisect(y,.01,1);
+[T,fx,ea,iter]=bisect(y,.01,1,.1);