function B=gaperture(rx,ry,x,y)
rx2=2*rx*rx;
ry2=2*ry*ry;
bx=exp(-(x.^2)/rx2);
by=exp(-(y.^2)/ry2);
B=by'*bx;