Hi. I'm not very familiar with Green's functions, but let me suggest a direction I think you need to go:
Consider the linear differential equation (of one variable):
[tex]Ly=f[/tex]
where the L is the differential operator like in your case, it's [itex]L=\frac{d^2}{dx^2}+4[/itex]. Then we can show the solution can be written in terms of a Green's function as:
[tex]y(x)=\int G(x,u)f(u)du[/tex]
where:
[tex]G(x,u)=\sum_{n=1}^{\infty} \frac{\phi_n(x)\phi_n^{*}(u)}{\lambda_n}[/tex]
where [itex]\phi_n[/itex] is a orthonormal set of eigenfunctions for the operator L, that is, normalized solutions to the equation:
[tex]y''+4y=\lambda_n y[/tex]
subject to the given boundary conditions. So, just need to find those huh? Also, keep in mind the conjugate ([itex]\phi^{*}[/itex]) of a real-valued function is just the function.
See: "Mathematics of Classical and Quantum Physics" by F. Byron and R. Fuller. Whole chapter on Green's functions.