effektives Gravitationspotential ploten

Invis · Anmeldungsdatum: 07.12.2013 Beiträge: 1

Hey Leute.

Ich habe ein kleines (zurzeit für mich ein Großes) Problem:
Ich versuche schon seit Tagen das effektive Potential in Matlab zu ploten, welches wir vorliegen haben, wenn wir uns den Merkur auf seiner Umlaufbahn um die Sonne anschauen.

Die allgemeine Definition des effektiven Garvitationspotentials ist gegeben durch:

mit

m = Masse des Merkurs, M = masse der Sonne, G = Gravitationskonstante, L Drehimpuls des Merkur.

Als erstes würde ich den Drehimpuls berechnen. Dazu habe ich mich an den Merkur-Daten http://nssdc.gsfc.nasa.gov/planetary/factsheet/mercuryfact.html bedient.
An dem Punkt der Periapsis gilt für die Bahngeschindigkeit mit

und

:

Daraus können wir nun mit Hilfe der maximalen Bahngeschwindigkeit

an der Periapsis

den Drehimpuls berechnen, der ja konstant ist:

Mein Problem ist jetzt, dass wenn ich nun das effetktive Potential ploten will, nicht den typischen Graphen bekomme, sondern solch ein "gedöns":
https://drive.google.com/file/d/0BymH6Rz7PTmmZTQzMnhtdHdubVU/edit?usp=sharing

Hier mal mein Matlab-Code:

jmd · Anmeldungsdatum: 28.10.2012 Beiträge: 577

unten steht a/r^2 ?

Außerdem
Die Einheiten derart zu verändern ist keine gute Idee
Warum nimmst du nicht die SI Einheiten und wandelst die erst ganz unten um?

x=r*1AE und V(x)