Green's function implementation

Hello,

How can I implement the Green's function to find the internal fields of an infinite slab

coditions:

1- one dimentional issue

2- inciednt wave comes from one side (thus, just 2 kernels (G11,G21))

I need to know also information about the intial, jump and boundary conditions and how it is implemented

providing code will be great

Thank you