For a zero-pressure gradient boundary layer simulations, compressible similarity solution to the boundary layer equations can be obtained using Dorodnitsyn transofrmation which consists of a system of ODEs. In this function, we solve this system using MATLAB ode solver and the bisection method. The user needs to specify the Mack number of the flow and the Wall temperature BC (or adiabatic if needed).

Due to small temperature variation, a linear dependence of viscosity on temperature is assumed for simplicity. For more info and help email me maslani@fsu.edu