Расчёт космических скоростей для поверхности планеты (HP-35s)

Порт программы для fx-9750G. Исходные данные - радиус планеты в метрах и средняя плотность в кг на м3.
Вычисляется ускорение свободного падения на поверхности, а так же 1-я и 2-я космические скорости. Всё размерности в системе СИ.
Ввод текстовых сообщений после установки флага 10: EQN RCL R RCL A RCL D . . . ENTER

K001 LBL K
K002 SF 10
K003 RADIUS(M)
K004 PSE
K005 INPUT R
K006 P(KG/M^3)
K007 PSE
K008 INPUT P
K009 4
K010 x
K011 Pi
K012 x
K013 RCL R
K014 3
K015 y^x
K016 x
K017 3
K018 /
K019 G(from Constants)
K020 x
K021 RCL R
K022 x^2
K023 /
K024 STO G
K025 VIEW G
K026 X=V1,Y=V2
K027 PSE
K028 RCL R
K029 x
K030 sqrt()
K031 RCL G
K032 RCL R
K033 x
K034 2
K035 x
K036 sqrt(x)
K037 x<->y
K038 CF 10
K039 STOP
K040 GTO K001
K041 RTN
---------
CK=3CD0
LN=154

Тестовый пример:

XEQ K ENTER

RADIUS(M)
R?

6378000 R/S

P(KG/M^3)
P?

5510 R/S

G =
                9.8230

R/S

X=V1,Y=V2
11193.8679
7915.2599