Update README.md

ltd13002 · Dec 14, 2017 · 310a840 · 310a840
1 parent 03a545d
commit 310a840
Showing 1 changed file with 19 additions and 1 deletion.
diff --git a/README.md b/README.md
@@ -176,4 +176,22 @@ end
 
 
 # Part G
-
+```matlab
+P = linspace(.001,.01,10);
+n = 20;
+T = zeros(1,length(P));
+wmax = zeros(1,length(P));
+for i = 1:length(P)
+    T(i) = tension_sol(P(i),n);
+    w = membrane_solution(T(i),P(i),n);
+    wmax(i) = max(w);
+end
+clf
+setDefaults
+x = wmax';
+y = P';
+Z=x.^3;
+a=Z\y;
+x_fcn=linspace(min(x),max(x));
+plot(x,y,'o',x_fcn,a*x_fcn.^3)
+```