- /
-
Hot air balloon
on 22 Oct 2021
- 42
- 192
- 1
- 0
- 262
l=@linspace;
p=@patch;
t=l(-pi,pi);
r=1+.37*sech(2.7*(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',[.3 .5 .8])
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,:),'y')
axis equal off