• Remix
  • Share
  • New Entry

on 10 Oct 2022
  • 3
  • 25
  • 0
  • 0
  • 262
f=chebfun2(0);
rng('default')
for k=1:10
x0=2*rand-1;
y0=2*rand-1;
f=f+chebfun2(@(x,y)exp(-10*((x-x0).^2+(y-y0).^2)));
end
plot(f)
hold on
r=0.2.*roots(gradient(f));
colormap(summer)
plot3(r(:,1),r(:,2),f(r(:,1),r(:,2)),'k.','markersize',20)
% Some lighting almost always helps...
light
zlim([0 4])
axis off
view([-8 30]);
Remix Tree