Найти тему
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)
Около минуты