A new method for computing a p-solution to parametric interval linear systems with affine-linear and nonlinear dependencies