I am trying to write a simple python code for a plot of intensity vs wavelength for a given temperature, T=200K. So far I have this...
import scipy as sp imp
Just want to point out that there seems to be an equivalent of what OP wants to do in astropy:
https://docs.astropy.org/en/stable/api/astropy.modeling.physical_models.BlackBody.html
Unfortunately, it is not very clear to me yet how to get wavelength vs frequency based expression.