t=[0:0.02:1];
n=size(t,2);
R=100;
eps=20;
h=pi/(2*(n-1));
y=zeros(n,1);
y1=zeros(n,1);
x=zeros(n,1);
x1=zeros(n,1);
y2=zeros(n,1);
for i=1:n
  x(i)=R*cos((i-1)*h);
  y(i)=R*sin((i-1)*h);
  y1(i)=0;
  x1(i)=eps+(i-1)*(R-eps)/(n-1);
end;
plot(x,y,(eps/R)*x,(eps/R)*y,x1,y1,y1,x1,'k');