clc;
clear;
r=1;
alp=0:0.1:2*pi;
a=5;
b=5;
x=a+r*cos(alp);
y=b+r*sin(alp);
fill(x,y,'r');
axis([-10 30 -10 30]);