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

Исходные данные - радиус планеты в метрах и средняя плотность в кг на м3.
Вычисляется ускорение свободного падения на поверхности, а так же 1-я и 2-я космические скорости. Всё размерности в системе СИ.

File attachments: 
AttachmentSize
File KOCMCKOP.RAR344 bytes

Comments

Контрольный пример. А можно, в лучших традициях ПМК, сопровождать программу контрольным примером? Чтобы проверять при портировании на другие калькуляторы.

---------------------------
Истина где-то рядом
www.litres.ru/vitaliy-samurov/dozvonitsya-do-devy/

А чем плоха Земля в качестве контрольного примера? Исходные данные
Радиус: 6378000 м
Плотность: 5510 кг/м3
Результат
Ускорение свободного падения = 9.818622367
1-я космическая = 7913.480489 м/с
2-я космическая = 11191.35143

Мои программируемые калькуляторы:
Б3-21, Б3-34, МК-61, МК-52, МК-85
CASIO: cfx-9850GB+, fx-9750G+, fx-9750GII, fx-9860G, Algebra fx-2.0, fx-5800P, fx-7400G+
HP: 50G, 48G, 35s
TI: Nspire-CAS, Voyage-200, 89Titanium
SHARP EL-9600G

Та же программа для HP-35s. Ничем не плоха. Спасибо :)

Та же программа для HP-35s. Не совсем понятно, откуда такое расхождение с Вашей программой и с Python-ом:

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 RTN
---------
CK=1FD2
LN=151

Тот же пример (в калькуляторе установлено округление до 4-х знаков после точки):
XEQ K
R?
6378000 R/S
5510 R/S
G = 9.8230
R/S
11193.8679
7915.2599

Python:

>>> from math import*
>>> R=6378000
>>> R
6378000
>>> P=5510
>>> P
5510
>>> M=P*4*3.1415*R**3/3
>>> M
5.987988315636265e+24
>>> G=6.67e-11*M/(R**2)
>>> G
9.818332790371999
>>> sqrt(G*R)
7913.363794050708
>>> sqrt(2*G*R)
11191.186401538722
>>> 

---------------------------
Истина где-то рядом
www.litres.ru/vitaliy-samurov/dozvonitsya-do-devy/

Расхождение, вероятно связано с точностью констант. Я забивал скорость света с точностью до метра, ну и в гравитационной постояной больше знаков мантиссы, чем 3.

Мои программируемые калькуляторы:
Б3-21, Б3-34, МК-61, МК-52, МК-85
CASIO: cfx-9850GB+, fx-9750G+, fx-9750GII, fx-9860G, Algebra fx-2.0, fx-5800P, fx-7400G+
HP: 50G, 48G, 35s
TI: Nspire-CAS, Voyage-200, 89Titanium
SHARP EL-9600G

Да, в константах дело. Ввёл в программу 6.67e-11, выдало ускорение св.падения в 9.818622

---------------------------
Истина где-то рядом
www.litres.ru/vitaliy-samurov/dozvonitsya-do-devy/

Формулы. Было бы просто отлично в самых лучших традициях приводить математику и, если требуется, блок-схему :)

+1 :)

---------------------------
Истина где-то рядом
www.litres.ru/vitaliy-samurov/dozvonitsya-do-devy/