Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
Part G
  • Loading branch information
nin13001 committed Dec 16, 2017
1 parent 7785078 commit ca8a423
Show file tree
Hide file tree
Showing 4 changed files with 49 additions and 0 deletions.
22 changes: 22 additions & 0 deletions PartG.m
@@ -0,0 +1,22 @@
No = linspace(.001,.01,10);
n = 20;

T = zeros(1,length(No));
wmax = zeros(1,length(No));

for i = 1:length(No)
T(i) = tension_sol(No(i),n);
w = membrane_solution(T(i),No(i),n);
wmax(i) = max(w);
end

clf
x = wmax';
y = No';
Z=x.^3;
a=Z\y;
x_fcn=linspace(min(x),max(x));
plot(x,y,'o',x_fcn,a*x_fcn.^3)
title('Pressure vs Maximum Deflection')
xlabel('Maximum Deflection (\muM)')
ylabel('Pressure (MPa)')
27 changes: 27 additions & 0 deletions README.md
Expand Up @@ -154,3 +154,30 @@ ea = zeros(1,length(n));
|30|0.0602|0.16%|
|35|0.0602|0.09%|
|40|0.0603|0.05%|

## Part G
```
No = linspace(.001,.01,10);
n = 20;
T = zeros(1,length(No));
wmax = zeros(1,length(No));
for i = 1:length(No)
T(i) = tension_sol(No(i),n);
w = membrane_solution(T(i),No(i),n);
wmax(i) = max(w);
end
clf
x = wmax';
y = No';
Z=x.^3;
a=Z\y;
x_fcn=linspace(min(x),max(x));
plot(x,y,'o',x_fcn,a*x_fcn.^3)
title('Pressure vs Maximum Deflection')
xlabel('Maximum Deflection (\muM)')
ylabel('Pressure (MPa)')
```
![](assets/README-bbfc19b7.png)
Binary file added assets/PartG.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file added assets/README-bbfc19b7.png
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.

0 comments on commit ca8a423

Please sign in to comment.