program series
   integer :: n
   real :: sum_series       ! sum is a keyword
   sum_series = 0.0
!  add the first 100 terms of a simple series
   do n = 1, 100
      sum_series = sum_series + 1.0/real(n)**2
      print *, n,sum_series
   end do
end program series