# load simple 3um fiber mode profile
loaddata("bragg_fiber_simple_3um");

# overlap offset vectors
npts=51;
x=linspace(-3e-6,3e-6,npts);
y=x;
overlap_XY=matrix(npts,npts);

for (i=1:npts) {
  for (j=1:npts) {
    temp_val = overlap("mode3","global_mode1",x(i),y(j),0);
    overlap_XY(i,j) = temp_val(1);
  }
  ?i;  
}
image(x*1e6,y*1e6,overlap_XY,"x (um)","y (um)","Overlap");