- /
-
Shhhhiii,......
on 10 Oct 2022
- 2
- 5
- 0
- 0
- 274
t = (0:0.001:2*pi)';
n = 0.1:0.015:0.75;
r = (abs(cos(2*t/4)).^n+abs(sin(2*t/4)).^n).^(-1./n);
x1 = 3.85*r.*cos(t);
y1 = 2*r.*sin(t);
h=plot(x1,y1,'r')
rotate(h,[0,90],-5)
axis equal off
hold on
patch(-0.7*y1(:,25),0.1*x1(:,25),'k','LineW',11)
xp = -1.2:0.1:1.2;
yp1 = xp.^4-1;
yp2 = xp.^6-1.5;
%patch([xp,flip(xp)],[yp1,flip(yp2)],'k')