Phasenraumkurven einer Kugel im Gugelhupf

Eine Kugel bewegt sich in einem Potential der From V=m g ((A(ρ - (Subscript[ρ, 0])))^2).  Für positives A erhält man

Clear[z]z[ρ_ /; ρ < 4] := (((ρ - 2))^2)z[ρ_ /; ρ ≥ 4] := (((4 - 2))^2)Plot3D[z[(Sqrt[(x^2) + (y^2)])], {x, -4, 4}, {y, -4, 4}, 
      PlotPoints -> 51, PlotRange -> All]

Clear[z]z[ρ_ /; ρ < 4] := -(((ρ - 2))^2)z[ρ_ /; ρ ≥ 4] := -(((4 - 2))^2)Plot3D[z[(Sqrt[(x^2) + (y^2)])], {x, -4, 4}, {y, -4, 4}, 
      PlotPoints -> 51, PlotRange -> All]

e = (p^2) + (l/(ρ^2)) + 10 (((ρ - 2))^2);

\(t = 
      Table[Evaluate[e /. l -> 5], {p, \(-9\), 9, .4}, {ρ, 0.3, 
          5, .15}]; 
ListPlot3D[t]\)ListContourPlot[t, Contours -> 40, ContourShading -> False]

e = (p^2) + (l/(ρ^2)) - 10 (((ρ - 2))^2);

\(t = 
      Table[Evaluate[e /. l -> 5], {p, \(-9\), 9, .4}, {ρ, 0.3, 
          5, .15}]; 
ListPlot3D[t]\)ListContourPlot[t, Contours -> 40, ContourShading -> False]