%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%% Display the Blasius flat plate velocity 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 velocity
plot(x_vals(:,2),eta_vals)

