Phasenraumkurven einer Perle an einem Sinusförmigen Draht

Eine Kugel bewegt sich im einem Potential der Form V=m g A Sin[kρ].

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

Die Energie hat für positive A dann die Form

e = (p^2) + (l/(ρ^2)) + 10  Sin[ρ];

Die Phasenraumkurven für einen Drehimpuls von l=5 ergeben sich aus den Niveaulinien der Energie

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