Back to Gallery
  • Remix
  • Share
  • New Entry
on 6 Oct 2022
  • 1
  • 12
  • 0
  • 0
  • 277
k=@(x)kron(x,[1;1;1]);
e=@(x)((sum(x.^2,2)).^.5-1).^.5+x(:,3).^4-1;
E=2e3;
x=[1,0,1];
V=[1,1,1];
s=E;
for j = 0:E
p=x(end,:);
d=e(p);
if d>-1/E
n=(e(k(p)+eye(3)/E).'-d)*E;
n=n/norm(n);
V=V+(-2*dot(n,V)*n);
end
x=[x;p+V/10];
s=[s;randi(E)];
end
scatter(x(:,2),x(:,1),s/1e2,s,'|');
colormap(turbo);
axis equal off;
Remix Tree