A numerical method for calculating unsteady two-dimensional laminar boundary layers