method='weightedbergman';
alpha=2;
xm=-4;
xM=4;
ym=-4;
yM=4;
npts=100;
deg=40;
x=linspace(xm,xM,npts);
y=linspace(ym,yM,npts);
[X,Y]=meshgrid(x,y);
pts=[X(:),Y(:)];
% true solution
Vtrue1=SiciakRealDisk(X,Y);
surf(X,Y,Vtrue1)