close('all')
 clear

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Load data
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 p =load('blasius_solutions.dat');
     x = p(:,1);
     y = p(:,2);
   eta = p(:,3);
     u = p(:,4);
     v = p(:,5);
    f  = p(:,6);
    f1 = p(:,7);
    f2 = p(:,8);
    nu = p(:,9);
 U_inf = p(:,10);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 1. u(eta) plot
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  figure(1)
  hold on
  plot(u,eta, 'k-')
  xlabel('u/Uinf')
  ylabel('\eta')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2. v(eta) plot
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  figure(2)
  hold on
  plot(v,eta, 'k-')
  xlabel('v/Uinf')
  ylabel('\eta')