function ii=plsph_public(rr,sinu,xx,yy)

% 2.997 Optical Engineering / George Barbastathis
% Massachusetts Institute of Technology
% 10/12/1999

sph=0; % replace with expression for spherical wave
pl=0;  % replace with expression for plane wave
ii=abs(sph+pl).^2;
imagesc(ii)