Yukterez hat geschrieben:
Nun, ich würde mal sagen im Gegensatz zu deinem sieht es zumindest aus. Man könnte es in etwa so beschreiben: man gibt Zeit, Startwinkel, Startgeschwindigkeit, Radius und Erdanziehung hinein und bekommt Position und Geschwindigkeit heraus.

,

Da deine "explizite" Lösung von Wolfram auch keine geschlossene Lösung ist (die es auch gar nicht gibt), und daher nur auf numerische Näherungswerte führt, kann man sich die ganze Wolframabschreiberei sparen und gleich eine Näherungslösung aus der Differentialgleichung verwenden. Für phi=0 als oberer Totpunkt.
 = 0)
Läßt sich einmal integrieren zu
 = C)
 = C)
})
Wird numerisch gelöst durch Ansatz von Startwerten und dann werden für konstante dt aktuelle winkel und omega berechnet.
ich hoffe, du besitzt noch einen Basic Compiler:
Code:
dim as single g,r ,C,omeganull,omega,phi,t,deltat,phialt
omeganull=0.004:g=9.81:r=0.981:phi=0:t=0:deltat=0.05
C=omeganull^2/2 + g/r
omega=omeganull
print " t","phi","omega","delta phi"
print
do
print t,phi,omega,phi-phialt
sleep 100
omega=sqr(2*C-cos(phi)*2*g/r)
phialt=phi
phi=phi+omega*deltat
t=t+deltat
loop until phi>2*3.1415
sleep
end
.
.