diff --git a/README.md b/README.md index d8cc217..25ffb94 100644 --- 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) +```