Abstract Since obtaining an analytic solution to some mathematical and physical problems is often very difficult, academics in recent years have focused their efforts on treating these problems using numerical methods. In science and engineering, systems of integral differential equations and their solutions are extremely important. The Taylor collocation method is described as a matrix approach for solving numerically Linear Differential Equations (LDE) by using truncated Taylor series. Integral equations are used to solve problems such as radiative transmission and the oscillation of a string, membrane, or axle. Differential equations can be used to tackle oscillating difficulties. To discover approximate solutions for linear systems of integral differential equations with variable coefficients in terms of Taylor polynomials, the collocation approach, which is offered for differential and integral equation solutions, will be developed. A system of LDE will be translated into matrix equations, and a new matrix equation will be generated in terms of the Taylor coefficients matrix by employing Taylor collocation points. The needed system will be converted to a linear algebraic equation system. Finding the Taylor coefficients will lead to the Taylor series technique.