A new, quasi-simultaneous method to calculate interacting boundary layers