c program to generate ellipse
         open(unit=10,file='ellipse.dat',status='unknown')
         x=0.0
         do i=1,100
         x=x+0.01
         y=2.0*(1.0-x**2/1.0)**0.5
         write(10,*)x,y
         enddo
         stop
         end