function z=logan(x,y)

pts=[x y];

p1=0;
q1=0;

fxy1=logan_shepp(pts,p1,q1);

z=fxy1;