function B=caperture(x,y,r)
r2=2*r*r;
z=exp(-0.5);
bx=exp(-(x.^2)/r2);
by=exp(-(y.^2)/r2);
A=by'*bx;
C=A-z;
B=ceil(C)+floor(C);