%Yeonjin Jung
%ME3165
%02/13/2018

clear;
clc;
b=.5;
n=1000;
dy=b/n;
L=.1;
rho=1;
u=zeros(n+1,5);
u(1,:)=0;
u(n+1,:)=0.03;
dpdx=[-1e-3 -.5e-3 0 .5e-3 1e-3];
vis=8.9e-4;

a=ones(n,5)*1;
b=ones(n,5)*(-2);
c=ones(n,5);
d=ones(n,length(dpdx))*dy^2*(1/vis).*dpdx;


for z=1:5
    d(n,z)=d(n,z)-u(n+1,z);

for k=2:n
    m=a(k,z)/(b(k-1,z));
    b(k,z)=b(k,z)-m*c(k-1,z);
    d(k,z)=d(k,z)-m*(d(k-1,z));

end
u(n,z)=d(n,z)/b(n,z);
for k=n-1:-1:5
    u(k,z)=(d(k,z)-c(k,z)*u(k+1,z))/b(k,z);

end
end
figure
plot(u,0:n)
title('Flow between Parallel pipes')
xlabel('v in m/s')
ylabel('nodes')
ren=rho*L/vis.*u;