We prove that if a vector-function f belongs to the Morrey space L 1,λ (Ω, R N ), with Ω ⊂ R n , n ≥ 3, N ≥ 2, λ ∈ [0, n − 2], then there exists a very weak solution u of the systemsuch that Du belongs to the space L q,n−q(n−λ−1) loc (Ω, R nN ) for any q ∈ [1, n n−1 [, provided the matrix of coefficients (A ij ) has L ∞ ∩ VMO entries.