#### code:

```
xx1 = [...];
yy1 = [...];
rr1 = [...];
theta=0:pi/50:2*pi;
for i = 1:length(xx1)
X = xx1(i)+rr1(i)*cos(theta);
Y = yy1(i)+rr1(i)*sin(theta);
plot(X,Y,['k','-'],xx1(i),yy1(i),['k','.']);
hold on;
axis([0 4 0 4])
axis square;
end
```

```
xx1 = [...];
yy1 = [...];
rr1 = [...];
theta=0:pi/50:2*pi;
for i = 1:length(xx1)
X = xx1(i)+rr1(i)*cos(theta);
Y = yy1(i)+rr1(i)*sin(theta);
plot(X,Y,['k','-'],xx1(i),yy1(i),['k','.']);
hold on;
axis([0 4 0 4])
axis square;
end
```

## Comments