diff --git a/README.md b/README.md index 9fbf9b9..2d6f30a 100644 --- a/README.md +++ b/README.md @@ -1 +1,30 @@ -# ME3255_FInalProject \ No newline at end of file +# 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 +``` diff --git a/membrane_solution3.m b/membrane_solution3.m new file mode 100644 index 0000000..ca1d139 --- /dev/null +++ b/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 \ No newline at end of file