clear all;
t=[0:.005:.5];
n=length(t);
a=1.5;
x1=linspace(-1,1);
for i=1:n
   tt=t(i);
   if tt==0 
     f1=sqrt(1-x1.^2);
   else
     f1=1/sqrt(2)/tt*real(sqrt(sqrt((1-2*tt*x1).^2+4*tt^2*(1-x1.^2))-1+2*tt*x1));
   end
   fill(f1,-x1,'b')
   axis('equal');
   hold on
   plot([0,0],[-1.5,1.5])
   hold off
   pause(0.001);
end;