• Remix
  • Share
  • New Entry

on 27 Oct 2021
  • 7
  • 24
  • 0
  • 0
  • 277
t=0:0.01:4;
k=(0:0.01:1)';
x=k.*t.*cos(2*pi*t);
y=k.*t.*sin(2*pi*t);
z=k.*flip(t);
figure(Color='k')
surf(z,y,x)
colormap autumn
shading interp
camlight
hold on
f=@(x,y,z)x^2+(y^2+z^2)/(.7-.1*x)^2-1;
fimplicit3(@(x,y,z)f((x+.5)/2,y/2,z/2),EdgeC='n',FaceC=[.96 .64 0.4]);
Warning: Function behaves unexpectedly on array inputs. To improve performance, properly vectorize your function to return an output with the same size and shape as the input arguments.
axis equal off
view(-90,90)
Remix Tree