  The scaling of the exact solution of a hyperbolic balance law generates a
family of scaled problems in which the source term does not depend on the
current solution. These problems are used to construct a sequence of solutions
whose limiting function solves the original hyperbolic problem. Thus this gives
rise to an iterative procedure. Its convergence is demonstrated both
theoretically and analytically. The analytical demonstration is in terms of a
local in time convergence and existence theorem in the $L^2$ framework for the
class of problems in which the source term $s(q)$ is bounded, with $s(0) = 0$,
is locally Lipschitz and belongs to $C^2(\mathbb{R}) \cap H^1 (\mathbb{R}) $. A
convex flux function, which is usual for existence and uniqueness for
conservation laws, is also needed. For the numerical demonstration, a set of
model equations is solved, where a conservative finite volume method using a
low-dissipation flux is implemented in the iteration stages. The error against
reference solutions is computed and compared with the accuracy of a
conventional first order approach in order to assess the gaining in accuracy of
the present procedure. Regarding the accuracy only a first order scheme is
explored because the development of a useful procedure is of interest in this
work, high-order accurate methods should increase the computational cost of the
global procedure. Numerical tests show that the present approach is a feasible
method of solution.
