submarine

Remix
Share
New Entry

Jan Studnicka
/
submarine

on 15 Oct 2021

f=@(x,y,z)x^2+1e2*(y^2+z^2)/(7-x)^2-1;

figure(Color='#113');

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;

X(274)=-.1;

Y(274)=1;

a=fill(X,Y,0);

hold on;

b=fill3(X,X-X,Y,0);

c=fimplicit3(f,EdgeC='n');

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.

set([a,b,c],FaceC='#CDE')

axis equal off;

view(-145,15);

light(Pos=[1 1 1]);

material metal;

xlim([-.9 1]);