Solve the following system:

X' = ([[1,-1],[1,-1]])X + ([1/t],[1/t]), X(1) = ([2],[-1])

using X = phi(t)*phi^(-1)(t_0)X_0 + phi(t)*integral[from t_0 to t](phi^(-1)(s)F(s)ds)