research/TCS/preliminary/trajectory/trajectory.f90

   1 ! Solve laplace's equation for a point potential at the origin
   2
   3 program main
   4
   5
   6         implicit none
   7
   8         real(8), allocatable :: V(:,:)
   9         real(8) :: width, dx, V0
  10         integer :: i, j
  11         i = ceiling(width/dx)
  12         allocate(V(-i:i, -i:i))
  13
  14         V(0,0) = V0
  15
  16
  17
  18 end program main
  19
  20