- /
-
Another Torpedo
on 21 Oct 2021
- 3
- 19
- 0
- 0
- 280
f=@(x,y,z)x^2+99*(y^2+z^2)/(x-1)^2-1;
figure(Color='#27b');
p=fimplicit(@(x,y)f(x,y,0));
X=p.XData;
Y=p.YData;
i=[220 402];
X(i)=1;
Y(i)=[.1 -.1];
a={Y,0,'EdgeC','n','FaceC','#bde'};
fill(X,a{:});
hold;
fill3(X,X-X,a{1:4});
fimplicit3(f,a{3:6});
axis equal off;
view(i-255)
camlight;
material metal
xlim([-1,.6]);