       function f=f(x,y)
            pnt=[x,y];
            fxy=logan_shepp(pnt,0,5);
       f=fxy;