Back to Gallery
  • Remix
  • Share
  • New Entry
on 6 Nov 2023
  • 9
  • 17
  • 1
  • 0
  • 478
drawframe(1);
Write your drawframe function below
function drawframe(f)
set(gcf,'color','k')
ax=newplot;
set(ax,'position',[0 0 1 1]);
axis(ax,'off','tight');
xlim(ax,[0 30])
daspect(ax,[1 1 1])
qs=linspace(-8,9,48);
q=qs(f);
if q<=8
c=[ q 0 0
-q 1 0];
else
d=9-q;
c=[ 8 1-d 0
-8 d 0];
end
ct=[15 10 10]+c;
nx=30;
ny=20;
nz=20;
a=.05;
b=1;
[y,x,z]=ndgrid(1:ny,1:nx,1:nz);
B=zeros(ny,nx,nz);
for i=1:size(ct,1)
B=B+b*exp(-a*((ct(i,1)-x).^2 + (ct(i,2)-y).^2 + (ct(i,3)-z).^2));
end
S=isosurface(B,.8);
p=patch(ax,S,'FaceC', '#f71', 'EdgeC', 'n','FaceL', 'g');
material(ax,[.6 .9 .3 2 .5])
camlight(ax)
end
Animation
Remix Tree