Fixing bugs in the boundary layer correction and RK2.

Fixing a bug in the initial value for the momentum thickness.
Replacing the np.trapz integration by an RK2 method.