@bbgodfrey has explained very well his approach in the answer to this question. But, In my answer, I am trying to solve both the equations simultaneously.

Eqn1 = f'''[x] + f[x] f''[x] + 4 - (f'[x])^2 == 0

Eqn2 = T''[x] + Pr f[x] T'[x] == 0

BC1 = f[0] == 0;

BC2 =  f'[0] == 0;

BC3 = f'[inf1] == 2;

BC4 = T'[0] == -1;

param1 = {Pr -> 3.97};

inf1 = 5;

Sol1 = NDSolve[{Eqn1, Eqn2, BC1, BC2, BC3, BC4, BC5} /. param1, {f, 
   T}, {x, 0, inf1}, Method -> {"Shooting", 
"StartingInitialConditions" -> {f[0] == 0, f'[0] == 0, f''[0] == 3.48,
 T[0] == 0, T'[0] == -10}}]

Plot[{f'[x] /. Sol1, T[x] /. Sol1}, {x, 0, inf1}, PlotRange -> All, 
 PlotStyle -> {Black, Red}, Frame -> True, 
 FrameStyle -> Directive[Black, Bold, 12], PlotRange -> All, 
 Axes -> False]

@bbgodfrey has explained very well his approach in the answer to this question. But, In my answer, I am trying to solve both the equations simultaneously.

Eqn1 = f'''[x] + f[x] f''[x] + 4 - (f'[x])^2 == 0

Eqn2 = T''[x] + Pr f[x] T'[x] == 0

BC1 = f[0] == 0;

BC2 =  f'[0] == 0;

BC3 = f'[inf1] == 2;

BC4 = T'[0] == -1;

BC5 = T[inf1] == 0;

param1 = {Pr -> 3.97};

inf1 = 5;

Sol1 = NDSolve[{Eqn1, Eqn2, BC1, BC2, BC3, BC4, BC5} /. param1, {f, 
   T}, {x, 0, inf1}, Method -> {"Shooting", 
"StartingInitialConditions" -> {f[0] == 0, f'[0] == 0, f''[0] == 3.48,
 T[0] == 0, T'[0] == -10}}]

Plot[{f'[x] /. Sol1, T[x] /. Sol1}, {x, 0, inf1}, PlotRange -> All, 
 PlotStyle -> {Black, Red}, Frame -> True, 
 FrameStyle -> Directive[Black, Bold, 12], PlotRange -> All, 
 Axes -> False]

I hope this helps.