potential1ds.f90

Ground state of the time-independent Gross-Pitaevskii equation

real(kind(1.d0)) function potentialV0(z)
  !
  ! User supplied additional potential
  !
  ! This function takes on input the coordinates z and returns the value
  ! of the potential at this point *** exclusive of the harmonic trap ***
  !
  implicit none
  integer, parameter :: dp=kind(1.d0)
  
  real(dp), intent(in) :: z
  
  real(dp), parameter :: gamma = 2.d0

  ! quartic potential
  potentialV0 = gamma*z**4
 
  return
end function potentialV0