%%HP: T(3)A(D)F(.);
@ PLNK by Nelson W. Wallace
DIR
  BBOD2
    \<< "BBOD2"
"T in K, \Gl's in \Gmm"
{
"  :T: 
 :\Gl1: 
 :\Gl2: "
{ 1 0 } } INPUT
OBJ\-> DTAG '1_\Gm'
\->UNIT '\Gl2' STO DTAG
'1_\Gm' \->UNIT '\Gl1'
STO DTAG '1_K'
\->UNIT 'T' STO 4 FIX
BBLST TMENU CLLCD
"" "" ""
"\.S W  / W / \.S dWdT"
    \>>
  PLOT
    \<< "PLOT" 'WW'
RCL STEQ 1 35 XRNG
'\Gl1' INDEP ERASE
AUTO DRAW PICT RCL
\-> wwp
      \<< 'WP' RCL
STEQ ERASE AUTO
DRAW PICT RCL \-> wpp
        \<< wwp {
# 0d # 0d } wpp GOR
PICT STO 0 FIX {
(1,0) "WAVELENGTH"
"BLACKBODY" } AXES
LABEL DRAX { }
PVIEW
        \>>
      \>>
    \>>
  WW 'C1/(\Gl1^5*(EXP
(C2/(\Gl1*T))-1))'
  WP 'C3/(\Gl1^4*(EXP
(C2/(\Gl1*T))-1))'
  \Gl '1_\Gm'
  T '300_K'
  BBLST { { "WATTS"
    \<< '\.S(\Gl1,\Gl2,C1/(
\Gl^5*(EXP(C2/(\Gl*T))-
1)),\Gl)' \->NUM 6 FIX
    \>> } { "PHOTN1"
    \<< '\.S(\Gl1,\Gl2,C3/(
\Gl^4*(EXP(C2/(\Gl*T))-
1)),\Gl)' \->NUM 4 SCI
    \>> } { "Watts2 "
    \<< 'C1/(\Gl1^5*(
EXP(C2/(\Gl1*T))-1))'
EVAL
    \>> } { "Photn2 "
    \<< 'C3/(\Gl1^4*(
EXP(C2/(\Gl1*T))-1))'
EVAL
    \>> } { "Watts3"
    \<< '\.S(\Gl1,\Gl2,C1*
C2*EXP(C2/(\Gl*T))/(\Gl
^6*T^2*(EXP(C2/(\Gl*T
))-1)^2),\Gl)' \->NUM
    \>> } { "Photn3"
    \<< '\.S(\Gl1,\Gl2,C3*
C2*EXP(C2/(\Gl*T))/(\Gl
^5*T^2*(EXP(C2/(\Gl*T
))-1)^2),\Gl)' \->NUM
    \>> } }
  C1 '37415_W*\Gm^4/
cm^2'
  C2 '14387.9_\Gm*K'
  C3 '1.88365E23_\Gm^
3/(s*cm^2)'
  \Gl1 '35.85_\Gm'
  \Gl2 '12_\Gm'
  PPAR {
(1,-2.55163531432E16)
(35,1.70109020955E17)
\Gl1 0 { (1,0)
"WAVELENGTH"
"BLACKBODY" }
FUNCTION Y }
  EQ 'C3/(\Gl1^4*(EXP
(C2/(\Gl1*T))-1))'
  IERR '
.641028178082_W/cm^
2'
END
