scatter3(x1,y1,z1,'MarkerEdgeColor','b');
hold on;
for k=[10]
 F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
 zt = F([x1(k+1:k+10),y1(k+1:k+10)]);
 scatter3(x1(k+1:k+10),y1(k+1:k+10),zt, 'r','fill');
end
hold on;
for k=[60]
 F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
 zt = F([x1(k+1:k+10),y1(k+1:k+10)]);
 scatter3(x1(k+1:k+10),y1(k+1:k+10),zt, 'r','fill');
end
hold on;
for k=[120]
 F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
 zt = F([x1(k+1:k+10),y1(k+1:k+10)]);
 scatter3(x1(k+1:k+10),y1(k+1:k+10),zt, 'm','fill');
end
hold on;
% for k=[50]
%  F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
%  zt = F([x1(k+1:k+10),y1(k+1:k+10)]);
%  scatter3(x1(k+1:k+10),y1(k+1:k+10),zt, 'c','fill');
% end
% hold on;
% for k=[250]
%  F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
%  zt = F([x1(k+1:k+10),y1(k+1:k+10)]);
%  scatter3(x1(k+1:k+10),y1(k+1:k+10),zt, 'k','fill');
% end
% hold on;
% for k=[80]
%  F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
%  zt = F([x1(k+1:k+10),y1(k+1:k+10)]);
%  scatter3(x1(k+1:k+10),y1(k+1:k+10),zt, 'r','fill');
% end
% hold on;
% for k=[350]
%  F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
%  zt = F([x1(k+1:k+10),y1(k+1:k+10)]);
%  scatter3(x1(k+1:k+10),y1(k+1:k+10),zt, 'r','fill');
% end
% hold on;
% for k=[360]
%  F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
%  zt = F([x1(k+1:k+20),y1(k+1:k+20)]);
%  scatter3(x1(k+1:k+20),y1(k+1:k+20),zt, 'g','fill');
% end
% hold on;
% 
% for k=[380]
%  F = scatteredInterpolant(x1(1:k),y1(1:k),z1(1:k), 'linear');
%  zt = F([x1(k+1:k+10),y1(k+1:k+10)]);
%  scatter3(x1(k+1:k+10),y1(k+1:k+10),zt, 'g','fill');
% end
% hold on;