Solving numerically nonlinear systems of balance laws by multivariate sigmoidal functions approximation