We prove that each (vector-valued) function in Sobolev space on a compact set K, which in the interior K0 of K satisfies a system of differential equations, can be approximated by solutions in a neighbourhood of K plus sums of potentials of measures supported on the boundary of K. We discuss the particular case where, for all compact sets K, one can dispense with potentials in such approximations