%!PS-Adobe-3.0 EPSF-3.0 TG030-7N.EPS %%BoundingBox: 70 82 785 580 %START PDFDE011.EPS /pdfmark07 where {pop} {userdict /pdfmark07 /cleartomark load put} ifelse /languagelevel where {pop languagelevel} {1} ifelse 2 lt { userdict (<<) cvn ([) cvn load put userdict (>>) cvn (]) cvn load put} if [ /Title (PostScript pictures: http://130.149.60.45/~farbmetrik/TG03/) /Author (compare K. Richter "Computergrafik ...": ISBN 3-8007-1775-1) /Subject (goto: http://www.ps.bam.de or http://130.149.60.45/~farbmetrik) /Keywords (image reproduction, colour devices, colour vision) /Creator (klaus.richter@mac.com) /CreationDate (D:2013060112000) /ModDate (D:2013060112000) /DOCINFO pdfmark07 [ /View [ /FitB ] /DOCVIEW pdfmark07 %END PDFDE011 %SETCMYK_OLV.PS /setcmyk_olvcolor %transfer cmyk to olv {/ncolor exch def /ycolor exch def /mcolor exch def /ccolor exch def ncolor 0 eq {1 ccolor sub 1 mcolor sub 1 ycolor sub setrgbcolor } {ccolor 0 eq mcolor 0 eq and ycolor 0 eq and {%only black n 1 ncolor sub dup dup setrgbcolor } {%cmy and black n 1 ccolor sub 1 ncolor sub mul 1 mcolor sub 1 ncolor sub mul 1 ycolor sub 1 ncolor sub mul setrgbcolor } ifelse } ifelse } bind def /FS {findfont exch scalefont setfont} bind def /MM {72 25.4 div mul} def /str {8 string } bind def /languagelevel where {pop languagelevel} {1} ifelse /PSL12 exch def /dictende {counttomark 2 idiv dup dict begin {def} repeat pop currentdict end} bind def /Times-Roman findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse }forall /Encoding ISOLatin1Encoding def currentdict end /Times-ISOL1 exch definefont pop /Times-Italic findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse }forall /Encoding ISOLatin1Encoding def currentdict end /TimesI-ISOL1 exch definefont pop /Times-Bold findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse }forall /Encoding ISOLatin1Encoding def currentdict end /TimesB-ISOL1 exch definefont pop /Times-BoldItalic findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse }forall /Encoding ISOLatin1Encoding def currentdict end /TimesBI-ISOL1 exch definefont pop /Courier-Bold findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse }forall /Encoding ISOLatin1Encoding def currentdict end /CourierB-ISOL1 exch definefont pop /TS {170 /Times-ISOL1 FS} bind def /TK {200 /Times-ISOL1 FS} bind def /TM {270 /Times-ISOL1 FS} bind def /TG {320 /Times-ISOL1 FS} bind def /TIS {170 /TimesI-ISOL1 FS} bind def /TIK {200 /TimesI-ISOL1 FS} bind def /TIM {270 /TimesI-ISOL1 FS} bind def /TIG {320 /TimesI-ISOL1 FS} bind def /TBS {170 /TimesB-ISOL1 FS} bind def /TBK {200 /TimesB-ISOL1 FS} bind def /TBM {270 /TimesB-ISOL1 FS} bind def /TBG {320 /TimesB-ISOL1 FS} bind def /TBIS {170 /TimesBI-ISOL1 FS} bind def /TBIK {200 /TimesBI-ISOL1 FS} bind def /TBIM {270 /TimesBI-ISOL1 FS} bind def /TBIG {320 /TimesBI-ISOL1 FS} bind def /CBKK {200 /CourierB-ISOL1 FS} bind def /CBK {250 /CourierB-ISOL1 FS} bind def /CBM {300 /CourierB-ISOL1 FS} bind def /CBG {350 /CourierB-ISOL1 FS} bind def /TSS {170 /Symbol FS} bind def /TSK {200 /Symbol FS} bind def /TSM {270 /Symbol FS} bind def /TSG {320 /Symbol FS} bind def /F89 0.8859 def % achromatic (D65), (D50), (C00) colours /XLWD65100 95.04 def /YLWD65100 100.00 def /ZLWD65100 108.88 def /XLWD65089 95.04 F89 mul def /YLWD65089 100.00 F89 mul def /ZLWD65089 108.88 F89 mul def /XLWD50100 96.42 def /YLWD50100 100.00 def /ZLWD50100 82.51 def /XLWD50089 96.42 F89 mul def /YLWD50089 100.00 F89 mul def /ZLWD50089 82.51 F89 mul def /XLWC00100 94.94 def /YLWC00100 100.00 def /ZLWC00100 108.76 def /XLWC00089 94.94 F89 mul def /YLWC00089 100.00 F89 mul def /ZLWC00089 108.76 F89 mul def /XYZC00100 303 array def /LABC00100 303 array def /YABC00100 303 array def /XYZC00089 303 array def /LABC00089 303 array def /YABC00089 303 array def %'TK0X'NCS_0099_ALLCD.LIS %NCS HUE CIRCLE 80 STEPs FROM NCS_CIE.TXT FOR ILLUMINANT C %FROM SWEDISH STANDARD SS019104 (1998), ANGLE ORDER %EXTRACT OF FILE 'TK0X'NCS_CIE_04.LIS %ELEMENTARY CIRCLE R->Y->G->B->R /NCS_0099_ALLCD 909 array def %80x9 + 21x9 steps, XYZ, xy, Lab* /NCS_0099_ALLCD [%BEG NCS_0099_ALLCD D=DECENDING=ANGLE ORDER % X Y Z x y L* a* b* (0099-R00B) 27.17 11.60 1.39 0.6765 0.2889 40.57 85.51 50.80 %calculated L*a*b* values of NCS seem very strange, XYZ seem ok. (0099-Y95R) 29.76 13.21 0.35 0.6869 0.3049 43.08 159.10 101.25 (0099-Y90R) 32.26 15.00 0.00 0.6826 0.3174 45.63 156.61 115.28 (0099-Y85R) 34.42 16.88 0.00 0.6709 0.3291 48.11 151.48 123.16 (0099-Y80R) 36.29 18.80 0.02 0.6585 0.3411 50.45 144.92 130.31 (0099-Y75R) 38.16 20.74 0.09 0.6470 0.3516 52.66 138.87 136.31 (0099-Y70R) 40.00 22.70 0.00 0.6380 0.3620 54.76 133.09 144.62 (0099-Y65R) 41.70 24.67 0.00 0.6283 0.3717 56.75 126.98 151.10 (0099-Y60R) 43.23 26.70 0.00 0.6182 0.3818 58.70 119.97 157.45 (0099-Y55R) 44.64 28.85 0.05 0.6070 0.3923 60.65 111.95 162.93 (0099-Y50R) 46.17 31.20 0.07 0.5963 0.4029 62.68 103.79 169.37 (0099-Y45R) 48.06 33.87 0.07 0.5861 0.4130 64.86 96.02 176.49 (0099-Y40R) 50.43 37.00 0.05 0.5765 0.4230 67.28 88.51 184.89 (0099-Y35R) 53.27 40.69 0.09 0.5664 0.4326 69.96 80.76 192.97 (0099-Y30R) 56.32 44.80 0.08 0.5565 0.4427 72.76 72.43 202.58 (0099-Y25R) 59.37 49.31 0.12 0.5457 0.4532 75.64 62.90 211.48 (0099-Y20R) 63.00 54.80 0.08 0.5344 0.4649 78.93 52.07 223.06 (0099-Y15R) 67.28 61.55 0.17 0.5216 0.4771 82.67 39.54 234.08 (0099-Y10R) 70.67 68.00 0.31 0.5085 0.4893 86.01 25.93 243.03 (0099-Y05R) 71.81 72.44 0.50 0.4961 0.5004 88.18 12.38 247.27 (0099-Y00R) 70.77 74.50 0.70 0.4848 0.5104 89.16 -0.16 247.59 (0099-G95Y) 68.01 74.35 0.83 0.4750 0.5193 89.09 -11.22 245.52 (0099-G90Y) 64.67 73.10 0.84 0.4666 0.5274 88.50 -20.72 243.37 (0099-G85Y) 61.65 71.78 0.91 0.4589 0.5343 87.86 -28.94 240.26 (0099-G80Y) 58.80 70.30 0.83 0.4525 0.5411 87.14 -36.12 239.03 (0099-G75Y) 55.79 68.37 0.59 0.4472 0.5480 86.19 -42.49 239.21 (0099-G70Y) 52.61 66.00 0.57 0.4414 0.5538 85.00 -48.20 235.57 (0099-G65Y) 49.23 63.24 0.58 0.4355 0.5594 83.57 -53.65 230.55 (0099-G60Y) 45.54 60.00 0.65 0.4288 0.5650 81.84 -59.08 223.73 (0099-G55Y) 41.45 56.24 0.52 0.4220 0.5726 79.75 -64.94 218.67 (0099-G50Y) 37.13 52.20 0.42 0.4137 0.5816 77.40 -71.59 212.40 (0099-G45Y) 32.82 48.21 0.41 0.4030 0.5920 74.96 -79.51 204.41 (0099-G40Y) 28.62 44.50 0.37 0.3894 0.6056 72.56 -89.52 197.11 (0099-G35Y) 24.60 41.26 0.25 0.3722 0.6241 70.36 -102.50 191.72 (0099-G30Y) 20.89 38.50 0.32 0.3498 0.6449 68.39 -118.15 184.03 (0099-G25Y) 17.45 36.19 0.62 0.3216 0.6669 66.66 -136.75 173.57 (0099-G20Y) 13.89 34.30 0.61 0.2846 0.7030 65.20 -163.20 169.00 (0099-G15Y) 10.03 32.77 0.58 0.2313 0.7553 63.97 -202.12 165.30 (0099-G10Y) 6.69 31.50 2.15 0.1659 0.7808 62.93 -133.71 82.03 (0099-G05Y) 4.63 30.42 6.00 0.1128 0.7411 62.01 -153.65 58.41 (0099-G00Y) 3.71 29.50 11.65 0.0827 0.6575 61.22 -163.22 38.18 (0099-B95G) 3.73 28.75 18.43 0.0732 0.5647 60.56 -160.11 21.36 (0099-B90G) 4.59 28.20 24.53 0.0801 0.4919 60.07 -145.79 9.45 (0099-B85G) 5.24 27.88 28.51 0.0850 0.4524 59.78 -136.36 2.71 (0099-B80G) 5.41 27.70 31.32 0.0840 0.4299 59.62 -133.58 -1.65 (0099-B75G) 5.69 27.58 34.02 0.0846 0.4099 59.51 -129.82 -5.53 (0099-B70G) 6.20 27.50 36.72 0.0880 0.3905 59.43 -123.92 -9.16 (0099-B65G) 6.63 27.46 39.53 0.0900 0.3730 59.40 -119.19 -12.68 (0099-B60G) 7.10 27.40 42.36 0.0924 0.3565 59.34 -114.13 -16.10 (0099-B55G) 7.53 27.27 45.26 0.0940 0.3407 59.22 -109.51 -19.56 (0099-B50G) 7.95 27.00 48.06 0.0958 0.3252 58.97 -104.51 -23.02 (0099-B45G) 8.36 26.57 50.74 0.0975 0.3102 58.57 -99.11 -26.48 (0099-B40G) 8.75 26.00 53.31 0.0994 0.2952 58.04 -93.34 -29.99 (0099-B35G) 9.08 25.35 56.11 0.1003 0.2800 57.41 -87.85 -33.77 (0099-B30G) 9.43 24.60 59.00 0.1014 0.2644 56.68 -81.80 -37.74 (0099-B25G) 9.82 23.75 62.04 0.1027 0.2484 55.84 -75.05 -41.96 (0099-B20G) 10.33 22.80 65.04 0.1052 0.2323 54.87 -66.88 -46.26 (0099-B15G) 10.86 21.78 68.63 0.1072 0.2151 53.79 -58.23 -51.15 (0099-B10G) 11.52 20.70 72.73 0.1098 0.1972 52.62 -48.33 -56.52 (0099-B05G) 12.33 19.54 77.38 0.1129 0.1789 51.31 -37.03 -62.42 (0099-B00G) 13.30 18.30 82.68 0.1164 0.1601 49.86 -24.30 -68.92 (0099-R95B) 14.42 16.98 89.40 0.1194 0.1406 48.23 -10.19 -76.53 (0099-R90B) 15.91 15.70 96.79 0.1239 0.1223 46.58 5.84 -84.42 (0099-R85B) 17.79 14.58 102.91 0.1315 0.1078 45.05 22.85 -91.01 (0099-R80B) 19.96 13.60 106.75 0.1422 0.0969 43.65 40.05 -95.84 (0099-R75B) 22.17 12.71 106.96 0.1563 0.0896 42.32 56.38 -98.26 (0099-R70B) 24.54 11.90 106.15 0.1721 0.0835 41.06 72.44 -99.94 (0099-R65B) 27.40 11.17 105.72 0.1899 0.0774 39.87 89.50 -101.73 (0099-R60B) 30.76 10.50 102.12 0.2145 0.0732 38.73 107.40 -101.42 (0099-R55B) 33.65 9.87 91.81 0.2486 0.0729 37.61 122.64 -96.52 (0099-R50B) 33.49 9.30 73.19 0.2888 0.0802 36.56 126.63 -84.59 (0099-R45B) 30.17 8.83 52.30 0.3305 0.0967 35.65 118.44 -67.57 (0099-R40B) 26.93 8.40 36.38 0.3755 0.1171 34.80 109.42 -51.20 (0099-R35B) 24.21 8.00 24.48 0.4271 0.1411 33.98 101.53 -35.43 (0099-R30B) 22.28 7.80 16.38 0.4796 0.1679 33.56 94.68 -20.92 (0099-R25B) 21.64 7.95 11.54 0.5262 0.1933 33.88 90.34 -8.64 (0099-R20B) 21.86 8.30 8.35 0.5676 0.2156 34.60 88.22 2.27 (0099-R15B) 22.34 8.70 6.17 0.6003 0.2337 35.40 87.04 11.78 (0099-R10B) 23.12 9.30 4.22 0.6310 0.2538 36.56 85.59 22.91 (0099-R05B) 24.79 10.27 2.67 0.6571 0.2722 38.32 85.33 35.54 %(0099-R00B) 27.17 11.60 1.39 0.6765 0.2889 40.57 85.51 50.80 (9900-W99N) 0.74 0.78 0.85 0.3127 0.3290 7.05 0.00 0.00 (9500-W95N) 2.40 2.53 2.75 0.3127 0.3290 18.05 0.00 0.00 (9000-W90N) 4.18 4.40 4.79 0.3127 0.3290 24.95 0.00 0.00 (8500-W85N) 6.09 6.41 6.98 0.3127 0.3290 30.42 0.00 0.00 (8000-W80N) 8.14 8.57 9.33 0.3127 0.3290 35.14 0.00 0.00 (7500-W75N) 10.36 10.90 11.86 0.3127 0.3290 39.41 0.00 0.00 (7000-W70N) 12.75 13.42 14.61 0.3127 0.3290 43.38 0.00 0.00 (6500-W65N) 15.35 16.15 17.58 0.3127 0.3290 47.17 0.00 0.00 (6000-W60N) 18.18 19.13 20.83 0.3127 0.3290 50.84 0.00 0.00 (5500-W55N) 21.28 22.39 24.37 0.3127 0.3290 54.43 0.00 0.00 (5000-W50N) 24.68 25.96 28.27 0.3127 0.3290 58.00 0.00 0.00 (4500-W45N) 28.43 29.91 32.57 0.3127 0.3290 61.58 0.00 0.00 (4000-W40N) 32.59 34.29 37.33 0.3127 0.3290 65.19 0.00 0.00 (3500-W35N) 37.23 39.17 42.65 0.3127 0.3290 68.87 0.00 0.00 (3000-W30N) 42.44 44.65 48.61 0.3127 0.3290 72.66 0.00 0.00 (2500-W25N) 48.32 50.84 55.36 0.3127 0.3290 76.58 0.00 0.00 (2000-W20N) 55.04 57.91 63.05 0.3127 0.3290 80.69 0.00 0.00 (1500-W15N) 62.76 66.03 71.89 0.3127 0.3290 85.01 0.00 0.00 (1000-W10N) 71.74 75.48 82.18 0.3127 0.3290 89.62 0.00 0.00 (0500-W05N) 82.31 86.60 94.29 0.3127 0.3290 94.57 0.00 0.00 (0000-W00N) 94.94 99.89 108.76 0.3127 0.3290 99.96 0.00 0.00 ] def %END NCS_0099_ALLCD /DecodeXYZ* {dup 6 29 div ge {dup dup mul mul} {4 29 div sub 108 841 div mul} ifelse} bind def /cvishow {cvi 6 string cvs show} def /cvsshow1 {10 mul cvi 0.1 mul 7 string cvs show} def /cvsshow2 {100 mul cvi 0.01 mul 7 string cvs show} def /cvsshow3 {1000 mul cvi 0.001 mul 7 string cvs show} def /cvishowb {cvi 6 string cvs show ( ) show} def /cvsshow1b {10 mul cvi 0.1 mul 7 string cvs show ( ) show} def /cvsshow2b {100 mul cvi 0.01 mul 7 string cvs show ( ) show} def /cvsshow3b {1000 mul cvi 0.001 mul 7 string cvs show ( ) show} def /cvishowr {0.5 add cvi 6 string cvs show} def /cvsshow1r {0.05 add 10 mul cvi 0.1 mul 7 string cvs show} def /cvsshow2r {0.005 add 100 mul cvi 0.01 mul 7 string cvs show} def /cvsshow3r {0.0005 add 1000 mul cvi 0.001 mul 7 string cvs show} def /cvi100 {100 mul cvi /xi exch def xi 0 lt {/vorz -1 def} {/vorz 1 def} ifelse /xia xi abs def vorz -1 eq {(\255) show} {(0) show} ifelse xia 0 ge xia 9 le {(0000) show} if xia 10 ge xia 99 le and {(000) show} if xia 100 ge xia 999 le and {(00) show} if xia 1000 ge xia 9999 le and {(0) show} if xia cvishow } bind def /cvi1000 {1000 mul cvi /xi exch def xi 0 ge xi 9 le {(000) show} if xi 10 ge xi 99 le and {(00) show} if xi 100 ge xi 999 le and {(0) show} if xi cvishow } bind def /rec %x, y width heigth {/heigth exch def /width exch def moveto width 0 rlineto 0 heigth rlineto width neg 0 rlineto closepath } bind def /colrecfi %x y width heigth c m y k {setcmykcolor rec fill} bind def /colrecst %x y width heigth c m y k {setcmykcolor rec stroke} bind def %%EndProlog gsave %lanindL2.EPS START 20000505 /lanind 0 def /lantex [(G) (E) (S) (F) (I) (J) (M)] def /showde {0 lanind eq {show} {pop} ifelse} bind def /showen {1 lanind eq {show} {pop} ifelse} bind def /showes {2 lanind eq {show} {pop} ifelse} bind def /showfr {3 lanind eq {show} {pop} ifelse} bind def /showit {4 lanind eq {show} {pop} ifelse} bind def /showjp {5 lanind eq {show} {pop} ifelse} bind def /showea {1 lanind le {show} {pop} ifelse} bind def /lanindf where {pop /lanind1 lanindf def /lanind2 lanindf def /lanindd laninddf def} {/lanind1 0 def /lanind2 0 def} ifelse /colormf where {pop /colorm1 colormf def /colorm2 colormf def /colormd colormdf def} {/colorm1 0 def} ifelse /deintpf where {pop /deintp1 deintpf def /deintP2 deintpf def /deintpd deintpdf def} {/deintp1 0 def} ifelse /xcolorf where {pop /xcolor1 xcolorf def /xcolor2 xcolorf def /xcolord xcolordf def} {/xcolor1 0 def} ifelse /xchartf where {pop /xchart1 xchartf def /xchart2 xchartf def /xchartd xchartdf def /xchartm xchart2f xchart1f sub 1 add def} {/xchart1 0 def /xchart2 0 def /xchartm 1 def} ifelse /xchart3f where {pop /xchart3 xchart3f def} {/xchart3 0 def} ifelse /xchart4f where {pop /xchart4 xchart4f def} {/xchart4 0 def} ifelse /pchartf where {pop /pchart1 pchartf def /pchart2 pchartf def /pchartd pchartdf def} {/pchart1 3 def} ifelse /colsepf where {pop /colsep1 colsepf def /colseP2 colsepf def /colsepd colsepdf def} {/colsep1 0 def} ifelse /pmetamf where {pop /pmetam1 pmetamf def /pmetam2 pmetamf def /pmetamd pmetamdf def} {/pmetam1 0 def} ifelse %either defaul values for xchart=0 or values for xchart=1 /lanind lanind1 def % /colorm colorm1 def % /deintp deintp1 def % /xcolor xcolor1 def % /xchart xchart1 def % /pchart pchart1 def % /colsep colsep1 def % /pmetam pmetam1 def % colorm 0 eq deintp 0 eq and {/Txx (d) def} if %colorm=0, deintp=0 colorm 0 eq deintp 1 eq and {/Txx (e) def} if %colorm=0, deintp=1 colorm 1 eq deintp 0 eq and {/Txx (dd) def} if %colorm=1, deintp=0 colorm 1 eq deintp 1 eq and {/Txx (de) def} if %colorm=1, deintp=1 xchart 0 eq {/Txx (-) def} if %always independent of intended output gsave %LANINDL2 END /SD65 80 array def %380 bis 775 mit deltalambda=05 /SD50 80 array def %380 bis 775 mit deltalambda=05 /SP40 80 array def %380 bis 775 mit deltalambda=05 /SA00 80 array def %380 bis 775 mit deltalambda=05 /SE00 80 array def %380 bis 775 mit deltalambda=05 /SC00 80 array def %380 bis 775 mit deltalambda=05 /SP00 80 array def %380 bis 775 mit deltalambda=05 /SQ00 80 array def %380 bis 775 mit deltalambda=05 /SDAK 80 array def %380 bis 775 mit deltalambda=05 /RA 80 array def %380 bis 770 nm /F2 240 array def %3 x 80 Daten fuer Bereich 380 bis 770nm /F2N 240 array def %normiert auf Y=100 fuer D65 und 10nm %C %C LITERATURQUELLEN LICHTARTEN %C SA,SC,SD65, CIE-DOKUMENT NR.15 FARBMESSUNG, 1971 %C TABELLE 1.1.1, SEITE 81 %C SD55,SD75, CIE-DOKUMENT NR.15 FARBMESSUNG, 1971 %C TABELLE 1.1.4, SEITE 85 %C LINEAR INTERPOLIERT VON 5 ZU 5NM %C SG,SP,SXE, DIN 5033, BLATT 7 ,SEITE 8-9 %C AUSGABE 1966 /SD65 %80 data 5 nm [50.0 52.3 54.6 68.7 82.8 87.1 91.5 92.5 93.4 90.1 86.7 95.8 104.9 110.9 117.0 117.4 117.8 116.3 114.9 115.4 115.9 112.4 108.8 109.1 109.4 108.6 107.8 106.3 104.8 106.2 107.7 106.0 104.4 104.2 104.0 102.0 100.0 98.2 96.3 96.1 95.8 92.2 88.7 89.3 90.0 89.8 89.6 88.6 87.7 85.5 83.3 83.5 83.7 81.9 80.0 80.1 80.2 81.2 82.3 80.3 78.3 74.0 69.7 70.7 71.6 73.0 74.3 68.0 61.6 65.7 69.9 72.5 75.1 69.3 63.6 55.0 46.4 56.6 66.8 65.1] def /SD50 %80 data 5 nm CIE 15.3, Seite 30; aus JE540-1N.PS [24.5 27.2 29.9 39.6 49.3 52.9 56.5 58.3 60.0 58.9 57.8 66.3 74.8 81.0 87.2 88.9 90.6 91.0 91.4 93.2 95.1 93.5 92.0 93.8 95.7 96.2 96.6 96.8 97.1 99.6 102.1 101.4 100.8 101.5 102.3 101.2 100.0 98.9 97.7 98.3 98.9 96.2 93.5 95.6 97.7 98.5 99.3 99.2 99.0 97.4 95.7 97.3 98.9 97.3 95.7 96.9 98.2 100.6 103.0 101.1 99.1 93.3 87.4 89.5 91.6 92.2 92.9 84.9 76.9 81.7 86.5 89.5 92.6 85.4 78.2 68.0 57.7 70.3 82.9 0.0] def %0.0 not used /SA00 %80 data 5 nm CIE 15.3 Seite 30; aus JE540-1N.PS [9.80 10.90 12.09 13.35 14.71 16.15 17.68 19.29 20.99 22.79 24.67 26.64 28.70 30.85 33.09 35.41 37.81 40.30 42.87 45.52 48.24 51.04 53.91 56.85 59.86 62.93 66.06 69.25 72.50 75.79 79.13 82.52 85.95 89.41 92.91 96.44 100.00 103.58 107.18 110.80 114.44 118.08 121.73 125.39 129.04 132.70 136.35 139.99 143.62 147.24 150.84 154.42 157.98 161.52 165.03 168.51 171.96 175.38 178.77 182.12 185.43 188.70 191.93 195.12 198.26 201.36 204.41 207.41 210.36 213.27 216.12 218.92 221.67 224.36 227.00 229.59 232.12 234.59 237.01 239.37 0.0] def %0.0 not used % AUS BILDPLA.FOR aus vrbam/wrpos:KRPROG.FOR, 5987-6023 5.10.92 %C PROGRAMM ZUR BERECHNUNG DER RELATIVEN SPEKTRALEN STRAHLUNGSVERTEI- %C LUNG EINES HOHLRAUMSTRAHLERS ALS FUNKTION DER WELLENLAENGE LAMDA %C BEI GEGEBENER TEMPERATUR T %C VERWENDETE VARIABLE %C T TEMPERATUR %C C2 PLANCKSCHE KONSTANTE (=1.4388 CM K (1967)) %C LAMDA(I) WELLENLAENGE %C S(I) STRAHLUNGSVERTEILUNG DES HOHLRAUMSTRAHLERS %C %C LITERATUR OFFIZIELLE EMPFEHLUNGEN DER CIE,DEUTSCHE VERSION, %C CA. 1969, SEITE 21 UND TABELLE 1.1.1 %C %SDP4 80 array def %380 bis 775 mit deltalambda=05 % DO 50 I=1,54 % 50 S(I)=100. % 1*(560./FLOAT(LAMDA(I)))**5 % 2*(EXP(C2/(560.*10.**(-7)*T))-1.) % 3/(EXP(C2/(FLOAT(LAMDA(I))*10.**(-7)*T))-1.) %new CIE 15.3:2004, Eq. (3.1) page 2 %S(I) = [100{560/lambda}**5] [exp{(c2*10**7)/(TNC*560 )} - 1] % /[exp{(c2*10**7)/(TNC*lambda)} - 1] % exp = e {...} = 2.71 {...} %/TNC 2848 def %NLA /TNC 4000 def /EE 2.718281 def %e=mac pocket calculator /C2 1.4350000 10 7 exp mul def /lamd 80 array def 0 1 79 {%i=0,79 /i exch def /lambda 380 i 5 mul add def lamd i lambda put /E1 560 lambda div 5 exp def /E4 C2 TNC 560 mul div def /E5 C2 TNC lambda mul div def /E6 EE E4 exp def /E7 EE E5 exp def SP40 i E1 100 mul E6 1 sub E7 1 sub div mul put } for %i=0,79 /xex 10 2 exp def /yex EE 2 exp def /SE00 %80 data 5 nm [100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0] def %0.0 not used /SC00 %80 data 5 nm [ 33.00 39.92 47.40 55.17 63.30 71.81 80.60 89.53 98.10 105.80 112.40 117.75 121.50 123.45 124.00 123.60 123.10 123.30 123.80 124.09 123.90 122.92 120.70 116.90 112.10 106.98 102.30 98.81 96.90 96.78 98.00 99.94 102.10 103.95 105.20 105.67 105.30 104.11 102.30 100.15 97.80 95.43 93.20 91.22 89.70 88.83 88.40 88.19 88.10 88.06 88.00 87.86 87.80 87.99 88.20 88.20 87.90 87.22 86.30 85.30 84.00 82.21 80.20 78.24 76.30 74.36 72.40 70.40 68.30 66.30 64.40 62.80 61.50 60.20 59.20 58.50 58.10 58.00 58.20 58.50 ] def %0.0 not used 0 1 79 {/i exch def %i=0,79 /lamda 380 i 5 mul add def SP00 i SE00 i get lamda 560 div mul put SQ00 i SE00 i get 1120 lamda sub 560 div mul put } for %i=0,79 %C LITERATURQUELLEN SPEKTRALWERTE %C F2 CIE-DOKUMENT NR.15 FARBMESSUNG, 1971 %C TABELLE 2.1, SEITE 93-102 /F2 %3x80 data [.001368 .002236 .004243 .007650 .01431 .02319 .04351 .07763 .1344 .2148 .2839 .3285 .3483 .3481 .3362 .3187 .2908 .2511 .1954 .1421 .09564 .05795 .03201 .01470 .0049 .0024 .0093 .0291 .06327 .1096 .1655 .2257 .2904 .3597 .4334 .5121 .5945 .6784 .7621 .8425 .9163 .9786 1.0263 1.0567 1.0622 1.0456 1.0026 .9384 .8544 .7514 .6424 .5419 .4479 .3608 .2835 .2187 .1649 .1212 .0874 .0636 .04677 .03290 .02270 .01584 .01136 .00811 .00579 .004109 .002899 .002049 .001440 .0009999 .0006901 .000476 .0003323 .0002348 .0001662 .0001174 .00008308 .00005871 .000039 .000064 .00012 .000217 .000396 .00064 .00121 .00218 .004 .0073 .0116 .01684 .023 .0298 .0380 .048 .06 .0739 .09098 .1126 .139 .1693 .208 .2586 .323 .4073 .503 .6082 .71 .7932 .862 .9149 .954 .9803 .995 1.0 .995 .9786 .952 .9154 .870 .8163 .757 .6949 .631 .5668 .503 .4412 .381 .321 .265 .217 .175 .1382 .107 .0816 .061 .04458 .032 .0232 .017 .01192 .00821 .005723 .004102 .002929 .002091 .001484 .001047 .00074 .00052 .0003611 .0002492 .0001719 .00012 .0000848 .00006 .0000424 .00003 .0000212 .00645 .01055 .02005 .03621 .06785 .1102 .2074 .3713 .6456 1.0391 1.3856 1.623 1.7471 1.7826 1.7721 1.7441 1.6692 1.5281 1.2876 1.0419 .8130 .6162 .4652 .3533 .2720 .2123 .1582 .1117 .07825 .05725 .04216 .02984 .0203 .0134 .00875 .00575 .0039 .00275 .0021 .0018 .00165 .0014 .0011 .0010 .0008 .0006 .00034 .00024 .00019 .0001 .00005 .00003 .00002 .00001 % 26*0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 %10 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 %20 0.0 0.0 0.0 0.0 0.0 0.0 %26 ] def /F10 %3x80 data [.00016 .0006624 .002362 .007242 .01911 .0434 .08474 .1406 .2045 .2647 .3147 .3577 .3837 .3867 .3707 .343 .3023 .2541 .1956 .1323 .08051 .04107 .01617 .005132 .003816 .01544 .03746 .07136 .1177 .173 .2365 .3042 .3768 .4516 .5298 .6161 .7052 .7938 .8787 .9512 1.0142 1.0743 1.1185 1.1343 1.124 1.0891 1.0305 .9507 .8563 .7549 .6475 .5351 .4316 .3437 .2683 .2043 .1526 .1122 .08126 .05793 .04085 .02862 .019941 .01384 .0100 .006605 .004552 .003145 .002175 .001506 .001045 .0007274 .0005083 .0003564 .000251 .0001777 .0001264 .00009015 .00006453 .00004634 .00001736 .00007156 .0002534 .0007685 .002004 .004509 .008756 .01446 .02139 .0295 .03868 .0496 .06208 .0747 .08946 .1063 .1282 .1528 .1852 .2199 .2536 .2977 .3391 .3954 .4608 .5314 .6067 .6857 .7618 .8233 .8752 .9238 .962 .9822 .9918 .9991 .9973 .9824 .9555 .9152 .8689 .8256 .7774 .7203 .6583 .5939 .528 .4618 .3981 .3395 .2835 .2282 .1798 .1402 .1076 .08119 .06028 .0441 .0318 .0226 .0159 .01113 .007749 .0053751 .0037177 .002565 .001768 .001222 .0008462 .0005864 .0004074 .000284 .0001987 .0001395 .00009843 .00006982 .00004974 .00003554 .00002549 .00001834 .0007048 .002928 .01048 .03234 .08601 .1971 .3894 .6568 .9725 1.2825 1.5535 1.7985 1.9673 2.0273 1.9948 1.9007 1.7454 1.5549 1.3176 1.0302 .7721 .5701 .4152 .3024 .2185 .1592 .112 .08225 .06071 .04305 .03045 .02058 .01368 .007918 .003988 .001091 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 %10 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 %20 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 %30 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 %40 0.0 0.0 0.0 0.0 %44 ] def 0 1 100 {/i exch def %i=0,1,101 /i30 i 9 mul 1 add def /i31 i30 1 add def /i32 i30 2 add def /j30 i 3 mul def /j31 j30 1 add def /j32 j30 2 add def XYZC00100 j30 NCS_0099_ALLCD i30 get put XYZC00100 j31 NCS_0099_ALLCD i31 get put XYZC00100 j32 NCS_0099_ALLCD i32 get put /XQ XYZC00100 j30 get XLWC00100 div def /YQ XYZC00100 j31 get YLWC00100 div def /ZQ XYZC00100 j32 get ZLWC00100 div def XQ 0 lt {/XQ 0.00000001 def} if YQ 0 lt {/YQ 0.00000001 def} if ZQ 0 lt {/ZQ 0.00000001 def} if LABC00100 j30 YQ 0.008856 lt {903.3 YQ mul} {YQ 0.33333333 exp 116 mul 16 sub} ifelse put LABC00100 j31 XQ 0.008856 lt {7.787 XQ mul 16 116 div add} {XQ 0.33333333 exp} ifelse YQ 0.008856 lt {7.787 YQ mul 16 116 div add} {YQ 0.33333333 exp} ifelse sub 500 mul put LABC00100 j32 YQ 0.008856 lt {7.787 YQ mul 16 116 div add} {YQ 0.33333333 exp} ifelse ZQ 0.008856 lt {7.787 ZQ mul 16 116 div add} {ZQ 0.33333333 exp} ifelse sub 200 mul put XYZC00089 j30 NCS_0099_ALLCD i30 get F89 mul put XYZC00089 j31 NCS_0099_ALLCD i31 get F89 mul put XYZC00089 j32 NCS_0099_ALLCD i32 get F89 mul put /XQ XYZC00089 j30 get XLWC00100 div def /YQ XYZC00089 j31 get YLWC00100 div def /ZQ XYZC00089 j32 get ZLWC00100 div def XQ 0 lt {/XQ 0.00000001 def} if YQ 0 lt {/YQ 0.00000001 def} if ZQ 0 lt {/ZQ 0.00000001 def} if LABC00089 j30 YQ 0.008856 lt {903.3 YQ mul} {YQ 0.33333333 exp 116 mul 16 sub} ifelse put LABC00089 j31 XQ 0.008856 lt {7.787 XQ mul 16 116 div add} {XQ 0.33333333 exp} ifelse YQ 0.008856 lt {7.787 YQ mul 16 116 div add} {YQ 0.33333333 exp} ifelse sub 500 mul put LABC00089 j32 YQ 0.008856 lt {7.787 YQ mul 16 116 div add} {YQ 0.33333333 exp} ifelse ZQ 0.008856 lt {7.787 ZQ mul 16 116 div add} {ZQ 0.33333333 exp} ifelse sub 200 mul put } for %i=0,1,101 72 90 translate 0.01 MM dup scale 20 setlinewidth 1 1 1 setrgbcolor 0 0 moveto 25000 0 rlineto 0 17000 rlineto 25000 neg 0 rlineto closepath fill 0 0 0 setrgbcolor 0 0 moveto 25000 0 rlineto 0 17000 rlineto 25000 neg 0 rlineto closepath stroke TK 100 -180 moveto (TG030-7) show xchart 0 eq {(N) show} if (, NCS colours) showen (, NCS\255Farben) showde TK (, normalized:) showen (, normiert:) showde TIK ( Yn = Yw) show TK ( = 100) show /ymax 16600 def /idy 206 def /xcolorl xcolor def %0 1 1 {/xcolorl exch def %xcolorl=0,1 %xcolorl 1 eq {0 1 79 {/i exch def %i=0,79 % F2 i F10 i get put % } for %i=0,79 % } if /YKSUM 0 def %For CIE Illuminant E 0 1 79 {/i exch def %i=0,1,79 /YKSUM YKSUM F2 80 i add get add def } for %i=0,1,79 /F2N 240 array def 0 1 79 {/i exch def %i=0,1,79 F2N 00 i add F2 00 i add get YKSUM div 100 mul put F2N 80 i add F2 80 i add get YKSUM div 100 mul put F2N 160 i add F2 160 i add get YKSUM div 100 mul put } for %i=0,1,79 /FXN 80 array def %FXN 0 = F2N imax /FYN 80 array def /FZN 80 array def %defines FXN, FYN, FZN 0 to 79 0 1 79 {/i exch def %i=0,79 FXN i F2N 000 i add get put FYN i F2N 080 i add get put FZN i F2N 160 i add get put } for %i=0,79 %special wavelength range with complementary colours for E? /FXD 63 array def /FYD 63 array def /FZD 63 array def /DOW 63 array def 00 1 60 {/i exch def %i=0,60 %range 400..700nm FXD i FXN 005 i add get put FYD i FYN 005 i add get put FZD i FZN 005 i add get put DOW i 400 i 5 mul add put %range 400 to 700 } for %i=0,60 %for interpolation FXD 61 FXD 60 get put FYD 61 FYD 60 get put FZD 61 FZD 60 get put FXD 62 FXD 60 get put FYD 62 FYD 60 get put FZD 62 FZD 60 get put /FAD 63 array def /FBD 63 array def /FCD 63 array def /proc_ABC_KNP_FWERTiX {%BEG proc_ABC_KNP_FWERT %input FAD,FBD,FCD, KNP %output FWERTiX/Y/Z (i=1,4) kNP 1 ge {/FWERT1X FAD kNP 1 sub get def /FWERT1Y FBD kNP 1 sub get def /FWERT1Z FCD kNP 1 sub get def} {/FWERT1X FAD 0 get def %again 400nm for i=0 /FWERT1Y FBD 0 get def /FWERT1Z FCD 0 get def} ifelse /FWERT2X FAD kNP get def /FWERT2Y FBD kNP get def /FWERT2Z FCD kNP get def /FWERT3X FAD kNP 1 add get def /FWERT3Y FBD kNP 1 add get def /FWERT3Z FCD kNP 1 add get def /FWERT4X FAD kNP 2 add get def /FWERT4Y FBD kNP 2 add get def /FWERT4Z FCD kNP 2 add get def } bind def %END proc_ABC_KNP_FWERTiX %********************************************************* /proc_XIE_TN1 {%BEG proc_XIE_TN1 %input XIE, FWERTiX/Y/Z (i=1,4) %output TN1 /A1N XIE 1 XIE sub mul XIE 2 sub mul 6 div def /A2N XIE 1 add XIE 1 sub mul XIE 2 sub mul 0.5 mul def /A3N XIE 1 add XIE mul 2 XIE sub mul 0.5 mul def /A4N XIE 1 add XIE mul XIE 1 sub mul 6 div def /FNX A1N FWERT1X mul A2N FWERT2X mul add A3N FWERT3X mul add A4N FWERT4X mul add def /FNY A1N FWERT1Y mul A2N FWERT2Y mul add A3N FWERT3Y mul add A4N FWERT4Y mul add def /FNZ A1N FWERT1Z mul A2N FWERT2Z mul add A3N FWERT3Z mul add A4N FWERT4Z mul add def /TNX FNX FXU1 mul FNY FXU2 mul add FNZ FXU3 mul add def } bind def %END proc_XIE_TN1 xchart1 1 xchart2 {/xchart exch def %BEG xchart=0,1 %for 80 and 21 steps xchart 0 eq {/IM0 00 def /IM1 079 def} if xchart 1 eq {/IM0 80 def /IM1 100 def} if /xchartl 5 def %for CIE illuminant C gsave %Start-transformation 1 dup dup setrgbcolor 0 0 moveto 25000 0 rlineto 0 17000 rlineto 25000 neg 0 rlineto closepath fill 0 dup dup setrgbcolor 0 0 moveto 25000 0 rlineto 0 17000 rlineto 25000 neg 0 rlineto closepath stroke %calculation of XW,YW,ZW for illuminant D65, D50, ... 0 1 79 {/i exch def %i=0,79 SDAK i xchartl 0 eq {SD65 i get} if xchartl 1 eq {SD50 i get} if xchartl 2 eq {SP40 i get} if xchartl 3 eq {SA00 i get} if xchartl 4 eq {SE00 i get} if xchartl 5 eq {SC00 i get} if xchartl 6 eq {SP00 i get} if xchartl 7 eq {SQ00 i get} if put } for %i=0,79 /YKSUM 0 def 0 1 79 {/i exch def %i=0,1,79 /YKSUM YKSUM SDAK i get F2 80 i add get mul add def } for %i=0,1,79 0 1 79 {/i exch def %i=0,1,79 F2N 00 i add SDAK i get F2 00 i add get mul YKSUM div 100 mul put F2N 80 i add SDAK i get F2 80 i add get mul YKSUM div 100 mul put F2N 160 i add SDAK i get F2 160 i add get mul YKSUM div 100 mul put } for %i=0,1,79 %XW,YW,ZW for white /XW 0 def /YW 0 def /ZW 0 def 0 1 079 {/i exch def %i=0,79 /XW XW F2N 00 i add get add def /YW YW F2N 80 i add get add def /ZW ZW F2N 160 i add get add def } for %i=0,79 /SUMW XW YW add ZW add def %for D65, D50, P40, A00, E00 /XANW XW SUMW div def /YANW YW SUMW div def /ZANW ZW SUMW div def /AANW XW YW 0.0001 add div def /BANW ZW YW 0.0001 add div 0.4 mul neg def 0 1 62 {/i exch def %i=0,62 FAD i FXD i get put FBD i FYD i get put FCD i FZD i get put } for %i=0,62 /ix0 100 def /ix1 ix0 1300 add def /ixd 600 def /iymax 16550 def /idy 206 def TBIK ix0 iymax idy add 50 add moveto (Cod.) show ix1 ixd 0 mul add iymax idy add 50 add moveto TBK (i) show ix1 ixd 1 mul add 030 sub 0 moveto 0 17000 rlineto stroke ix1 ixd 1 mul add iymax idy add 50 add moveto TBK ([) show TBIK (X, Y, Z, x, y) show TBK (]) show TBS 0 -90 rmoveto (100) show 0 90 rmoveto TBK ix1 ixd 6 mul add 030 sub 0 moveto 0 17000 rlineto stroke ix1 ixd 6 mul add iymax idy add 50 add moveto TBK ([) show TBIK (L*, a*, b*, C*) show TBS -80 -90 rmoveto (ab) show 0 90 rmoveto TBIK (, h) show TBS 0 -90 rmoveto (ab) show 0 90 rmoveto TBIK TBIK (, a', b', c') show TBS 0 -90 rmoveto (ab) show 0 90 rmoveto TBIK TBK (]) show TBS 0 -90 rmoveto (100) show 0 90 rmoveto TBK ix1 ixd 14 mul add 030 sub 0 moveto 0 17000 rlineto stroke ix1 ixd 14 mul add iymax idy add 50 add moveto TBK ([) show TBIK (Y, A, B, C) show TBS -80 -90 rmoveto (AB) show 0 90 rmoveto TBIK (, h) show TBS 0 -90 rmoveto (AB) show 0 90 rmoveto TBIK TBIK (, a, b, c) show TBS 0 -90 rmoveto (AB) show 0 90 rmoveto TBIK TBK (]) show TBS 0 -90 rmoveto (100) show 0 90 rmoveto TBK ix1 ixd 22 mul add 030 sub 0 moveto 0 17000 rlineto stroke ix1 ixd 22 mul add iymax idy add 50 add moveto TBK ([) show TBIK (i) show TBS 0 -90 rmoveto (d) show 0 90 rmoveto TBK TBK (, ) show TBIK (i) show TBS 0 -90 rmoveto (c) show 0 90 rmoveto TBK TBK (, ) show TSK (l) show TBS 0 -90 rmoveto (d) show 0 90 rmoveto TBK TBK (, ) show TSK (l) show TBS 0 -90 rmoveto (c) show 0 90 rmoveto TBK TBK (]) show TBS 0 -90 rmoveto (100) show 0 90 rmoveto TBK ix1 ixd 26 mul add 030 sub 0 moveto 0 17000 rlineto stroke ix1 ixd 26 mul add iymax idy add 50 add moveto TBK ([) show TBIK (X, Y, Z) show TBK (]) show TBS 0 -90 rmoveto (89) show 0 90 rmoveto TBK ix1 ixd 29 mul add 030 sub 0 moveto 0 17000 rlineto stroke ix1 ixd 29 mul add iymax idy add 50 add moveto TBK ([) show TBIK (L*, a*, b*, C*) show TBS -80 -90 rmoveto (ab) show 0 90 rmoveto TBIK (, h) show TBS 0 -90 rmoveto (ab) show 0 90 rmoveto TBIK TBK (]) show TBS 0 -90 rmoveto (89) show 0 90 rmoveto TBK ix1 ixd 34 mul add 030 sub 0 moveto 0 17000 rlineto stroke ix1 ixd 34 mul add iymax idy add 50 add moveto TBK ([) show TBIK (Y, A, B, C) show TBS -80 -90 rmoveto (AB) show 0 90 rmoveto TBIK (, h) show TBS 0 -90 rmoveto (AB) show 0 90 rmoveto TBIK TBK (]) show TBS 0 -90 rmoveto (89) show 0 90 rmoveto TBK /D13 1 3 div def /a2C00 1 XLWC00100 0.0001 add div D13 exp def /b2C00 1 ZLWC00100 0.0001 add div D13 exp neg def TK %20000 15000 moveto a2C00 cvsshow3r (,) b2C00 cvsshow3r /aWC00100 XLWC00100 YLWC00100 0.0001 add div def /bWC00100 ZLWC00100 YLWC00100 0.0001 add div 0.4 mul neg def /aWC00089 XLWC00100 YLWC00100 0.0001 add div def /bWC00089 ZLWC00100 YLWC00100 0.0001 add div 0.4 mul neg def %20000 14000 moveto aWC00100 cvsshow3r (,) bWC00100 cvsshow3r /a'WC00100 XLWC00100 YLWC00100 0.0001 add div D13 exp a2C00 mul def /b'WC00100 ZLWC00100 YLWC00100 0.0001 add div D13 exp b2C00 mul def /a'WC00089 XLWC00100 YLWC00100 0.0001 add div D13 exp a2C00 mul def /b'WC00089 ZLWC00100 YLWC00100 0.0001 add div D13 exp b2C00 mul def %20000 13000 moveto a'WC00100 cvsshow3r (,) b'WC00100 cvsshow3r xchart 1 eq {/iymax iymax 80 idy mul add def} if IM0 1 IM1 {/i exch def %i=0,IM1 /i30 i 9 mul def ix0 ixd 0 mul add iymax i idy mul sub moveto NCS_0099_ALLCD i30 get show ix1 ixd 0 mul add iymax i idy mul sub moveto i 9 le {(0) show} if i cvishow /j30 i 3 mul def /j31 j30 1 add def /j32 j30 2 add def /XC00100 XYZC00100 j30 get def /YC00100 XYZC00100 j31 get def /ZC00100 XYZC00100 j32 get def /SUM XC00100 YC00100 add ZC00100 add def /xC00100 XC00100 SUM 0.0001 add div def /yC00100 YC00100 SUM 0.0001 add div def ix1 ixd 1 mul add iymax i idy mul sub moveto XC00100 cvsshow1r ix1 ixd 2 mul add iymax i idy mul sub moveto YC00100 cvsshow1r ix1 ixd 3 mul add iymax i idy mul sub moveto ZC00100 cvsshow1r ix1 ixd 4 mul add iymax i idy mul sub moveto xC00100 cvsshow3r ix1 ixd 5 mul add iymax i idy mul sub moveto yC00100 cvsshow3r /a'C00100 XC00100 YC00100 0.0001 add div D13 exp a2C00 mul def /b'C00100 ZC00100 YC00100 0.0001 add div D13 exp b2C00 mul def /c'C00100 a'C00100 a'WC00100 sub dup mul b'C00100 b'WC00100 sub dup mul add 0.0001 add sqrt def /L*C00100 LABC00100 j30 get def /a*C00100 LABC00100 j31 get def /b*C00100 LABC00100 j32 get def /C*C00100 a*C00100 dup mul b*C00100 dup mul add 0.0001 add sqrt def /h*C00100 b*C00100 a*C00100 0.0001 add atan def ix1 ixd 06 mul add iymax i idy mul sub moveto L*C00100 cvsshow1r ix1 ixd 07 mul add iymax i idy mul sub moveto a*C00100 cvsshow1r ix1 ixd 08 mul add iymax i idy mul sub moveto b*C00100 cvsshow1r ix1 ixd 09 mul add iymax i idy mul sub moveto C*C00100 cvsshow1r ix1 ixd 10 mul add iymax i idy mul sub moveto h*C00100 cvsshow1r ix1 ixd 11 mul add iymax i idy mul sub moveto a'C00100 cvsshow3r ix1 ixd 12 mul add iymax i idy mul sub moveto b'C00100 cvsshow3r ix1 ixd 13 mul add iymax i idy mul sub moveto c'C00100 cvsshow3r /aC00100 XC00100 YC00100 0.0001 add div def /bC00100 ZC00100 YC00100 0.0001 add div 0.4 mul neg def /cC00100 aC00100 aWC00100 sub dup mul bC00100 bWC00100 sub dup mul add 0.0001 add sqrt def /AC00100 XC00100 YC00100 0.0001 add div XLWC00100 YLWC00100 0.0001 add div sub YC00100 mul def /BC00100 ZC00100 YC00100 0.0001 add div ZLWC00100 YLWC00100 0.0001 add div sub YC00100 mul 0.4 mul neg def /CC00100 AC00100 dup mul BC00100 dup mul add 0.0001 add sqrt def /hC00100 BC00100 AC00100 0.0001 add atan def ix1 ixd 14 mul add iymax i idy mul sub moveto YC00100 cvsshow1r ix1 ixd 15 mul add iymax i idy mul sub moveto AC00100 cvsshow1r ix1 ixd 16 mul add iymax i idy mul sub moveto BC00100 cvsshow1r ix1 ixd 17 mul add iymax i idy mul sub moveto CC00100 cvsshow1r ix1 ixd 18 mul add iymax i idy mul sub moveto hC00100 cvsshow1r ix1 ixd 19 mul add iymax i idy mul sub moveto aC00100 cvsshow3r ix1 ixd 20 mul add iymax i idy mul sub moveto bC00100 cvsshow3r ix1 ixd 21 mul add iymax i idy mul sub moveto cC00100 cvsshow3r %********************************************************* %calculation of lamdad and lambdac for X, Y, Z /FF1 XC00100 def /FF2 YC00100 def /FF3 ZC00100 def /U1 XLWC00100 def /U2 YLWC00100 def /U3 ZLWC00100 def /FXU1 FF2 U3 mul FF3 U2 mul sub def /FXU2 FF3 U1 mul FF1 U3 mul sub def /FXU3 FF1 U2 mul FF2 U1 mul sub def /INP -1 def /IPN -1 def %assumption jump from TN to TP from negativ to positiv 0 1 60 {/j exch def %j=1,60 /TN FAD j get FXU1 mul FBD j get FXU2 mul add FCD j get FXU3 mul add def /TP FAD j 1 add get FXU1 mul FBD j 1 add get FXU2 mul add FCD j 1 add get FXU3 mul add def TN 0 le TP 0 gt and {/INP j def exit} if } for %j=1,60 0 1 60 {/j exch def %j=1,60 /TN FAD j get FXU1 mul FBD j get FXU2 mul add FCD j get FXU3 mul add neg def /TP FAD j 1 add get FXU1 mul FBD j 1 add get FXU2 mul add FCD j 1 add get FXU3 mul add neg def TN 0 le TP 0 gt and {/IPN j def exit} if } for %j=1,60 TK ix1 ixd 22 mul add iymax i idy mul sub moveto INP cvishow ix1 ixd 23 mul add iymax i idy mul sub moveto IPN cvishow %********************************************************* /kNP INP def kNP 0 ge {%kNP>=0 proc_ABC_KNP_FWERTiX %input FAD,KNP, output FWERTiXYZ (i=1,4) /XIE1 0 def %start values TN=0 for XIE1, TP>0 for XIE2 /XIE2 1 def 0 1 9 {/mk exch def %mk=0,9 /XIE XIE1 XIE2 add 0.5 mul def proc_XIE_TN1 %input XIE, ouput TN1 TNX 0 le {/XIE1 XIE def /TN1 TNX def} {/XIE2 XIE def /TP1 TNX def} ifelse } for %mk=0,9 /XIE XIE1 XIE2 add 0.5 mul def /DOWXIE1 DOW INP get DOW INP 1 add get DOW INP get sub XIE mul add def } if %kNP>=0 %********************************************************* /kNP IPN def kNP 0 ge {%kNP>=0 proc_ABC_KNP_FWERTiX %input FAD,KNP, output FWERTiXYZ (i=1,4) /XIE1 0 def %start values TN=0 for XIE1, TP>0 for XIE2 /XIE2 1 def /FXU1 FXU1 neg def /FXU2 FXU2 neg def /FXU3 FXU3 neg def 0 1 9 {/mk exch def %mk=0,9 /XIE XIE1 XIE2 add 0.5 mul def proc_XIE_TN1 %input XIE, ouput TN1 TNX 0 le {/XIE1 XIE def /TN1 TNX def} {/XIE2 XIE def /TP1 TNX def} ifelse } for %mk=0,9 /XIE XIE1 XIE2 add 0.5 mul def /DOWXIE2 DOW IPN get DOW IPN 1 add get DOW IPN get sub XIE mul add def } if %kNP>=0 %********************************************************* %0 1 1 0 setcmyk_olvcolor INP 0 ge {%INP>=0,-1 ix1 ixd 24 mul add iymax i idy mul sub moveto % xpos x1 add xdel 9 mul add 400 add ymax ydel LP mul sub moveto DOWXIE1 cvishow /DOWEL1 DOWXIE1 def /DOWEL1c () def }%INP>=0 {%INP=-1 ix1 ixd 24 mul add iymax i idy mul sub moveto % xpos x1 add xdel 9 mul add 400 add ymax ydel LP mul sub moveto DOWXIE2 cvishow (c) show /DOWEL1 DOWXIE2 def /DOWEL1c (c) def } ifelse %INP>=0,-1 %1 0 1 0 setcmyk_olvcolor IPN 0 ge {%IPN>=0,-1 ix1 ixd 25 mul add iymax i idy mul sub moveto % xpos x1 add xdel 10 mul add 400 add ymax ydel LP mul sub moveto DOWXIE2 cvishow /DOWEL2 DOWXIE2 def /DOWEL2c () def }%IPN>=0 {%IPN=-1 ix1 ixd 25 mul add iymax i idy mul sub moveto % xpos x1 add xdel 10 mul add 400 add ymax ydel LP mul sub moveto DOWXIE1 cvishow (c) show /DOWEL2 DOWXIE1 def /DOWEL2c (c) def } ifelse %IPN>=0,-1 %********************************************************************** /XC00089 XYZC00089 j30 get def /YC00089 XYZC00089 j31 get def /ZC00089 XYZC00089 j32 get def ix1 ixd 26 mul add iymax i idy mul sub moveto XC00089 cvsshow1r ix1 ixd 27 mul add iymax i idy mul sub moveto YC00089 cvsshow1r ix1 ixd 28 mul add iymax i idy mul sub moveto ZC00089 cvsshow1r /a'C00089 XC00089 YC00089 0.0001 add div D13 exp a2C00 mul def /b'C00089 ZC00089 YC00089 0.0001 add div D13 exp b2C00 mul def /c'C00089 a'C00089 a'WC00089 sub dup mul b'C00089 b'WC00089 sub dup mul add 0.0001 add sqrt def /L*C00089 LABC00089 j30 get def /a*C00089 LABC00089 j31 get def /b*C00089 LABC00089 j32 get def /C*C00089 a*C00089 dup mul b*C00089 dup mul add 0.0001 add sqrt def /h*C00089 b*C00089 a*C00089 0.0001 add atan def ix1 ixd 29 mul add iymax i idy mul sub moveto L*C00089 cvsshow1r ix1 ixd 30 mul add iymax i idy mul sub moveto a*C00089 cvsshow1r ix1 ixd 31 mul add iymax i idy mul sub moveto b*C00089 cvsshow1r ix1 ixd 32 mul add iymax i idy mul sub moveto C*C00089 cvsshow1r ix1 ixd 33 mul add iymax i idy mul sub moveto h*C00089 cvsshow1r %ix1 ixd 11 mul add iymax i idy mul sub moveto a'C00089 cvsshow3r %ix1 ixd 12 mul add iymax i idy mul sub moveto b'C00089 cvsshow3r %ix1 ixd 13 mul add iymax i idy mul sub moveto c'C00089 cvsshow3r /aC00089 XC00089 YC00089 0.0001 add div def /bC00089 ZC00089 YC00089 0.0001 add div 0.4 mul neg def /cC00089 aC00089 aWC00089 sub dup mul bC00089 bWC00089 sub dup mul add 0.0001 add sqrt def /AC00089 XC00089 YC00089 0.0001 add div XLWC00089 YLWC00089 0.0001 add div sub YC00089 mul def /BC00089 ZC00089 YC00089 0.0001 add div ZLWC00089 YLWC00089 0.0001 add div sub YC00089 mul 0.4 mul neg def /CC00089 AC00089 dup mul BC00089 dup mul add 0.0001 add sqrt def /hC00089 BC00089 AC00089 0.0001 add atan def ix1 ixd 34 mul add iymax i idy mul sub moveto YC00089 cvsshow1r ix1 ixd 35 mul add iymax i idy mul sub moveto AC00089 cvsshow1r ix1 ixd 36 mul add iymax i idy mul sub moveto BC00089 cvsshow1r ix1 ixd 37 mul add iymax i idy mul sub moveto CC00089 cvsshow1r ix1 ixd 38 mul add iymax i idy mul sub moveto hC00089 cvsshow1r %ix1 ixd 19 mul add iymax i idy mul sub moveto aC00089 cvsshow3r %ix1 ixd 20 mul add iymax i idy mul sub moveto bC00089 cvsshow3r %ix1 ixd 21 mul add iymax i idy mul sub moveto cC00089 cvsshow3r } for %i=IM0,IM1 showpage grestore } for %end xchart=0,1 %} for %end xcolorl=0,1 %%Trailer