Programming Midterm

program midterm
  implicit none
  real:: Ax=0.0, Ay=0.0, Mxl=0.0, Mxr=0.0, Myl=0.0, Myr=0.0, SumA=0.0, SumMx=0.0, SumMy=0.0,xc=0.0,yc=0.0
  integer:: i=0, N

   print*, "Inout N here"
   read(*,*) N

   do i= 1,N-1
        Ax = Ax + 1 - real(i)/real(N)
        Ay = Ay + sqrt(1-(- 1+real(i)/real(N))**2)
        Mxl = Mxl + 0.5*(1-(- 1+real(i)/real(N))**2)
        Mxr =  Mxr + 0.5*(1-real(i)/real(N))**2
        Myl = Myl + (- 1+real(i)/real(N))*sqrt(1-(- 1+real(i)/real(N))**2)
        Myr = Myr + real(i)/real(N)*(1-real(i)/real(N))
   end do

   SumA =  (Ax + Ay + 1.0)/real(N)
   SumMx = ( Mxl+ Mxr+ 0.5 )/real(N)
   SumMy = ( Myl+ Myr)/real(N)

   xc = SumMy/SumA
   yc = SumMx/SumA

   print *, "xc=",xc,"yc=",yc
end program midterm