Finding the system output by convolution pretty new to this concept, so do bear with me.

A linear dynamic system is described by the following differential equation:

Transfer function H(s) is calculated to be  = 


I've already found the transfer function by assuming that the initial conditions are set to 0, getting the calculated stuff as shown above. I will then inverse Laplace it in order to find h(t), which turns out to be (exp^(-4t)t)

I will then use the formula, which is : 

This is what I got, and I was wondering if it is correct:

If so, why did the questions give me the initial conditions? What should I do with them? 
Also, how do I determine whether the system is stable? Unclear of this.

Thanks :)

As per EDIT:

I've converted them into differential equations where by: 

In this case, how do I find the TF? I'm supposed to divide by X(s) in both sides, but the equation on the right will have a problem won't it?

EDIT 2:
You can't compute the output by convolving the input with the system's impulse response if the initial conditions are non-zero. The reason is simple: there is no information about the initial conditions in the impulse response or, equivalently, in the system's transfer function.

What you can do is use the following properties of the unilateral Laplace transform:

$$\mathcal{L}\{y'(t)\}=sY(s)-y(0^+)\tag{1}$$

$$\mathcal{L}\{y''(t)\}=s^2Y(s)-sy(0^+)-y'(0^+)\tag{2}$$

If you use $(1)$ and $(2)$ when you transform the differential equation, you'll get an equation for the Laplace transform of $y(t)$ which takes into account the given initial conditions. The inverse Laplace transform of that expression gives you the desired system output for non-zero initial conditions.

As for stability, causal systems with a rational transfer function are stable if all poles (i.e., the zeros of the transfer function's denominator) are inside the left half-plane.