t=[0:0.02:1];
n=size(t,2);
h=1/(n-1);
R=100;
hx=2*R/(n-1);
ksi=80;
hy=ksi/(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+(i-1)*hx;
  y(i)=0;
  y1(i)=ksi;
  x1(i)=R;
  y2(i)=(i-1)*hy;
end;
plot(x,y,x,y1,x1,y2,-x1,y2,'k');