Blasius Boundary Layer
Note that for my third boundary condition, I used 10,000,000 as a placeholder for infinity...

ClearAll["Global`*"]
eqs={ 2f'''[η]+f[η] f''[η] ==0, f[0]==0,f'[0]==0,f'[10000000.0]==1};
f[η_]=f[η]/.NDSolve[eqs,f[η],{η,0,20}][[1]];
g0[η_]=D[f[η],η]; c=η/.FindRoot[g0[η]-0.99,{η,5}];
η[g_]=InverseFunction[g0][g];
Plot[η[g],{g,0,1},PlotRange->{{0,1.2},{0,8}},AspectRatio-> 1,PlotStyle->Directive[RGBColor[1,0,0],AbsoluteThickness[3]],ImageSize->180,Frame->True,FrameLabel->{"(SuperStar[u])","η"}]
Print["The coefficient for δ at 0.99 \!\(\*SubscriptBox[\(U\), \(\(0\)\(\\ \)\)]\) is C = ",c, " so (δ/x)=",c," ((Subscript[Re, x])^(-0.5))"]