1

partg.m

@@ -1,5 +1,16 @@
-for P = [0.001:0.01:10]
-    [pw_se,w] = SE_diff(T,P,n);
-    [root,fx,ea,iter] = bisect_final_project(@(T) SE_diff(T,P,40),-5,10,.1);
-    disp(root)
-end
+
+P = linspace(0.001,0.01,10);
+
+for i = 1:length(P)
+    func = @(T) SE_diff(T,P(i),10);
+    [root(i),fx,ea,iter] = bisect_final_project(func,0.001,1,.1);
+end
+
+for i = 1:length(root)
+    w = membrane_solution(root(i),P(i),10);
+    w1(:,i) = w; % displays w values as column
+end
+
+w_max = max(w1);
+v_cub = interp1(w_max,P,'pchip');
+plot(w_max,P)