function z=logan_shepp_product(x,y)

pts=[x y];

p1=0;
q1=0;
p2=0;
q2=0;

fxy1=logan_shepp(pts,p1,q1);
fxy2=logan_shepp(pts,p2,q2);

z=fxy1.*fxy2;