10 подписчиков
Исходный код:
m = 11*(5.9736*10**24)#Сколько звезда имеет солнечных масс
Rb = 617100000000#Радиус Бетельгейзе
M = (4+12+16)/2#Средняя молярная масса
P = (2 + 6 + 8)/2#Средний порядновый номер
N = (m/M*(6.02*10**23))*((P*2+P)*2)#Количество кварков
v = N * ((((10**(-19))/2)**3)*((4/3)*3.14))#Обьем N числа кварков в кубических метрах
G = 6.67430*10**(-11)#Гравитационная постоянная
_R = int((G*m)/(300000**2))#Радиус горизонта событий
_R *= 1000
p = m/v#Плотность, она тут не нужна
print(p)#Вывод плотности
R = ((3*v)/(4*3.14))*(1/3)#Радиус ядра
g = (G * m)/R#Ускорение свободного падения на поверхности ядра черной дыры
if R < _R and Rb > _R:
print('Радиус звезды Бетельгейзе, которую мы и брали в расчет равен: ' + str(Rb) + 'м' + '\nРадиус горизонта событий равен: ' + str(_R) + 'м' + '\nА радиус ядра черной дыры равен: ' + str(R) + 'м' + '\nУскорение свободного падения на поверхности равно: ' + str(g) + '\nГипотиза верна')
else:
print('Гипотеза не верна: радиус горизонта событий меньше радиуса ядра черной дыры')
print(_R)
print(R)
Около минуты
4 февраля 2022