I have been doing a simulation of Blasius equation: F'''+FF''/2=0 with F(eta) where eta is a similarity variable eta=(y-1)/(x^(1/2)) F'(0)=1 u(y=1,x<<1)=1 F(0)=0 v(y=1,x<<1)=0 F'(infinite)=0 u(y=infinite, x<<1)=0 You can observe that the BC's are different of the flow over flat plate. This is obtained for a mixing thickness in the near field of a two-dimensional jet, near the orifice of exhaust (are you agree?). in this page: http://www.rit.edu/~pnveme/Matlab_Course/Matlab_App_ODE.html where shooting method is employed in Matlab, it is said numerical methods (Runge Kutta, or the internal Matlab function ODE45) have normalized value F'(infinity)=1. I have programmed it in Matlab but now I don't know how consider the value 0 at infinity instead of 1. Could you help me?. Thanks.