This paper is concerned with weak solutions of elliptic equations and of systems of the form \(\begin{aligned} {\left\{ \begin{array}{ll} \displaystyle -\sum _{i=1}^{n}D_i \left( A^\alpha _i\left( x,Du(x)\right) \right) = 0, \ \ & \textrm{in }\ \Omega , \ \ \forall \ \alpha \in \left\{ 1, \cdots , m\right\} ,\\ u(x)=u_*(x), & \textrm{on }\ \partial \Omega . \end{array}\right. } \end{aligned}\) We show that, assuming some high degree of integrability of \(Du_*\) the gradient of the boundary datum, we can obtain bounds of the difference \( |u^\gamma -u^\gamma _* |_\infty \) , \(\gamma \in \left\{ 1, \cdots , m\right\} \) . In particular, when \(m=1\) , things are easier and we can get a better result.