%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%% Display the Blasius flat plate vorticity profile %%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% scaled wall shear
tau_wall=0.33206;

% integrate the ODE
[eta_vals, x_vals] = ode45(@blaseq, [0 10], [0; 0; tau_wall]);

% plot the vorticity
plot(x_vals(:,3)+0.000001,eta_vals)

