Clear[max, Pr, T, f, η, p];
max = 50;
Pr = 0.72;
pohl = NDSolve[{f'''[η] + 3 f[η] f''[η] -
     2 (f'[η])^2 + T[η] == 0,
   T''[η] + 2 Pr f[η] T'[η] == 0, f[0] == f'[0] == 0,
   f'[max] == 0, T[0] == 1, T[max] == 0}, {f, T}, {η, max}]

p4 = Plot[{Evaluate[f'[η] /. pohl]}, {η, 0, max},
  PlotRange -> All,
  PlotLabel ->
   Style[Framed["Hydrodynamic development is depicted on this plot"],
    10, Blue, Background -> Lighter[Yellow]], ImageSize -> Large,
  BaseStyle -> {FontWeight -> "Bold", FontSize -> 18},
  AxesLabel -> {"η", "f'[η]"}, PlotLegends -> "Expressions"]

Pr = 0.72;
pohl72 =
  NDSolve[{f'''[η] + 3 f[η] f''[η] - 2 (f'[η])^2 +
       T[η] == 0, T''[η] + 2 Pr f[η] T'[η] == 0,
    f[0] == f'[0] == 0, f'[max] == 0, T[0] == 1, T[max] == 0},
   {f, T}, {η, max}];

Pr = 0.6;
pohl = NDSolve[{f'''[η] + 3 f[η] f''[η] - 2 (f'[η])^2 + T[η] == 0,
   T''[η] + 2 Pr f[η] T'[η] == 0, f[0] == f'[0] == 0,
   f'[max] == 0, T[0] == 1, T[max] == 0},
  {f, T}, {η, max},
  Method -> {"Shooting",
    "StartingInitialConditions" ->
     Thread[{f[0], f'[0], f''[0], T[0], T'[0]} ==
       ({f[0], f'[0], f''[0], T[0], T'[0]} /. First@pohl72)]}]

Plot[{Evaluate[f'[η] /. pohl]}, {η, 0, max},
  PlotRange -> All,
  PlotLabel ->
   Style[Framed["Hydrodynamic development is depicted on this plot"],
    10, Blue, Background -> Lighter[Yellow]], ImageSize -> Large,
  BaseStyle -> {FontWeight -> "Bold", FontSize -> 18},
  AxesLabel -> {"η", "f'[η]"}, PlotLegends -> "Expressions"]

max = 50;
Pr = .72;
a = 0.7172594734816521`
b = -0.4344414944896132`
pohl = NDSolve[{f'''[[Eta]] + 3 f[[Eta]] f''[[Eta]] -
     2 (f'[[Eta]])^2 + T[[Eta]] == 0,
   T''[[Eta]] + 2 Pr f[[Eta]] T'[[Eta]] == 0, f[0] == f'[0] == 0,
   f''[0] == a, T[0] == 1, T'[0] == b}, {f, T}, {[Eta], max}]

b=.995b