A numerical procedure (integral method) for calculating the interacting turbulent boundary layer is set up. With this method, some free interactions with various upstream conditions are simulated in order to investigate the influence of upstream state on the interacting turbulent boundary layer. The results obtained by this numerical simulation can be summarized as follows; Free interaction of upstream unstabilized (or separated) turbulent boundary layer is subcritical regardless of its external Mach number, while free interaction of upstream stabilized turbulent boundary layer has two different characteristics (subcritical, supercritical) according to the external Mach number.