Is possible to simulate gravity?

[jesusmontera]

I have modified your sample 1D potential barrier with a negative potential to simulate gravity, but the particle with momentum does not get trapped. It reflects or cross the potential well. Is it possible to simulate gravity in 1D with your library?.

def GravitationalPotential(n,xcenter,minPot):
  pot= np.zeros(n,dtype=float)
  for x in range(n):
    d = abs(x-xcenter)
    if d ==0:
      pot[x]=minPot
    else:
      pot[x] = minPot / d**2
  return pot

V0 =-20
barrier = GravitationalPotential(128,64,V0)

[jesusmontera]

Oh my god, gravitational potential energy is r instead of r**2 from force. Now particle with m=1 AU and Kinetic 1 Mev can't move near a black hole (but wave oscilates). I suppose I must increase its Kinetic