@Obliv, actually I have done problems using variation of parameters for DE with constant coefficient, from
here , and they mentioned that one of the disadvantages of variation of parameter method is that complementary solution is must to use this method.