不要問怎麼算，
我用速算法，
正規算法請參考原文書的二階常係數 O.D.E. 之
undetermined coeffient 或是 variation of parameter，
當然你要用 Laplace transform 或是矩陣也可以解。