%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%% Display the Blasius flat plate wall shear %%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% flow velocity
U=1;

% length of the plate
l=1;

% density
rho=1;

% scaled wall shear
fpp0=0.33206;

% Renolds number based on the length
Re_l=100;

% viscosity
nu=U*l/Re_l;

% x-limits
x_lo=U^3*nu*(fpp0*rho)^2+0.00001;
x_hi=0.99999;

% x-values
x_vals=x_lo:(x_hi-x_lo)/100:1;

% plot
plot(x_vals,fpp0*rho*sqrt(U^3*nu./x_vals),'r')
hold on

% Renolds number based on the length
Re_l=1000;

% viscosity
nu=U*l/Re_l;

% x-limits
x_lo=U^3*nu*(fpp0*rho)^2+0.00001;
x_hi=0.99999;

% x-values
x_vals=x_lo:(x_hi-x_lo)/100:1;

% plot
plot(x_vals,fpp0*rho*sqrt(U^3*nu./x_vals),'g')

% Renolds number based on the length
Re_l=10000;

% viscosity
nu=U*l/Re_l;

% x-limits
x_lo=U^3*nu*(fpp0*rho)^2+0.00001;
x_hi=0.99999;

% x-values
x_vals=x_lo:(x_hi-x_lo)/100:1;

% plot
plot(x_vals,fpp0*rho*sqrt(U^3*nu./x_vals),'b')
hold off
