Menu "Macros"
	"BBRadiation"
End

Function BBRadiation()
	String temperature="1273", graphName, xWave, yWave
	Variable upLim=10000, tempNum
	
	Prompt temperature, "Temperature for Black Body Radiation in Kelvin"
	Prompt upLim, "Upper limit for Main Graph Wavelength in nm"
	DoPrompt "Enter Parameters", temperature, upLim
	
	
	xWave = "BBWaveLength" + temperature
	yWave = "BBIntensity" + temperature
	Make/D/N=(upLim) $xWave, $yWave
	Wave xWaveWav = $xWave
	Wave yWaveWav = $yWave
	tempNum = str2num(temperature)
	yWaveWav=(((2*pi*(299792458^2)*6.62606957E-34)/((x*(1E-9))^5))*(1/(e^(((6.62606957E-34)*299792458)/((x*(1E-9))*(1.3806488E-23)*tempNum))-1)))
	xWaveWav=x
	graphName = "BBRadiation_"+temperature
	Display/N=graphName yWaveWav vs xWaveWav
	Label left "Intensity (watts/m\\S2\\M)";DelayUpdate
	Label bottom "Wavelength (nm)"
	
	
End