program main
  implicit none
  
  call tabelle (-20, 20, 1);
end program

subroutine tabelle ( start, endw, step )
  implicit none
  integer start, endw, step
  integer k 
  
  if (step .LT. 1) then
    step = 1
  endif

  ! Tabellenkopf
  write (*,*) "+-------------------+"
  write (*,*) "|    x    |    y    |"
  write (*,*) "+-------------------+"
  
  ! x-y-Wertepaare ausgeben
  do k=start, endw, step
    write (*, '(1X,2H| ,I7,3H | , I7, 2H |)') k, (2*(k**2) + 4*k + 3)
  enddo
  
  write (*,*) "+-------------------+"
end subroutine