Text this: Numerical integration by the matrix method of boundary value problems for linear inhomogeneous ordinary differential equations of the third order with variable coefficients