Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
  • Loading branch information
ltd13002 committed Dec 14, 2017
2 parents 9675c3d + 9f45a13 commit eb672e6
Show file tree
Hide file tree
Showing 2 changed files with 24 additions and 1 deletion.
Binary file added Part G/Part g.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
25 changes: 24 additions & 1 deletion README.md
Expand Up @@ -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
Expand Down Expand Up @@ -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)
```

0 comments on commit eb672e6

Please sign in to comment.