The new methods accurately integrate forced and damped oscillators. A family of analytical functions is introduced known as T-functions which are dependent on three parameters. The solution is expressed as a series of T-functions calculating their coefficients by means of recurrences which involve the perturbation function. In the T-functions series method the perturbation parameter is the factor in the local truncation error. Furthermore, this method is zero-stable and convergent. An application of this method is exposed to resolve a physic IVP, modeled by means of forced and damped oscillators. The good behavior and precision of the methods, is evidenced by contrasting the results with other-reputed algorithms implemented in MAPLE
