Given y'=[matrix]y + [matrix of constants]. y(0)=[initial conditions] Would I treat the matrix with variation of parameters? Or should I ignore the constants on right hand side and find the general solution. Then when it comes to solving the initial conditions, then do the calculations with it? Thanks in advance.