;
pro prob08c
common state, state
sav_filename = 'prob08a_060_.sav'
restore, filename = sav_filename, /verbose
state = {S: 0.995d, kohler_a: kohler_a, kohler_b: kohler_b}
result = newton(rd * 4, 'droplet_radius', tolx = 1e-10, tolf = 1e-6, /double)
print, result, (-state.kohler_b/(state.s - 1))^(1./3.)
end
;
function droplet_radius, x
common state, state
return, state.S - 1.d - state.kohler_a / x + state.kohler_b / x^3
end