vpasolve() uses the MuPAD symbolic Engine's numeric:;solve routine. Depending on the form of the nonlinearity, numeric::polysysroots might be invoked, but more likely numeric::fsolve will be invoked. For that,
"numeric::fsolve implements a purely numerical Newton type root search with a working precision set by the environment variable DIGITS"