Clearscreen.
//ER = Error
//PT = Previous time
//TS = Time step
//VP = vessel position (in a global vector mode)
//TL = Terrain level = geolocation under predicted vessel position
If SHIP:PERIAPSIS < 100 {

Set ER to 1000.
Set PT to Time:Seconds.
Set TS to 1.

Until ER < 100 {
//Set VP to vessel position at PT + TS.
Set VP to POSITIONAT(SHIP,(PT + TS)).

//Set TL to VP:geolocation.
set TL to BODY:GEOPOSITIONOF(VP):TerrainHeight.
Set ER to ((VP - BODY:POSITION):MAG - BODY:RADIUS - TL).
Set PT to PT + TS.
Print ER. //For testing
Print PT + TS.
Clearscreen.
}.
Print (PT + TS) - time:Seconds.
}.