def test_b_wien():
 from .. import b_wien
 from ... import units as u
 t = (5778 * u.K)
 w = (b_wien / t).to(u.nm)
 assert (round(w.value) == 502)

def test_b_wien ( ) : from . . import b_wien from ... import units as u t = ( 5778 * u . K ) w = ( b_wien / t ) . to ( u . nm ) assert ( round ( w . value ) == 502 )
b_wien should give the correct peak wavelength for given blackbody temperature .