torpedo

Remix
Share
New Entry

Jan Studnicka
/
torpedo

on 18 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,'FaceC','#bde','EdgeC','n'};

fill(X,a{:});

hold;

Current plot held

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

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(150,30);

light(Po=[-2 1 1]);

xlim([-1,.6]);