Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
Part A
  • Loading branch information
nin13001 committed Dec 15, 2017
1 parent d38552a commit 7535481
Show file tree
Hide file tree
Showing 2 changed files with 54 additions and 1 deletion.
31 changes: 30 additions & 1 deletion README.md
@@ -1 +1,30 @@
# ME3255_FInalProject
# ME3255_FInalProject
I noticed that I was never assigned to a group so I did this project solo
## Part A

```
function [w] = membrane_solution3(T,P)
od = ones(8,1);
od(3:3:end) = 0;
k = -4*diag(ones(9,1))+diag(ones(9-3,1),3)+diag(ones(9-3,1),-3)+diag(od,1)+diag(od,-1);
y = -(10/4)^2*(P/T)*ones(9,1);
w = k\y;
grid on
[x,y] = meshgrid(0:10/4:10,0:10/4:10);
z = zeros(size(x));
z(2:end-1,2:end-1) = reshape(w,[3 3]);
surf(x,y,z)
title('Membrane Displacement')
xlabel('X (\muM)')
ylabel('Y (\muM)')
zlabel('Displacement (\muM)')
disp(w)
end
```
24 changes: 24 additions & 0 deletions membrane_solution3.m
@@ -0,0 +1,24 @@
function [w] = membrane_solution3(T,P)


od = ones(8,1);
od(3:3:end) = 0;
k = -4*diag(ones(9,1))+diag(ones(9-3,1),3)+diag(ones(9-3,1),-3)+diag(od,1)+diag(od,-1);


y = -(10/4)^2*(P/T)*ones(9,1);
w = k\y;

grid on
[x,y] = meshgrid(0:10/4:10,0:10/4:10);
z = zeros(size(x));
z(2:end-1,2:end-1) = reshape(w,[3 3]);
surf(x,y,z)
title('Membrane Displacement')
xlabel('X (\muM)')
ylabel('Y (\muM)')
zlabel('Displacement (\muM)')

disp(w)

end

0 comments on commit 7535481

Please sign in to comment.