! Solve laplace's equation for a point potential at the origin

program main

	
	implicit none

	real(8), allocatable :: V(:,:)
	real(8) :: width, dx, V0
	integer :: i, j
	i = ceiling(width/dx)
	allocate(V(-i:i, -i:i))

	V(0,0) = V0
	
	

end program main