Start from the boundary layer equations for laminar, steady, two-dimensional and incompressible flow over a flat plate to get the Blasius solution:

1. apply the similarity variables and derive the boundary layer transformed ordinary differential equation (ODE)
2. reduce the ODE into system of first order ODEs
3. solve the above system using RK-4th or Adams-Bashforth 4th order methods together with a shooting technique with the help of methods such as Newton-Raphson or bisection methods
4. write a Matlab code to implement the methods in part 3