function body(R)
% body(R) : contour dividing streamline for plate boundary layer
%           for a given Reynolds number R
%           edit m,n,h to change resolution/contour values
m=201;
n=101;
R=500;
a=i*1.2165*2^0.5/R^0.5;
dx=16./200;
dy=0.5/100;
for v1=1:m,
	for v2=1:n,
		x=-1.+v1*dx;
		y=v2*dy;
		z=x+i*y;
		x1(v2,v1)=x;
		y1(v2,v1)=y;
		w(v2,v1)=z-a*z^0.5;
	end
end
psi=imag(w);
phi=real(w);
h=[0.0 0.0];
contour(x1,y1,psi,h);
set(gca,'FontName', 'Times','FontSize',12);
xlabel('x','FontSize',12);
ylabel('y','FontSize',12);