2D Near to far field transformation in Matlab

Hi, I'm having some problems in implementing a 2D near to far field transformation with Matlab. I'm following Taflove's book, but I can't get it work!
I have few questions:

- Do I need both E and H tangential fields, or I can assume a PEC or PMC inside my surface and then use just one of them?

- I'm using the fft matlab function to get the fields in phasor form, is this correct?

- About the exponential term to correct the phase inside the integral, I'm using this on the first side (the left side, for the other three sides I'm using similar terms):
X0*cos(phi) + DY.*sin(phi)
where X0 is the x coordinate and DY goes from Y0 to Y1; is it correct?

- How can I implement the integration? trapz function is ok?

Thanks a lot to everybody