diff --git a/Part G/Part g.png b/Part G/Part g.png new file mode 100644 index 0000000..7f2d72c Binary files /dev/null and b/Part G/Part g.png differ diff --git a/README.md b/README.md index 76b1062..25ffb94 100644 --- a/README.md +++ b/README.md @@ -110,7 +110,8 @@ pw_se = pw-se; T = zeros(1,length(n)); ea = zeros(1,length(n)); for i = 1:length(n) - [T(i), ea(i)] = tension_sol(P,n(i)); + +[T(i), ea(i)] = tension_sol(P,n(i)); end ``` ```matlab @@ -172,3 +173,25 @@ end |30|0.0602|0.15%| |35|0.0602|0.09%| |40|0.0603|0.06%| + + +# 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) +```