- /
-
Parachute
on 22 Oct 2021
- 37
- 72
- 0
- 0
- 263
l=@linspace;
p=@patch;
t=l(2*pi,-2*pi);
r=1+.3*sech(3*(pi/2+t));
x=r.*cos(t);
y=r.*sin(t);
a=[.1;.1];
b=-[1.3;1.7];
c=.1*[-1 -1 1 1];
d=-[1.7 2 2 1.7];
% figure('color', [.6 1 1])
hold
plot([a -a],b,'k')
p(c,d,'k')
p(x,y,'r')
X=[x;y];
P=[.8 0;0 1]*X;
p(P(1,:),P(2,:),'g')
P=[.4 0;0 1]*X;
p(P(1,:),P(2,:),'b')
axis equal off