%!PS-Adobe-3.0 EPSF-3.0 SE900-7N
%%BoundingBox: 70 83 400 582

%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/SE90/)
  /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:2013090112000)
  /ModDate (D:2013090112000)
/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 findfont dup length dict  begin
{1 index /FID ne {def} {pop pop} ifelse }forall
/Encoding ISOLatin1Encoding def
currentdict end
/Courier-ISOL1 exch definefont pop

/Courier-Oblique findfont dup length dict  begin
{1 index /FID ne {def} {pop pop} ifelse }forall
/Encoding ISOLatin1Encoding def
currentdict end
/CourierI-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

/Courier-BoldOblique findfont dup length dict  begin
{1 index /FID ne {def} {pop pop} ifelse }forall
/Encoding ISOLatin1Encoding def
currentdict end
/CourierBI-ISOL1 exch definefont pop

/TG   {350 /Times-ISOL1 FS} bind def
/TM   {280 /Times-ISOL1 FS} bind def
/TK   {240 /Times-ISOL1 FS} bind def
/TS   {200 /Times-ISOL1 FS} bind def
/TL   {160 /Times-ISOL1 FS} bind def

/TIG  {350 /TimesI-ISOL1 FS} bind def
/TIM  {280 /TimesI-ISOL1 FS} bind def
/TIK  {240 /TimesI-ISOL1 FS} bind def
/TIS  {200 /TimesI-ISOL1 FS} bind def
/TIL  {160 /TimesI-ISOL1 FS} bind def

/TBG  {350 /TimesB-ISOL1 FS} bind def
/TBM  {280 /TimesB-ISOL1 FS} bind def
/TBK  {240 /TimesB-ISOL1 FS} bind def
/TBS  {200 /TimesB-ISOL1 FS} bind def
/TBL  {160 /TimesB-ISOL1 FS} bind def

/TBIG {350 /TimesBI-ISOL1 FS} bind def
/TBIM {280 /TimesBI-ISOL1 FS} bind def
/TBIK {240 /TimesBI-ISOL1 FS} bind def
/TBIS {200 /TimesBI-ISOL1 FS} bind def
/TBIL {160 /TimesBI-ISOL1 FS} bind def

/SG {350 /Symbol FS} bind def
/SM {280 /Symbol FS} bind def
/SK {240 /Symbol FS} bind def
/SS {200 /Symbol FS} bind def
/SL {160 /Symbol FS} 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 1 def /lanind2 1 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 /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
 
/cvishow {cvi 6 string cvs show ( ) show} def
/cvsshow1 {10 mul cvi 0.1 mul 9 string cvs show ( ) show} def
/cvsshow2 {100 mul cvi 0.01 mul 9 string cvs show ( ) show} def
/cvsshow3 {1000 mul cvi 0.001 mul 9 string cvs show ( ) show} def
/cvsshow4 {10000 mul cvi 0.0001 mul 9 string cvs show ( ) show} def
/cvsshow5 {100000 mul cvi 0.00001 mul 9 string cvs show ( ) show} def

/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

%calculation of transformation data
/FARBE 27 array def
/FARBE [(R) (m) ( 561_770) %00
        (Y) (m) ( 495_770) %01
        (G) (m) ( 475_573) %02
        (C) (m) ( 380_561) %03
        (B) (m) ( 380_495) %04
        (M) (m) ( 573_475) %05
        (N) (0) ( 380_770) %06
        (W) (0) ( 380_770) %07
        (W) (1) ( 380_770) %08
       ] def

%/FARBE [(R) (m) ( 565_770) %00
%        (Y) (m) ( 495_770) %01
%        (G) (m) ( 475_575) %02
%        (C) (m) ( 380_565) %03
%        (B) (m) ( 380_495) %04
%        (M) (m) ( 575_475) %05
%        (N) (0) ( 380_770) %06
%        (W) (0) ( 380_770) %07
%        (W) (1) ( 380_770) %08
%       ] def

/xshift 10 array def %max 6 device + 4 elementary
/yshift 10 array def
/colori 10 array def

/RX0 640 array def %8x80
/RX0 [

%(%Rm      00 565_770)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
0000 0000 0000 0000 0000 0000 0000 0500 1000 1000

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%(%Ym    01 495_770)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
0000 0000 0000 0500 1000 1000 1000 1000 1000 1000
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%(%Gm      02 475_575)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
0000 0000 0000 0000 0000 0000 0000 0000 0000 0500

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
1000 1000 1000 1000 1000 1000 1000 1000 1000 0500

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%(%Cm     04 380_565)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
1000 1000 1000 1000 1000 1000 1000 0500 0000 0000

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%(%Bm   05 380_495)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
1000 1000 1000 0500 0000 0000 0000 0000 0000 0000
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000

%(%Mm     06 575_475)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
1000 1000 1000 1000 1000 1000 1000 1000 1000 0500

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
0000 0000 0000 0000 0000 0000 0000 0000 0000 0000
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
0000 0000 0000 0000 0000 0000 0000 0000 0000 0500

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%(%N0(r=0,001) 09 380_770)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
0001 0001 0001 0001 0001 0001 0001 0001 0001 0001
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
0001 0001 0001 0001 0001 0001 0001 0001 0001 0001

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
0001 0001 0001 0001 0001 0001 0001 0001 0001 0001
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
0001 0001 0001 0001 0001 0001 0001 0001 0001 0001

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
0001 0001 0001 0001 0001 0001 0001 0001 0001 0001
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
0001 0001 0001 0001 0001 0001 0001 0001 0001 0001

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
0001 0001 0001 0001 0001 0001 0001 0001 0001 0001
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
0001 0001 0001 0001 0001 0001 0001 0001 0001 0001

%(%W1(r=1,000) 10 380_770)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%(%W1(r=1,000) 10 380_770)
%380 %385 %390 %395 %400 %405 %410 %415 %420 %425
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%430 %435 %440 %445 %450 %455 %460 %465 %470 %475
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%480 %485 %490 %495 %500 %505 %510 %515 %520 %525
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%530 %535 %540 %545 %550 %555 %560 %565 %570 %575
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%580 %585 %590 %595 %600 %605 %610 %615 %620 %625
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%630 %635 %640 %645 %650 %655 %660 %665 %670 %675
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

%680 %685 %690 %695 %700 %705 %710 %715 %720 %725
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
%730 %735 %740 %745 %750 %755 %760 %765 %770 %775
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000

] 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

20 setlinewidth
/xdel  1005 def
/ydel   393 def
/x0  050 def
/x1 -150 def

/TELi 8 array def
/TELi [(D65) (D50) (P40) (A00) (E00) (C00) (P00) (Q00)] def %7 illumin.

/Xi 79 array def %max 6 device colour + 4 elementary colours
/Yi 79 array def
/Zi 79 array def
/XANi 79 array def
/YANi 79 array def
/ZANi 79 array def

/xchart4 0 def %
1 1 1 {/xchart4 exch def %xchart4=0,2

/xchartl xchart def
%0 1 7 {/xchartl exch def %xchartl=0,7

/K 0 def %for normalization 1
%/K 1 def %for normalization 0,886

gsave

5 /Times-ISOL1 FS
/cvishow {cvi 6 string cvs show} def
75 85 moveto
lanind cvishow (-) show
colorm cvishow
deintp cvishow
xcolorl cvishow
xchart cvishow
pchart cvishow
colsep cvishow (-L) show pmetam cvishow

5 /Times-ISOL1 FS
340 85 moveto
(SE900-7) show
xchart 0 ge {(N) show} if
(_) show xchart 1 add cvishow
%xchart 0 eq {(N) show}
%            {deintp colorm 2 mul add cvishow} ifelse

72 90 translate

0.010 MM dup scale

20 setlinewidth
0 0 1 0 setcmyk_olvcolor
0 0 moveto 11500     0 rlineto 0 17250 rlineto
           11500 neg 0 rlineto closepath fill
0 0 0 1 setcmyk_olvcolor
0 0 moveto 11500     0 rlineto 0 17250 rlineto
           11500 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


/ymax 17000 def
TBM
50 ymax moveto
(CIE\255Daten f\374r alle Optimalfarben) showde
(CIE data for all optimal colours) showea
TBM
( mit Maximum (m)) showde
( of maximum (m)) showea

TBIM ( C) show 0 -60 rmoveto TBK (AB) show 0 60 rmoveto TBM
(, ) show TELi xchartl get show

( und) showde
( and) showea

TBIM ( Y) show TBK 0 -60 rmoveto (w) show
xcolorl 0 eq {() show}{(,10) show} ifelse 0 60 rmoveto
TBM K 0 eq {(=100) show}{(=88,6) show} ifelse
TBIM (, Y) show
TBK 0 -60 rmoveto (m) show 0 60 rmoveto
TBM
(=495_770) show

0 ymax 100 sub moveto 11500 0 rlineto stroke

/ymaxt ymax 400 sub def

/xpos 100 def
/VtextAB 17 array def

xchart4 0 eq {%xchart4=0
          %0  1  2    3    4    5    6    7    8    9   10
/VtextAB [() () (X)  (Y)  (Z)  (x)  (y)  (z)  (h) ( ) ( )] def
} if %xchart4=0


xchart4 1 eq {%xchart4=1
          %0  1  2    3    4    5    6    7    8   9   10
/VtextAB [() () (Y)  (A)  (B)  (C)  (a)  (b)  (h) ( ) ( )] def
} if %xchart4=1

xchart4 2 eq {%xchart4=2
          %0  1  2    3    4    5    6    7    8    9   10
/VtextAB [() () (L*) (a*) (b*) (C*) (a') (b') (h) ( ) ( )] def
} if %xchart4=2

TBIM
           0 1 09 {/i exch def %i=0,09
                   xpos i 2 ge {x0 add} if
                        xdel i mul add ymaxt moveto
                  i 08 le {VtextAB i get show} if
                  i 02 ge i 04 le and {TBK 0 -60 rmoveto 
                                        K 0 eq {(100) show}
                                               {(88.6) show} ifelse
                                           0  60 rmoveto TBIM} if
                  i 2 ge {%i>=1
                           TBK 0 -60 rmoveto
                           i 05 eq xchart4 1 eq and {(AB) show} if
                           i 05 eq xchart4 2 eq and {(ab) show} if
                           
                           i 08 eq xchart4 0 eq and {(xy) show} if
                           i 08 eq xchart4 1 eq and {(AB) show} if
                           i 08 eq xchart4 2 eq and {(ab) show} if
                           0 60 rmoveto TBIM
                          } if %i>=1
                  } for %i=0,09

xpos 000 add ymaxt moveto
TM (i) show
TBK 0 -60 rmoveto (1) show 0 60 rmoveto
SM (, l) show 
%xchart4 2 eq { TM (*) show TBIM} if
TBK 0 -60 rmoveto (1) show 0 60 rmoveto

xpos xdel add ymaxt moveto
TM (i) show
TBK 0 -60 rmoveto (2) show 0 60 rmoveto
SM (, l) show 
%xchart4 2 eq { TM (*) show TBIM} if
TBK 0 -60 rmoveto (2) show 0 60 rmoveto


xpos x1 add xdel 9 mul add ymaxt moveto
TM (i) show
TBK 0 -60 rmoveto (d) show 0 60 rmoveto
SM (, l) show 
%xchart4 2 eq { TM (*) show TBIM} if
TBK 0 -60 rmoveto (d) show 0 60 rmoveto


xpos x1 add xdel 10 mul add ymaxt moveto
TM (i) show
TBK 0 -60 rmoveto (c) show 0 60 rmoveto
SM (, l) show 
%xchart4 2 eq { TM (*) show TBIM} if
TBK 0 -60 rmoveto (c) show 0 60 rmoveto

xpos x1 add xdel 11 mul add 50 sub ymaxt moveto
(Code) show

/ymax ymax 800 sub def

%      FXU(1)=FF(2)*U(3)-FF(3)*U(2)
%      FXU(2)=FF(3)*U(1)-FF(1)*U(3)
%      FXU(3)=FF(1)*U(2)-FF(2)*U(1)

%      SUBROUTINE DOWEL (F,FXU,IMAX,IO)
%      DIMENSION F(99,3),FXU(3)
%      I=1
%      T2=F(1,1)*FXU(1)+F(1,2)*FXU(2)+F(1,3)*FXU(3)
%      IF(T2) 1,1,4
%    1 T1=T2
%      I=I+1
%      T2=F(I,1)*FXU(1)+F(I,2)*FXU(2)+F(I,3)*FXU(3)
%      IF(T2) 1,1,2
%    2 IF(T1) 3,3,1
%    3 IO=I-1
%      RETURN
%    4 T1=F(IMAX,1)*FXU(1)+F(IMAX,2)*FXU(2)+F(IMAX,3)*FXU(3)
%      IF(T1) 5,5,1
%    5 IO=IMAX
%      RETURN
%      END

%*********************************************************
/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

%*********************************************************

/WEL1 60 array def %irregular definition
/WEL2 60 array def %irregular definition
      %0   1   2   3   4   5   6   7   8   9
/WEL1 [405 435 450 460 465 470 475 480 485 490
       495 500 510 520 530 540 545 550 555 560
      ] def

0 1 19 {/k exch def %k=0,19 %wavelength 405..700nm

        /j  WEL1 k get 400 sub 5 idiv def
        /LP k def
        /LPP 20.5 k add def

       /FF1 FXD j get def %start at 405..700nm
       /FF2 FYD j get def
       /FF3 FZD j get def

       /U1 XW def
       /U2 YW def
       /U3 ZW 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 {/i exch def %i=1,60
               /TN FAD i get FXU1 mul
                   FBD i get FXU2 mul add
                   FCD i get FXU3 mul add def
               /TP FAD i 1 add get FXU1 mul
                   FBD i 1 add get FXU2 mul add
                   FCD i 1 add get FXU3 mul add def
               TN 0 le TP 0 gt and {/INP i def exit} if
              } for %i=1,60

       0 1 60 {/i exch def %i=1,60
               /TN FAD i get FXU1 mul
                   FBD i get FXU2 mul add
                   FCD i get FXU3 mul add neg def
               /TP FAD i 1 add get FXU1 mul
                   FBD i 1 add get FXU2 mul add
                   FCD i 1 add get FXU3 mul add neg def
               TN 0 le TP 0 gt and {/IPN i def exit} if
              } for %i=1,60

TM
xpos xdel  0 mul add ymax ydel LP mul sub moveto INP cvishow %i1
xpos xdel  1 mul add ymax ydel LP mul sub moveto IPN cvishow %i2

xpos xdel  0 mul add ymax ydel LPP mul sub moveto IPN cvishow %i1
xpos xdel  1 mul add ymax ydel LPP mul sub moveto INP cvishow %i2
%*********************************************************
/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
/XIENP XIE 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
/XIEPN XIE def

       } if %kNP>=0

%*********************************************************
%0 1 1 0 setcmyk_olvcolor
INP 0 ge {%INP>=0,-1
          xpos xdel 0 mul add 400 add ymax ydel LP mul sub moveto
          DOWXIE1 0.03 add cvishow %includes rounding failure
          xpos xdel 1 mul add 400 add ymax ydel LPP mul sub moveto
          DOWXIE1 0.03 add cvishow
          /DOWEL1 DOWXIE1 def /DOWEL1c () def
         }%INP>=0
         {%INP=-1
          xpos xdel 0 mul add 400 add ymax ydel LP mul sub moveto
          DOWXIE2 cvishow (c) show
          xpos xdel 1 mul add 400 add ymax ydel LPP 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
          xpos xdel 1 mul add 400 add ymax ydel LP mul sub moveto
          DOWXIE2 0.03 add cvishow
          xpos xdel 0 mul add 400 add ymax ydel LPP mul sub moveto
          DOWXIE2 0.03 add cvishow
          /DOWEL2 DOWXIE2 def /DOWEL2c () def
         }%IPN>=0
         {%IPN=-1
          xpos xdel 1 mul add 400 add ymax ydel LP mul sub moveto
          DOWXIE1 cvishow (c) show
          xpos xdel 0 mul add 400 add ymax ydel LPP mul sub moveto
          DOWXIE1 cvishow (c) show
          /DOWEL2 DOWXIE1 def /DOWEL2c (c) def
         } ifelse %IPN>=0,-1

0 0 0 1 setcmyk_olvcolor

%*********************************************************
0 1 79 {/i exch def %i=0,79
        RA i 0 put
       } for %i=0,79

INP -1 eq {/i1 0 def}                %380nm
          {/i1 INP 5 add def} ifelse %405..700nm

IPN -1 eq {/i2 79 def}               %770nm
          {/i2 IPN 5 add def} ifelse %405..700nm

i2 i1 ge {%i1<i2
          i1 1 i2 {/i exch def %i=i1,i2
                   RA i 1 put
                  } for %i=i1,i2
         }%i1<i2
         {%i1>i2
          i1 1 79 {/i exch def %i=i1,79
                   RA i 1 put
                  } for %i=i1,79
          0  1 i2 {/i exch def %i=0,i2
                   RA i 1 put
                  } for %i=0,i2
                  
         } ifelse %i1>i2

%correction
i2 i1 ge {%i1<i2
          RA i1 0.5 put
          RA i2 0.5 XIEPN 0.5 sub add put
         }%i1<i2
         {%i1>i2
          RA i1 0.5 put
          RA i2 0.5 XIEPN 0.5 sub add put
         } ifelse %i1>i2

%***************************************************
/X 0 def
/Y 0 def
/Z 0 def
          0 1 079 {/i exch def %i=0,79
                   /X X F2N  00 i add get RA i get mul add def
                   /Y Y F2N  80 i add get RA i get mul add def
                   /Z Z F2N 160 i add get RA i get mul add def
                  } for %i=0,79

K 1 eq {%K=1
        /X X 0.8859 mul def
        /Y Y 0.8859 mul def
        /Z Z 0.8859 mul def
       } if %K=1

/SUM X Y add Z add def
/XAN X SUM 0.0001 add div def
/YAN Y SUM 0.0001 add div def
/ZAN Z SUM 0.0001 add div def

%*************************************************
/xpos 100 def

0 1 1 {/Icom exch def %Icom=0,1

Icom 1 eq {%Icom=1
           /X XW X sub def
           /Y YW Y sub def
           /Z ZW Z sub def
           /SUM X Y add Z add def
           /XAN X SUM 0.0001 add div def
           /YAN Y SUM 0.0001 add div def
           /ZAN Z SUM 0.0001 add div def
           /LP 20.5 LP add def
          } if %Icom=1

xpos x1 add xdel 11 mul add ymax ydel LP mul sub moveto
Icom 0 eq {%Icom=0,1
           k  0 eq {(Cm) show} if
           k  6 eq {(Gm) show} if
           k 10 eq {(Ym) show} if
          }%Icom=0
          {%Icom=1
           k  0 eq {(Rm) show} if
           k  6 eq {(Mm) show} if
           k 10 eq {(Bm) show} if
          } ifelse %Icom=0,1

xchart4 0 eq {%xchart4=0

/Hxy YAN YANW sub XAN XANW sub 0.0001 add atan def

xpos x0 add xdel  2 mul add ymax ydel LP mul sub moveto X cvsshow2
xpos x0 add xdel  3 mul add ymax ydel LP mul sub moveto Y cvsshow2
xpos x0 add xdel  4 mul add ymax ydel LP mul sub moveto Z cvsshow2
xpos x0 add xdel  5 mul add ymax ydel LP mul sub moveto XAN cvsshow4
xpos x0 add xdel  6 mul add ymax ydel LP mul sub moveto YAN cvsshow4
xpos x0 add xdel  7 mul add ymax ydel LP mul sub moveto ZAN cvsshow4

xpos x0 add xdel  8 mul add ymax ydel LP mul sub moveto Hxy cvsshow1

} if %xchart4=0

xchart4 1 eq {%xchart4=1

/AAN X Y 0.0001 add div def
/BAN Z Y 0.0001 add div 0.4 mul neg def
/CAN AAN dup mul BAN dup mul add 0.0001 add sqrt def

/A AAN AANW sub Y mul def
/B BAN BANW sub Y mul def
/C A dup mul B dup mul add 0.0001 add sqrt def
/HAB B A 0.0001 add atan def

xpos x0 add xdel  2 mul add ymax ydel LP mul sub moveto Y cvsshow2
xpos x0 add xdel  3 mul add ymax ydel LP mul sub moveto A cvsshow2
xpos x0 add xdel  4 mul add ymax ydel LP mul sub moveto B cvsshow2
xpos x0 add xdel  5 mul add ymax ydel LP mul sub moveto C cvsshow2

xpos x0 add xdel  6 mul add ymax ydel LP mul sub moveto AAN cvsshow4
xpos x0 add xdel  7 mul add ymax ydel LP mul sub moveto BAN cvsshow4

xpos x0 add xdel  8 mul add ymax ydel LP mul sub moveto HAB cvsshow1

} if %xchart4=1

xchart4 2 eq {%xchart4=2

/D13 1 3 div def
/D841_108 841 108 div def
/D16_116 16 116 div def
/D24_116 24 116 div 3 exp def

X XW div D24_116 gt {/FXXN X XW div 0.0001 add D13 exp def}
                     {/FXXN D841_108 X XW div mul D16_116 add def} ifelse

Y YW div D24_116 gt {/FYYN Y YW div 0.0001 add D13 exp def}
                     {/FYYN D841_108 Y YW div mul D16_116 add def} ifelse

Z ZW div D24_116 gt {/FZZN Z ZW div 0.0001 add D13 exp def}
                     {/FZZN D841_108 Z YW div mul D16_116 add def} ifelse

/CIELS 116. FYYN mul 16. sub def
/CIEAS 500. FXXN FYYN sub mul def
/CIEBS 200. FYYN FZZN sub mul def
/CIECS CIEAS dup mul CIEBS dup mul add 0.5 exp def
/CIEHS CIEBS CIEAS 0.000001 add atan def

/a' X Y 0.0001 add div D13 exp  0.2191  mul def
/b' Z Y 0.0001 add div D13 exp -0.08376 mul def
/c' a' dup mul b' dup mul add 0.0001 add sqrt def

xpos x0 add xdel 2 mul add ymax ydel LP mul sub moveto CIELS cvsshow2
xpos x0 add xdel 3 mul add ymax ydel LP mul sub moveto CIEAS cvsshow2
xpos x0 add xdel 4 mul add ymax ydel LP mul sub moveto CIEBS cvsshow2
xpos x0 add xdel 5 mul add ymax ydel LP mul sub moveto CIECS cvsshow2

xpos x0 add xdel 6 mul add ymax ydel LP mul sub moveto a' cvsshow4
xpos x0 add xdel 7 mul add ymax ydel LP mul sub moveto b' cvsshow4

xpos x0 add xdel 8 mul add ymax ydel LP mul sub moveto CIEHS cvsshow1

} if %xchart4=2

%*********************************************************
%calculation of lamdad and lambdac for X, Y, Z
       /FF1 X def
       /FF2 Y def
       /FF3 Z def

       /U1 XW def
       /U2 YW def
       /U3 ZW 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 {/i exch def %i=1,60
               /TN FAD i get FXU1 mul
                   FBD i get FXU2 mul add
                   FCD i get FXU3 mul add def
               /TP FAD i 1 add get FXU1 mul
                   FBD i 1 add get FXU2 mul add
                   FCD i 1 add get FXU3 mul add def
               TN 0 le TP 0 gt and {/INP i def exit} if
              } for %i=1,60

       0 1 60 {/i exch def %i=1,60
               /TN FAD i get FXU1 mul
                   FBD i get FXU2 mul add
                   FCD i get FXU3 mul add neg def
               /TP FAD i 1 add get FXU1 mul
                   FBD i 1 add get FXU2 mul add
                   FCD i 1 add get FXU3 mul add neg def
               TN 0 le TP 0 gt and {/IPN i def exit} if
              } for %i=1,60

TM
xpos x1 add xdel  9 mul add ymax ydel LP mul sub moveto INP cvishow
xpos x1 add xdel 10 mul add ymax ydel LP 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
          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
          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
          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
          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

0 0 0 1 setcmyk_olvcolor

Icom 0 eq {%Icom=0,1
           Xi k X put
           Yi k Y put
           Zi k Z put
           XANi k XAN put
           YANi k YAN put
           ZANi k ZAN put
          }%Icom=0
          {%Icom=1
           /kP 20 k add def
           Xi kP X put
           Yi kP Y put
           Zi kP Z put
           XANi kP XAN put
           YANi kP YAN put
           ZANi kP ZAN put
          } ifelse %Icom=1

} for %k=0,19

} for %Icom=0,1

%**********************************************************************
%XW,YW,ZW and XK,YK,ZK 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

K 0 eq {%K=0
        /XK XW def
        /YK YW def
        /ZK ZW def
       }%K=0
       {%K=1
        /XK XW 0.8859 mul def
        /YK YW 0.8859 mul def
        /ZK ZW 0.8859 mul def
       } ifelse %K=1

/SUMK XK YK add ZK add def %for D65, D50, P40, A00, E00
/XANK XK SUMK div def
/YANK YK SUMK div def
/ZANK ZK SUMK div def

/LP LP 1.5 add def
xpos xdel 0 mul add 400 add ymax ydel LP mul sub moveto (380) show
xpos xdel 1 mul add 400 add ymax ydel LP mul sub moveto (770) show

xchart4 0 eq {%xchart4=0

/HxyK YANK YANK sub XANK XANK sub 0.0001 add atan def

xpos x0 add xdel  2 mul add ymax ydel LP mul sub moveto XK cvsshow2
xpos x0 add xdel  3 mul add ymax ydel LP mul sub moveto YK cvsshow2
xpos x0 add xdel  4 mul add ymax ydel LP mul sub moveto ZK cvsshow2
xpos x0 add xdel  5 mul add ymax ydel LP mul sub moveto XANK cvsshow4
xpos x0 add xdel  6 mul add ymax ydel LP mul sub moveto YANK cvsshow4
xpos x0 add xdel  7 mul add ymax ydel LP mul sub moveto ZANK cvsshow4

xpos x0 add xdel  8 mul add ymax ydel LP mul sub moveto HxyK cvsshow1

} if %xchart4=0

xchart4 1 eq {%xchart4=1

/AANK XK YK 0.0001 add div def
/BANK ZK YK 0.0001 add div 0.4 mul neg def
/CANK AANK dup mul BANK dup mul add 0.0001 add sqrt def

/AK AANK AANK sub YK mul def
/BK BANK BANK sub YK mul def
/CK AK dup mul BK dup mul add 0.0001 add sqrt def
/HABK BK AK 0.0001 add atan def

xpos x0 add xdel  2 mul add ymax ydel LP mul sub moveto YK cvsshow2
xpos x0 add xdel  3 mul add ymax ydel LP mul sub moveto AK cvsshow2
xpos x0 add xdel  4 mul add ymax ydel LP mul sub moveto BK cvsshow2
xpos x0 add xdel  5 mul add ymax ydel LP mul sub moveto CK cvsshow2

xpos x0 add xdel  6 mul add ymax ydel LP mul sub moveto AANK cvsshow4
xpos x0 add xdel  7 mul add ymax ydel LP mul sub moveto BANK cvsshow4

xpos x0 add xdel  8 mul add ymax ydel LP mul sub moveto HABK cvsshow1

} if %xchart4=1

xchart4 2 eq {%xchart4=2

/D13 1 3 div def
/D841_108 841 108 div def
/D16_116 16 116 div def
/D24_116 24 116 div 3 exp def

XK XW div D24_116 gt {/FXXNK XK XW div 0.0001 add D13 exp def}
                     {/FXXNK D841_108 XK XW div mul D16_116 add def} ifelse

YK YW div D24_116 gt {/FYYNK YK YW div 0.0001 add D13 exp def}
                     {/FYYNK D841_108 YK YW div mul D16_116 add def} ifelse

ZK ZW div D24_116 gt {/FZZNK ZK ZW div 0.0001 add D13 exp def}
                     {/FZZNK D841_108 ZK YW div mul D16_116 add def} ifelse

/CIELSK 116. FYYNK mul 16. sub def
/CIEASK 500. FXXNK FYYNK sub mul def
/CIEBSK 200. FYYNK FZZNK sub mul def
/CIECSK CIEASK dup mul CIEBSK dup mul add 0.5 exp def
/CIEHSK CIEBSK CIEASK 0.000001 add atan def

/a'K XK YK 0.0001 add div D13 exp  0.2191  mul def
/b'K ZK YK 0.0001 add div D13 exp -0.08376 mul def
/c'K a'K dup mul b'K dup mul add 0.0001 add sqrt def

xpos x0 add xdel 2 mul add ymax ydel LP mul sub moveto CIELSK cvsshow2
xpos x0 add xdel 3 mul add ymax ydel LP mul sub moveto CIEASK cvsshow2
xpos x0 add xdel 4 mul add ymax ydel LP mul sub moveto CIEBSK cvsshow2
xpos x0 add xdel 5 mul add ymax ydel LP mul sub moveto CIECSK cvsshow2

xpos x0 add xdel 6 mul add ymax ydel LP mul sub moveto a'K cvsshow4
xpos x0 add xdel 7 mul add ymax ydel LP mul sub moveto b'K cvsshow4

xpos x0 add xdel 8 mul add ymax ydel LP mul sub moveto CIEHSK cvsshow1

} if %xchart4=2


showpage
grestore

%} for %xchartl=0,7

} for %xchart4=0,2

%} for %xcolorl=0,1

%%Trailer