%	Developing Couette flow, 0<y<L
clear all
clf reset
y=[0: 0.01: 1];
t=[.01 .1 1];
%	This is kinematic viscosity times time
v=y;
nsum=5;
for j=1:length(t)
   v=erfc((1-y)./sqrt(4*t(j)))-erfc((1+y)./sqrt(4*t(j)));
   for n=1:nsum
   	v=v+erfc(((2*n+1)-y)./sqrt(4*t(j)))-erfc(((2*n+1)+y)./sqrt(4*t(j)));
   	plot(y,v,'k')
   	xlabel('y')
   	ylabel('Velocity')
   	title(['Developing Couette Flow; t= ',num2str(t)])
      hold on
   end
end
load trans.dat
JMAX=21;
del=1/(JMAX-1);
x=linspace(0,1,21);
xp=2.*(x-0.5);
time=[.0025 .025 .25];
for j=1:length(time)
   plot(xp,trans(j,:),'r+')
   axis([0 1 0 1])
   hold on
end

