OPTION BASE 1
DEFDBL A-H, J-W: DEFSNG X-Z: DEFINT I
DIM y(81): DIM x(81): DIM z(81): DIM scot(81)
y(1) = .00003959#: x(1) = .001379: z(1) = .0065: scot(1) = .000589
y(2) = .00006395#: x(2) = .002225: z(2) = .0105: scot(2) = .001108
y(3) = .00012: x(3) = .00425: z(3) = .0201: scot(3) = .002209
y(4) = .000216: x(4) = .00765: z(4) = .0362: scot(4) = .00453
y(5) = .000397: x(5) = .01431: z(5) = .0679: scot(5) = .00929
y(6) = .000643: x(6) = .02318: z(6) = .1102: scot(6) = .01852
y(7) = .00121: x(7) = .04352: z(7) = .2074: scot(7) = .03484
y(8) = .00218: x(8) = .07763: z(8) = .3713: scot(8) = .0604
y(9) = .004: x(9) = .1344: z(9) = .6456: scot(9) = .0966
y(10) = .0073: x(10) = .2148: z(10) = 1.0391: scot(10) = .1436
y(11) = .0116: x(11) = .2839: z(11) = 1.3858: scot(11) = .1998
y(12) = .0168: x(12) = .3285: z(12) = 1.623: scot(12) = .2625
y(13) = .023: x(13) = .3483: z(13) = 1.7471: scot(13) = .3281
y(14) = .0298: x(14) = .3481: z(14) = 1.7826: scot(14) = .3931
y(15) = .038: x(15) = .3362: z(15) = 1.7721: scot(15) = .455
y(16) = .048: x(16) = .3187: z(16) = 1.7441: scot(16) = .513
y(17) = .06: x(17) = .2908: z(17) = 1.6692: scot(17) = .567
y(18) = .0739: x(18) = .2511: z(18) = 1.5281: scot(18) = .62
y(19) = .091: x(19) = .1954: z(19) = 1.2876: scot(19) = .676
y(20) = .1126: x(20) = .1421: z(20) = 1.0419: scot(20) = .734
y(21) = .139: x(21) = .0956: z(21) = .813: scot(21) = .793
y(22) = .1693: x(22) = .058: z(22) = .6162: scot(22) = .851
y(23) = .208: x(23) = .032: z(23) = .4652: scot(23) = .904
y(24) = .2586: x(24) = .0147: z(24) = .3533: scot(24) = .949
y(25) = .323: x(25) = .0049: z(25) = .272: scot(25) = .982
y(26) = .4073: x(26) = .0024: z(26) = .2123: scot(26) = .998
y(27) = .503: x(27) = .0093: z(27) = .1582: scot(27) = .997
y(28) = .6082: x(28) = .0291: z(28) = .1117: scot(28) = .975
y(29) = .71: x(29) = .0633: z(29) = .0782: scot(29) = .935
y(30) = .7932: x(30) = .1096: z(30) = .0573: scot(30) = .88
y(31) = .862: x(31) = .1655: z(31) = .0422: scot(31) = .811
y(32) = .9149: x(32) = .2257: z(32) = .0298: scot(32) = .733
y(33) = .954: x(33) = .2904: z(33) = .0203: scot(33) = .65
y(34) = .9803: x(34) = .3597: z(34) = .0134: scot(34) = .564
y(35) = .995: x(35) = .4334: z(35) = .0087: scot(35) = .481
y(36) = 1.0002: x(36) = .5121: z(36) = .0057: scot(36) = .402
y(37) = .995: x(37) = .5945: z(37) = .0039: scot(37) = .3288
y(38) = .9786: x(38) = .6784: z(38) = .0027: scot(38) = .2639
y(39) = .952: x(39) = .7621: z(39) = .0021: scot(39) = .2076
y(40) = .9154: x(40) = .8425: z(40) = .0018: scot(40) = .1602
y(41) = .87: x(41) = .9163: z(41) = .0017: scot(41) = .1212
y(42) = .8163: x(42) = .9786: z(42) = .0014: scot(42) = .0899
y(43) = .757: x(43) = 1.0263: z(43) = .0011: scot(43) = .0655
y(44) = .6949: x(44) = 1.0567: z(44) = .001: scot(44) = .0469
y(45) = .631: x(45) = 1.0622: z(45) = .0008: scot(45) = .03315
y(46) = .5668: x(46) = 1.0456: z(46) = .0006: scot(46) = .02312
y(47) = .503: x(47) = 1.0026: z(47) = .00033: scot(47) = .01593
y(48) = .4412: x(48) = .9384: z(48) = .00023: scot(48) = .01008
y(49) = .381: x(49) = .8544: z(49) = .00017: scot(49) = .00737
y(50) = .321: x(50) = .7514: z(50) = .0001: scot(50) = .00497
y(51) = .26499: x(51) = .6424: z(51) = .00004: scot(51) = .003335
y(52) = .21699: x(52) = .5419: z(52) = .00002: scot(52) = .002235
y(53) = .17499: x(53) = .4479: z(53) = .00001: scot(53) = .001497
y(54) = .13823: x(54) = .3608: scot(54) = .001005
y(55) = .107: x(55) = .2835: scot(55) = .000677
y(56) = .08159: x(56) = .2187: scot(56) = .000459
y(57) = .06099: x(57) = .1649: scot(57) = .0003129
y(58) = .04458: x(58) = .1212: scot(58) = .0002146
y(59) = .032: x(59) = .0874: scot(59) = .000148
y(60) = .0232: x(60) = .0636: scot(60) = .0001026
y(61) = .01701: x(61) = .0468: scot(61) = .0000715
y(62) = .01192: x(62) = .0329: scot(62) = .0000501
y(63) = .008209: x(63) = .0227: scot(63) = .00003533#
y(64) = .005708: x(64) = .0158: scot(64) = .00002501#
y(65) = .004116: scot(65) = .0000178
y(66) = .002925: scot(66) = .00001273#
y(67) = .002095: scot(67) = .00000914#
y(68) = .00148: scot(68) = .0000066
y(69) = .001045: scot(69) = .00000478#
y(70) = .00072: scot(70) = .000003482#
y(71) = .0005: scot(71) = .000002546#
y(72) = .00036: scot(72) = .00000187#
y(73) = .00025: scot(73) = .000001379#
y(74) = .000175: scot(74) = .000001022#
y(75) = .00012: scot(75) = .00000076#
y(76) = .000083: scot(76) = .000000567#
y(77) = .00006: scot(77) = .000000425#
y(78) = .000041: scot(78) = .0000003196#
y(79) = .00003: scot(79) = .0000002413#
y(80) = .000021: scot(80) = .0000001829#
y(81) = .000014: scot(81) = .000000139#

FOR i = 54 TO 81: z(i) = 0: NEXT i
FOR i = 65 TO 81: x(i) = y(i) * 2.7689: NEXT i
10 SCREEN 0: CLS : INPUT "TEMP. K ( 350 TO 2E7 )"; t
IF t < 350 THEN 10
IF t > 20000000# THEN 10
ct = .01438 / t
l = 0: b = 0: r = 0: sco = 0
FOR i = 1 TO 81
u = .000000375# + .000000005# * i
lum = 1.27244D-19 / u / u / u / u / u / (EXP(ct / u) - 1)
l = l + y(i) * lum
b = b + z(i) * lum
r = r + x(i) * lum
sco = sco + scot(i) * lum
NEXT i
zpc = l * .00000031831#: PRINT : PRINT "CANDLES / cm^2 ="; zpc;
scc = sco * .00000079228#
zpsc = scc / zpc: PRINT "                    s/p ratio ="; zpsc
PRINT : PRINT "PEAK WAVELENGTH ="; USING "#####.##"; 2897760 / t;
PRINT " nm"
tt = r + l + b
zp = r / tt: PRINT : PRINT "X= "; USING "#.####"; zp; : PRINT "    ";
zp = l / tt: PRINT "Y= "; USING "#.####"; zp; : PRINT "    ";
zp = b / tt: PRINT "Z= "; USING "#.####"; zp
pt = .000000056696# * t * t * t * t: zp = pt
ct = .01438384# / t
PRINT : PRINT "TOTAL POWER="; zp * .0001; "W / cm^2"
p = 0
FOR u = .00000001# TO .00000020002# STEP .00000000005#
py = 1.86857D-24 / u / u / u / u / u / (EXP(ct / u) - 1)
p = p + py
IF u < .00000001002# THEN p = p - .5 * py
IF u > .00000019998# THEN p = p - .5 * py
NEXT u
PRINT : zp = p / pt: PRINT "% vacuum UV( 10 -  200 nm) ="; zp
p = 0
FOR u = .0000002 TO .00000028602# STEP .00000000005#
py = 1.86857D-24 / u / u / u / u / u / (EXP(ct / u) - 1)
p = p + py
IF u < .00000020002# THEN p = p - .5 * py
IF u > .00000028598# THEN p = p - .5 * py
NEXT u
zp = p / pt: PRINT "% UV-C    ( 200 -  286 nm) ="; zp
p = 0
FOR u = .000000286# TO .00000032002# STEP .00000000005#
py = 1.86857D-24 / u / u / u / u / u / (EXP(ct / u) - 1)
p = p + py
IF u < .00000028602# THEN p = p - .5 * py
IF u > .00000031998# THEN p = p - .5 * py
NEXT u
zp = p / pt: PRINT "% UV-B    ( 286 -  320 nm) ="; zp
p = 0
FOR u = .00000032# TO .00000040002# STEP .00000000005#
py = 1.86857D-24 / u / u / u / u / u / (EXP(ct / u) - 1)
p = p + py
IF u < .00000032002# THEN p = p - .5 * py
IF u > .00000039998# THEN p = p - .5 * py
NEXT u
zp = p / pt: PRINT "% UV-A    ( 320 -  400 nm) ="; zp
p = 0
FOR u = .0000004 TO .00000070002# STEP .00000000005#
py = 1.86857D-24 / u / u / u / u / u / (EXP(ct / u) - 1)
p = p + py
IF u < .00000040002# THEN p = p - .5 * py
IF u > .00000069998# THEN p = p - .5 * py
NEXT u
zp = p / pt: PRINT "% visible ( 400 -  700 nm) ="; zp
p = 0
FOR u = .0000007 TO .00000150002# STEP .00000000005#
py = 1.86857D-24 / u / u / u / u / u / (EXP(ct / u) - 1)
p = p + py
IF u < .00000070002# THEN p = p - .5 * py
IF u > .00000149998# THEN p = p - .5 * py
NEXT u
zp = p / pt: PRINT "% IR-A    ( 700 - 1500 nm) ="; zp
p = 0
FOR u = .0000015 TO .0000070002# STEP .0000000005#
py = 1.86857D-23 / u / u / u / u / u / (EXP(ct / u) - 1)
p = p + py
IF u < .0000015002# THEN p = p - .5 * py
IF u > .0000069998# THEN p = p - .5 * py
NEXT u
zp = p / pt: PRINT "% IR-B    (1500 - 7000 nm) ="; zp
p = 0
FOR u = .000007 TO .00100002# STEP .00000005#
py = 1.86857D-21 / u / u / u / u / u / (EXP(ct / u) - 1)
p = p + py
IF u < .00000702# THEN p = p - .5 * py
IF u > .00099998# THEN p = p - .5 * py
NEXT u
zp = p / pt: PRINT "% IR-C    (.007 -  1.0 mm) ="; zp

zp = 31415.92 * zpc / pt
PRINT : PRINT "LUMENS PER WATT = "; zp
PRINT : PRINT "hit X to end,     any other key to try another temp."
30 a$ = INKEY$: IF a$ = "" THEN 30
a$ = UCASE$(a$)
IF (a$) = "X" THEN END ELSE 10