Another Torpedo

Remix
Share
New Entry

Jan Studnicka
/
Another Torpedo

on 21 Oct 2021

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));

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.

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;

Current plot held

fill3(X,X-X,a{1:4});

fimplicit3(f,a{3:6});

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(i-255)

camlight;

material metal

xlim([-1,.6]);