import numpy as np
from sympy.interactive import init_printing

init_printing()

w = 1.4 # kW/m^2
c = 3e8 # m/s

M_He4 = 4.0039 # gram/mol
M_H   = 1.0080 # gram/mol

dM = 4 * M_H - M_He4 # gram/mol
print('dM = {0:.4f} gram/mol'.format(dM))

# distance Sun-Earth
R = 1.5e11               # meters
W = 4 * np.pi * R**2 * w # kW
q = W/c**2               # tonnes/s
print('q = {:e} tonnes'.format(q))

q_H = 4 * q * M_H / dM # tonnes/s
print('result: {:e} tonnes/s'.format(q_H))

710deebc1d28cd8460ad6b196d94b3a5fc633595
<mask token>

<mask token>
init_printing()
<mask token>
print('dM = {0:.4f} gram/mol'.format(dM))
<mask token>
print('q = {:e} tonnes'.format(q))
<mask token>
print('result: {:e} tonnes/s'.format(q_H))

<mask token>
init_printing()
w = 1.4
c = 300000000.0
M_He4 = 4.0039
M_H = 1.008
dM = 4 * M_H - M_He4
print('dM = {0:.4f} gram/mol'.format(dM))
R = 150000000000.0
W = 4 * np.pi * R ** 2 * w
q = W / c ** 2
print('q = {:e} tonnes'.format(q))
q_H = 4 * q * M_H / dM
print('result: {:e} tonnes/s'.format(q_H))

import numpy as np
from sympy.interactive import init_printing
init_printing()
w = 1.4
c = 300000000.0
M_He4 = 4.0039
M_H = 1.008
dM = 4 * M_H - M_He4
print('dM = {0:.4f} gram/mol'.format(dM))
R = 150000000000.0
W = 4 * np.pi * R ** 2 * w
q = W / c ** 2
print('q = {:e} tonnes'.format(q))
q_H = 4 * q * M_H / dM
print('result: {:e} tonnes/s'.format(q_H))

import numpy as np
from sympy.interactive import init_printing

init_printing()

w = 1.4 # kW/m^2
c = 3e8 # m/s

M_He4 = 4.0039 # gram/mol
M_H   = 1.0080 # gram/mol

dM = 4 * M_H - M_He4 # gram/mol
print('dM = {0:.4f} gram/mol'.format(dM))

# distance Sun-Earth
R = 1.5e11               # meters
W = 4 * np.pi * R**2 * w # kW
q = W/c**2               # tonnes/s
print('q = {:e} tonnes'.format(q))

q_H = 4 * q * M_H / dM # tonnes/s
print('result: {:e} tonnes/s'.format(q_H))
