solve x' = -x + 4y + e^t, y' = x-y+1