function out = p(x,y)
      d1 = sqrt( (x+1).^2 + y.^2 );
      d2 = sqrt( (x-1).^2 + y.^2 );

      l= .5; k1 = 1; k2 = 2;

    out = k1*(d1-l).*(x+1)./d1 + k2*(d2-l).*(x-1)./d2  ;