The Planck function radiance is B_lambda (T) = [2hc^2/lambda^5] [1/e^(hc/lambda k_b T) - 1] in units of Js^-1 m^-2 sr^-1 m^-1, where the wavelength lambda is in meters [m], the temperature T is in Kelvin [K], c = 2.998 times 10^8 ms^-1 is the speed of light, h = 6.626 times 10^-34 Js is Planck's constant, and k_b = 1.38 times 10^-23 J K^-1 is Boltzmann's constant. In terms of frequency the Planck function radiance is B_v (T) = [2hv^2/c^2] [1/e^(hv/k_b T) - 1] in units of J s^-1 m^-2 sr^-1 Hz^-1, where the frequency v is in Hertz [Hz]. a) Write a computer program in Matlab (or some other appropriate data analysis software package) called PlanckFcnWavelength_LastnameFirstname.m that computes the Planck function in 0.01 mu m steps from 0.01 mu m to 1000.00 mu m for a temperature of T_Earth = 255 K and a temperature T_Sun = 5917 K. Write all of your results to a file such that the first column contains index values that go from 1 to 100000, the second column contains the 100000 wavelengths in units of meters, the third column contains the 100000 solar radiance values in the units above and the fourth column contains the 100000 terrestrial radiance values in the units above. For columns 2 through 4 make sure your output is in exponential notation with 16 values and 8 decimal points. Call the data file containing these four columns of information PlanckFcnWavelength_EarthSun_data. Make a single plot with two curves, one curve of the solar radiance versus wavelength and the other of the terrestrial radiance versus wavelength