- /
-
Share you an overcoat for keeping warm
on 24 Oct 2022
- 6
- 16
- 2
- 0
- 279
t = 0:3:360;
c = cosd(t).^3;
s = sind(t).^3;
r1 = 1:0.25:10;
r2=1:0.1:1;
r3=3.5:0.25:4;
f=@(r,l)bsxfun(@times,r.',l-1);
a = f(r1,c);
b = f(r1,s);
d = f(r2,c);
e = f(r2,s);
g = f(r3,c);
h = f(r3,s);
plot(a,b,'ro','MarkerFaceColor','y')
hold on
plot(d,e,'ko','MarkerFaceColor','k')
plot(g,h,'bo')
axis equal off