The WKB analysis for higher order differential equations had not been developed until recent progress of the so-called exact WKB method, so it was difficult to find WKB solutions for multistate nonadiabatic models in general, whereas enormous efforts were made especially for the study of generalized Landau-Zener models. Under such circumstances Aoki-Kawai-Takei provided a closed expression for the $S$ matrix of a time-dependent multistate nonadiabatic system, in which diabatic potential curves are given as polynomial functions of time and they all intersect with each other. Their result is completely explicit, and the model covers much wider classes of systems than multistate nonadiabatic models studied so far. Here we numerically show how their exact WKB formula works as long as crossing points of diabatic potential curves are well separated. The model studied includes the cases where diabatic potential functions are nonlinear, not like generalized Landau-Zener models. We also discuss the situations where virtual turning points and new Stokes curves, new components in the Stokes geometry appearing only for higher order differential equations, come into play. In particular, we examine how the presence of new Stokes curves affects final transition probabilities using a couple of concrete examples.