Numerical solution of linear differential equations with discontinuous coefficients and Henstock integral

In this article we consider the problem of approximative solution of linear differential equations $y^{\mathrm{\prime }}+p(x)y=q(x)$ with discontinuous coefficients $p$ and $q$. We assume that coefficients of such equation are Henstock integrable functions. To find the approximative solution we change the original Cauchy problem to another problem with piecewise-constant coefficients. The sharp solution of this new problems is the approximative solution of the original Cauchy problem. We find the degree approximation in terms of modulus of continuity ${\omega }_{\delta }(P),{\omega }_{\delta }(Q)$, where $P$ and $Q$ are $f$-primitive for coefficients $p$ and $q$.