variation of parameters method