我们要解的方程是f'' + (1/2)*f*f, f(0) = 0, f(0) = 0, f'(inf) = 1。这个方程是在流体边界层理论中产生的，必须用数值方法求解。我们将这个问题重新定义为一个一阶ode系统:y = [f;f ';f] = [y (1);y (2);y(3) = y' = [f';f”;f”]= [y (2);y (3); -(1/2)*y(1)*y(3)] with y(1)(0) = 0, y(2)(0) = 0, y(2)(inf) = 1. This new system of equations may then be solved numerically using the shooting method. For a description of this numerical method, please follow the given link:https://en.wikipedia.org/wiki/Shooting_method．