Stoke-em LawΒΆ

import math
rhos = 6000 #kg/m^3
rhoo = 850 # kg/m^3
gee = 9.8
dsp = 0.012
Cd = 5.0
muo = 0.1
velocity = math.sqrt(  ((rhos*gee-rhoo*gee)*(4/3)*dsp)/(Cd*rhoo))
print(velocity)
re=velocity*dsp/muo
print(re)
0.4358952923379111
0.052307435080549324