In this research we introduce a new variation of parameters for systems of linear and nonlinear ordinary differential equations. We use known mathematical methods and techniques including Gronwall’s inequality and fixed point theory to obtain boundedness on all solutions and stability results on the zero solution.