- /
-
Heidi in the Sun
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]);