hold on
N = 1000;
xv = linspace(min(x), max(x), N);
yv = linspace(min(y), max(y), N);
[Xm,Ym] = ndgrid(xv, yv);
Zm = griddata(x, y, F1, Xm, Ym);
figure
surf(Xm, Ym, Zm)
shading interp
grid on
N = 7.5;
xv = linspace(min(xr), max(xr), N);
yv = linspace(min(yr), max(yr), N);
[Xr,Yr] = ndgrid(xv, yv);
Zr = griddata(xr, yr, F1r, Xr, Yr);
figure
plot3(Xr, Yr, Zr,'o','Color','g')
grid on
hold off