Near Field Power Density of Parabolic Reflector

Dear All,

Do you have any information or ideas about how can I plot in MATLAB the power density in the near field of parabolic reflector?

In order to get this plot:

Thanks in advance!