%!PS-Adobe-3.0 EPSF-3.0 http://farbe.li.tu-berlin.de/hep8/hep80-3n
%%BoundingBox: 70 85 421 335

%START PDFDE011.EPS
/pdfmark03 where {pop} {userdict /pdfmark03 /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: farbe.li.tu-berlin.de/hep8/hep8.HTM)
 /Author (compare K. Richter "Computergrafik ...": ISBN 3-8007-1775-1)
 /Subject (goto: http://farbe.li.tu-berlin.de or http://color.li.tu-berlin.de)
 /Keywords (image reproduction, colour devices)
 /Creator (klaus.richter@mac.com)
 /CreationDate (D:2024100112000)
 /ModDate (D:20241001112000)
/DOCINFO pdfmark03
[ /View [ /Fit ]
/DOCVIEW pdfmark03
%END  PDFDE011

/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

/FS {findfont exch scalefont setfont} bind def
/MM {72 25.4 div mul} def /str {8 string } bind def

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

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

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

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

/SS {160 /Symbol FS} bind def
/SL {200 /Symbol FS} bind def
/SK {250 /Symbol FS} bind def
/SM {300 /Symbol FS} bind def
/SG {350 /Symbol FS} bind def

/CS {160 /Courier-ISOL1 FS} bind def
/CL {200 /Courier-ISOL1 FS} bind def
/CK {250 /Courier-ISOL1 FS} bind def
/CM {300 /Courier-ISOL1 FS} bind def
/CG {350 /Courier-ISOL1 FS} bind def

/CBS {160 /CourierB-ISOL1 FS} bind def
/CBL {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

/nGs {350 /Times-ISOL1 FS  show} bind def
/kGs {350 /TimesI-ISOL1 FS  show} bind def
/bGs {350 /TimesB-ISOL1 FS  show} bind def
/jGs {350 /TimesBI-ISOL1 FS  show} bind def
/sGs {350 /Symbol FS  show} bind def
/iGs {300 /Times-ISOL1 FS 0 -80 rmoveto show 0   80 rmoveto} bind def
/eGs {300 /Times-ISOL1 FS 0 200 rmoveto show 0 -200 rmoveto} bind def
/ibGb {300 /TimesB-ISOL1 FS 0 -80 rmoveto show 0   80 rmoveto} bind def
/ebGb {300 /TimesB-ISOL1 FS 0 200 rmoveto show 0 -200 rmoveto} bind def
/ipG {300 /Times-ISOL1 FS 50 50 rmoveto (\267) show 50 -50 rmoveto} bind def

%20% kleiner
/nMs {300 /Times-ISOL1 FS  show TM} bind def
/kMs {300 /TimesI-ISOL1 FS  show TM} bind def
/bMs {300 /TimesB-ISOL1 FS  show TM} bind def
/jMs {300 /TimesBI-ISOL1 FS  show TM} bind def
/sMs {300 /Symbol FS  show TM} bind def
/iMs {250 /Times-ISOL1 FS 0 -60 rmoveto show 0 60 rmoveto TM} bind def
/eMs {250 /Times-ISOL1 FS 0 160 rmoveto show 0 -160 rmoveto TM} bind def
/ibMs {250 /TimesB-ISOL1 FS 0 -60 rmoveto show 0 60 rmoveto TM} bind def
/ebMs {250 /TimesB-ISOL1 FS 0 160 rmoveto show 0 -160 rmoveto TM} bind def
/ipM {250 /Times-ISOL1 FS 40 40 rmoveto (\267) show 40 -40 rmoveto TM} bind def

%40% kleiner
/nKs {250 /Times-ISOL1 FS  show TK} bind def
/kKs {250 /TimesI-ISOL1 FS  show TK} bind def
/bKs {250 /TimesB-ISOL1 FS  show TK} bind def
/jKs {250 /TimesBI-ISOL1 FS  show TK} bind def
/sKs {250 /Symbol FS  show TK} bind def
/iKs {200 /Times-ISOL1 FS 0 -50 rmoveto show 0 50 rmoveto TK} bind def
/eKs {200 /Times-ISOL1 FS 0 130 rmoveto show 0 -130 rmoveto TK} bind def
/ibKs {200 /TimesB-ISOL1 FS 0 -50 rmoveto show 0 50 rmoveto TK} bind def
/ebKs {200 /TimesB-ISOL1 FS 0 130 rmoveto show 0 -130 rmoveto TK} bind def
/ipK {200 /Times-ISOL1 FS 30 30 rmoveto (\267) show 30 -30 rmoveto TK} bind def

%60% kleiner
/nLs {200 /Times-ISOL1 FS  show TL} bind def
/kLs {200 /TimesI-ISOL1 FS  show TL} bind def
/bLs {200 /TimesB-ISOL1 FS  show TL} bind def
/jLs {200 /TimesBI-ISOL1 FS  show TL} bind def
/sLs {200 /Symbol FS  show TL} bind def
/iLs {160 /Times-ISOL1 FS 0 -40 rmoveto show 0 40 rmoveto TL} bind def
/eLs {160 /Times-ISOL1 FS 0 110 rmoveto show 0 -110 rmoveto TL} bind def
/ibLs {160 /TimesB-ISOL1 FS 0 -40 rmoveto show 0 40 rmoveto TL} bind def
/ebLs {160 /TimesB-ISOL1 FS 0 110 rmoveto show 0 -110 rmoveto TL} bind def
/ipL {160 /Times-ISOL1 FS 25 25 rmoveto (\267) show 25 -25 rmoveto TL} bind def
/jbLs {160 /TimesBI-ISOL1 FS 0 110 rmoveto show 0 -110 rmoveto TL} bind def

%80% smaller
/nSs {160 /Times-ISOL1 FS  show TS} bind def
/kSs {160 /TimesI-ISOL1 FS  show TS} bind def
/bSs {160 /TimesB-ISOL1 FS  show TS} bind def
/jSs {160 /TimesBI-ISOL1 FS  show TS} bind def
/sSs {160 /Symbol FS  show TS} bind def
/iSs {130 /Times-ISOL1 FS 0 -30 rmoveto show 0 30 rmoveto TS} bind def
/eSs {130 /Times-ISOL1 FS 0 80 rmoveto show 0 -80 rmoveto TS} bind def
/ibSs {130 /TimesB-ISOL1 FS 0 -30 rmoveto show 0 30 rmoveto TS} bind def
/ebSs {130 /TimesB-ISOL1 FS 0 80 rmoveto show 0 -80 rmoveto TS} bind def
/ipS {130 /Times-ISOL1 FS 20 20 rmoveto (\267) show 20 -20 rmoveto TS} bind def
/jbSs {130 /TimesBI-ISOL1 FS 0 80 rmoveto show 0 -80 rmoveto TS} bind def

/20rm {20 0 rmoveto} def

/cvishow {cvi 10 string cvs show} def
/cvsshow1 {10 mul cvi 0.1 mul 10 string cvs show} def
/cvsshow2 {100 mul cvi 0.01 mul 10 string cvs show} def
/cvsshow3 {1000 mul cvi 0.001 mul 10 string cvs show} def
/cvsshow4 {10000 mul cvi 0.0001 mul 10 string cvs show} def

/cvsshow1x {/nxx exch def                              %example nxx=99.1/99.0
             nxx 0 lt {(-) show}
                      {() show} ifelse
             nxx 10 mul cvi abs /nxi exch def          %nxi=991/990
             nxi 10 idiv /nxa exch def                 %nxa=99
             nxi nxa 10 mul sub /nxb exch def          %nxb=1/0
             nxa cvishow (,) show                      %nxa=99
             nxb cvishow                               %nxb=1/0
           } def
/cvsshow2x {/nxx exch def                              %example nxx=99.12/99,02/99,00
             nxx 0 lt {(-) show}
                      {() show} ifelse
             nxx 100 mul cvi abs /nxi exch def         %nxi=9912/9902/00
             nxi 100 idiv /nxa exch def                %nxa=99
             nxi nxa 100 mul sub /nxb exch def         %nxb=12/02/00
             nxa cvishow (,) show                      %nxb=99,
             nxb 10 ge {nxb cvishow} if                %nxb=12
             nxb  1 ge
             nxb  9 le and {(0) show nxb cvishow} if   %nxb=02
             nxb  0 eq {(00) show} if                  %nxb=00
           } def
/cvsshow3x {/nxx exch def                              %example nxx=99.123/99.012/99.001/99.000
             nxx 0 lt {(-) show}
             {() show} ifelse
             nxx 1000 mul cvi abs /nxi exch def        %nxi=99123/99012/99001/99000
             nxi 1000 idiv /nxa exch def               %nxa=99
             nxi nxa 1000 mul sub /nxb exch def        %nxb=123/012/001/000
             nxa cvishow (,) show                      %nxb=99,
             nxb 100 ge {nxb cvishow} if               %nxb=123/012/001/000
             nxb  10 ge
             nxb  99 le and {(0) show nxb cvishow} if  %nxb=012
             nxb   1 ge
             nxb   9 le and {(00) show nxb cvishow} if %nxb=001
             nxb   0 eq {(000) show} if                %nxb=000
           } def

/cvsshow4x {/nxx exch def                              %example nxx=99.123/99.0
             nxx 0 lt {(-) show}
                      {() show} ifelse
             nxx 10000 mul cvi abs /nxi exch def        %nxi=99123/99012/99001/9
             nxi 10000 idiv /nxa exch def               %nxa=99
             nxi nxa 10000 mul sub /nxb exch def        %nxb=123/012/001/000
             nxa cvishow (,) show                       %nxb=99,
             nxb 1000 ge {nxb cvishow} if               %nxb=123/012/001/000
             nxb  100 ge
             nxb  999 le and {(0) show nxb cvishow} if  %nxb=012
             nxb  10 ge
             nxb  99 le and {(00) show nxb cvishow} if  %nxb=012
             nxb   1 ge
             nxb   9 le and {(000) show nxb cvishow} if %nxb=001
             nxb   0 eq {(0000) show} if                %nxb=000
           } def

/cvsshow4s {/nxx exch def                              %example nxx=-0,1234
             nxx 0 lt {(-0,) show}
                      {(0,) show} ifelse
             /nxi nxx 10000 mul cvi abs def            %nxi=1234
             nxi 1000 ge {nxi cvishow} if              %nxb=123/012/001/000
             nxi 100 ge
             nxi 999 le and {(0) show nxi cvishow} if      %nxb=123/012/001/000
             nxi  10 ge
             nxi  99 le and {(00) show nxi cvishow} if  %nxb=012
             nxi   1 ge
             nxi   9 le and {(000) show nxi cvishow} if %nxb=001
             nxi   0 eq {(0000) show} if                %nxb=000
           } def

%XCHA01.PS BEG
/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 r g b
  {setrgbcolor rec fill} bind def

/colrecst %x y width heigth r g b
  {setrgbcolor rec stroke} bind def

/rem %x, y width heigth
  {/heigth exch 0.5 mul def /width exch 0.5 mul def
   /yleftb exch heigth 0.5 mul add def
   /xleftb exch width  0.5 mul add def
   xleftb yleftb
   moveto width 0 rlineto 0 heigth rlineto width neg 0 rlineto
   closepath } bind def

/colremfi %x y width heigth r g b
  {setrgbcolor rem fill} bind def

/colremst %x y width heigth r g b
  {setrgbcolor rem stroke} bind def

/tfr {1.0 0.0 0.0 setrgbcolor} bind def %Reproduktionsfarben
/tfg {0.0 1.0 0.0 setrgbcolor} bind def
/tfb {0.0 0.0 1.0 setrgbcolor} bind def
/tfc {0.0 1.0 1.0 setrgbcolor} bind def
/tfm {1.0 0.0 1.0 setrgbcolor} bind def
/tfy {1.0 1.0 0.0 setrgbcolor} bind def

/tfw {1.00 1.00 1.00 setrgbcolor} bind def %Graureihe
/tfh {0.75 0.75 0.75 setrgbcolor} bind def
/tfz {0.50 0.50 0.50 setrgbcolor} bind def
/tfd {0.25 0.25 0.25 setrgbcolor} bind def
/tfn {0.00 0.00 0.00 setrgbcolor} bind def

/YnW 100 def %HAULAB, CIELAB

/Yi   101 array def %101(W)+3*101(R,G,B)
/dYi  101 array def
/L*i  101 array def
/logL*i 101 array def
/X00i 101 array def %log(Xi)

/Y00i 101 array def %L*i, log(L*i)
/Y0ui 101 array def %L*i/Lu, log(L*i/L*u)
/Y10i 101 array def %dYi, log(dYi)
/Y1ui 101 array def %dYi/dYu, log(dYi/dYu)
/Y20i 101 array def %dYi/Yi, log(dYi/Yi) sensitivity
/Y2ui 101 array def %(dYi/dYu)/(Yi/Yu), log[(dYi/dYu)/(Yi/Yu)]
/Y30i 101 array def %Yi/dYi, log(Yi/dYi) contrast
/Y3ui 101 array def %(Yi/Yu)/(dYi/dYu), log[(Yi/Yu)/(dYi/dYu)]

/Yx0i 101 array def %one of four Y00i, Y10i, Y20i, Y20i
/Yxui 101 array def %one of four Y0ui, Y1ui, Y2ui, Y2ui

/MULX 1000 def
/MULY 1000 def

/ifunc 0 def %0:HAULAB, 1:CIELAB 2:IECsRGB 3:TUBsRGB
ifunc 0 eq {/i1 002 def}        %HauLAB
           {/i1 001 def} ifelse %all others
/i2 100 def

%data for HAULAB IECsRGB, TUBsRGB
/c32 3.2258 def %HAULAB n=0.31
/e10D32 1.0 3.2258 div def
/e20D32 2.2258 3.2258 div def

/c24 2.4 def %IECsRGB
/e10D24 1.0 2.4 div def
/e14D24 1.4 2.4 div def

/c30 3.0 def %CIELAB
/e10D30 1.0 3.0 div def
/e20D30 2.0 3.0 div def

/c23 2.3 def %TUBsRGB
/e10D23 1.0 10 ln div def
/e13D23 10 ln 1 sub 10 ln div def

/econst 2.71828182 def
/W2 2 sqrt def
/FL 0.0001 def

/x00t 0400 def %xpos for BEG equations
/x01t 1900 def %xpos for shift equations
/x00e 5250 def %xpos for Num equations

%STOP0
%***************************************************
ifunc 0 eq {%ifunc=0 BEG HAULAB
%2: 1/3.2258=0.3100
%for Phi=120 in he60/he60-3a.eps
%L* =134.60*(Y/Yn)**(1/3.2258)-34.60
%   =134.60*(Yu/Yn)**0.31*(Y/Yu)**0.31-34.60
%   =134.60*(18/100)**0.31*(Y/Yu)**0.31-34.60

%L*u=134.60*0.5876        *(Y/Yu)**0.31-34.60
%   =79.09                             -34.60
%   =45.39
%
%s*(Yu/Yn)**0.31=r*(Yu/Yu)**0.31=1
%r=s*(Yu/Yn)**0.31
% =134.60*(0.18)**0.31
% =134.60.5876
% =79.09
%
/Yn 100 def
/L*u 50 def
/Yu L*u 34.60 add 134.60 div 3.2258 exp 100 mul def
/dYu Yn e10D32 exp 134.60 div 3.2258 mul Yu e20D32 exp mul def
/iu 18 def

/aCIE 3.2258 134.60 div Yn e20D32 exp mul def
/bCIE aCIE iu e20D30 exp mul def

/cCIE 3.2258 134.60 div Yn e10D32 exp mul def
/dCIE cCIE iu e20D30 exp mul def

/eCIE 134.60 3.2258 div Yn e20D32 exp mul def
/fCIE eCIE iu e20D32 exp mul def

%*************************************************
/Yi_L*i_dYi_H_0 {%BEG Yi_L*i_dYi_H_0 H=HAULAB
i1 1 100 {/i exch def %i=1,100
         Yi  i i YnW mul 100. div put
         L*i i Yi i get YnW div e10D32 exp 134.60 mul 34.60 sub put
         dYi i Yi i get YnW div e20D32 exp c32 mul
         100 mul 134.60 div put
         logL*i i L*i i get log put
        } for %i=1,100
%STOPA
} bind def %END Yi_L*i_dYi_H_0 C=HAULAB

} if %ifunc=0 END HAULAB

%STOP
%***************************************************
ifunc 1 eq {%ifunc=1 BEG CIELAB
%2: 1/2,4=0.41667
%L*u=116*(Yu/Yn)**(1/3)-16
%   =116*(0.18)**(1/3) -16
%   =116*0.5656-16
%   =65.50-16=49,50
%
%(L*u+16)/116=(Yu/Yn)**(1/3)
%Yu=Yn*(L*u+16)/116)**3
%Yu=100*(65.50/116)**3
%  =100*0,5647**3
%  =100*0,1800=18.00
%Yn=100, Yu=18 L*u=49,50
%
%s*(Yu/Yn)**(1/3)=r*(Yu/Yu)**(1/3)=1
%r=s*(Yu/Yn)**(1/3)
% =116*(0.18)**(1/3)
% =116+0.5656
% = 65.50
%
/Yn 100 def
/L*u 50 def
/Yu L*u 16 add 116 div 3 exp 100 mul def
/dYu Yn e10D30 exp 116 div 3 mul Yu e20D30 exp mul def
/iu 18 def

/aCIE 3 116 div Yn e20D30 exp mul def
/bCIE aCIE iu e20D30 exp mul def

/cCIE 3 116 div Yn e10D30 exp mul def
/dCIE cCIE iu e20D30 exp mul def

/eCIE 116 3 div Yn e20D30 exp mul def
/fCIE eCIE iu e20D30 exp mul def

%*************************************************
/Yi_L*i_dYi_C_0 {%BEG Yi_L*i_dYi_C_0 C=CIELAB
i1 1 100 {/i exch def %i=1,100
         Yi  i i YnW mul 100. div put
         L*i i Yi i get YnW div e10D30 exp 116 mul 16 sub put
         dYi i Yi i get YnW div e20D30 exp c30 mul
         100 mul 116 div put
        } for %i=1,100
} bind def %END Yi_L*i_dYi_C_0 C=CIELAB

} if %ifunc=1 END CIELAB

%***************************************************
ifunc 2 eq {%ifunc=2 BEG IECsRGB

%L*=100(Y/Yn)**(1/2.4)
%  =100*(Yu/Yn)**(1/2.4)*(Y/Yu)**(1/2,4)
%  =g                   *(Y/Yu)**(1/2,4)
% g=100*(18/100)**(1/2,4)=100*(0,18)**0,4166
% g=48,95
 
%L*u=100(Yu/Yn)**(1/2.4)
%L*/L*u=(Y/Yu)**(1/2.4)
%log[L*/L*u]=(1/2,4)*log(Y/Yu)=0,4166*log(Y/Yu)
%ln [L*/L*n]=2.30258*0.4166*log(Y/Yu)=0,9593
%oder
%log[L*/L*u]=(1/2,3)*log(Y/Yu)=0,4347*log(Y/Yu)
%ln [L*/L*u]=2.30258*0.4347*log(Y/Yu)=1,001*log(Y/Yu)
%
%/Yn 100 def
%/L*u 50 def

%1: not used 50=100(Yu/100)**(1/2.4)
%0.5**(2.4)=(Yu/100)
%Yu=100*0.5**2.4=18.94

%2: 1/2,4=0.41667
%L*u=100(Yu/Yn)**(1/2.4)
%L*u=100(18/100)**(1/2.4)=48.94
%Yu=Yn*(Lu/100)**2.4
%Yu=100(Lu/100)**2.4
%Yn=100, Yu=18 L*u=48.95

%L*=100*(Y/Yn)**(1/2,4)
%dL*/dY=100*(1/2,4)*(1/Yn)*(Y/Yn)**(-1,4/2,4)
%fuer dL*=1:
%dY =(2,4*Yn)/100        *(Y/Yn)**(1,4/2,4)
%   =a                   *(Y/Yn)**(1,4/2,4)
%   =a*(Yu/Yn)**(1,4/2,4)*(Y/Yu)**(1,4/2,4)
%   =b                   *(Y/Yu)**(1,4/2,4)
%   =2,4*(Y/Yn)**(1.4/2,4)
%   =2,4*(1/Yn)**(1,4/2,4)*(Y)**(1,4/2,4)
%dYu=2,4*(Yu/Yn)**(1,4/2,4)
%   =2,4*(1/Yn)**(1,4/2,4)*(Yu)**(1,4/2,4)
%dY/dYu=(Y/Yu)**(1,4/2,4)

%dY/Y=2,4*(Y/Yn)**(1,4/2,4)*(1/Y)
%    =2,4*(1/Yn)**(1,4/2,4)*(Y)**(1,4/2,4)*(1/Y)
%    =2,4*(1/Yn)**(1,4/2,4)*(Y)**(-1/2,4)
%    =c*Y**(-1/2,4)

%a=(2,4*Yn)/100
% =2,4
%b=a*(Yu/Yn)**(1,4/2,4)
% =2,4*(18/100)**(1,4/2,4)
% =2,4*(18/100)**(0,583333)
% =2,4*0,36777=0,8862
%c=2,4(1/Yn)**(1,4/2,4)
% =2,4*0.01**0,5833=0,1635

/Yn 100 def
/L*u 50 def
/Yu L*u 100 div 2.4 exp 100 mul def
/dYu Yu 100 div e14D24 exp 2.4 mul def
/iu 18 def

/aCIE 2.4 def
/bCIE iu 100 div e14D24 exp aCIE mul def

/cCIE 2.4 100 div Yn e10D24 exp mul def
/dCIE cCIE iu e14D24 exp mul def

/eCIE 100 2.4 div Yn e14D24 exp mul def
/fCIE eCIE iu e14D24 exp mul def

%*************************************************
/Yi_L*i_dYi_I_0 {%BEG Yi_L*i_dYi_I_0 I=IECsRGB
i1 1 100 {/i exch def %i=1,100
         Yi  i i YnW mul 100. div put
         L*i i Yi i get YnW div e10D24 exp 100 mul put
         dYi i Yi i get YnW div e14D24 exp c24 mul
         100 mul 100 div put
        } for %i=1,100
} bind def %END Yi_L*i_dYi_I_0 I=IECsRGB

} if %ifunc=2 END IECsRGB

%***************************************************
ifunc 3 eq {%ifunc=3 BEG TUBsRGB

%L*=100(Y/Yn)**(1/2.3)
%L*u=100(Yu/Yn)**(1/2.3)
%/Yn 100 def
%/L*u 50 def

%1: not used 50=100(Yu/100)**(1/2.3)
%0.5**(2.3)=(Yu/100)
%Yu=100*0.5**2.3=20.31

%2: 1/2,3=0.4348
%L*u=100(Yu/Yn)**(1/2.3)
%L*u=100(18/100)**(1/2.3)=47.45
%Yu=Yn*(Lu/100)**2.3
%Yu=100(Lu/100)**2.3
%Yn=100, Yu=18 L*u=47.45

%L*=100*(Y/Yn)**(1/2,3)
%dL*/dY=100*(1/2,3)*1/Yn)*(Y/Yn)**(-1,3/2,3)
%dY=(2.3*Yn)/100*(Y/Yn)**(1,3/2,3)

%dYu=2.3*(Yu/Yn)**(1,3/2,3)

/Yn 100 def
/L*u 50 def

/Yu L*u 100 div 2.3 exp 100 mul def
/dYu Yu 100 div e13D23 exp 2.3 mul def
/iu 18 def

/aCIE 2.3 100 div Yn e13D23 exp mul def
/bCIE aCIE iu e13D23 exp mul def

/cCIE 2.3 100 div Yn e10D23 exp mul def
/dCIE cCIE iu e13D23 exp mul def

/eCIE 100 2.3 div Yn e13D23 exp mul def
/fCIE eCIE iu e13D23 exp mul def

%*************************************************
/Yi_L*i_dYi_T_0 {%BEG Yi_L*i_dYi_T_0 T=TUBsRGB
i1 1 100 {/i exch def %i=1,100
         Yi  i i YnW mul 100. div put
         L*i i Yi i get YnW div e10D23 exp 100 mul put
         dYi i Yi i get YnW div e13D23 exp c23 mul put
        } for %i=1,100
} bind def %END Yi_L*i_dYi_T_0 T=TUBsRGB

} if %ifunc=3 END TUBsRGB

%***********************************************
/ioute 0 def
/cero_line {%BEG cero_line
%1. log line = cero line
tfr
-2. MULX mul 1. log MULY mul moveto
 0. MULX mul 1. log MULY mul lineto stroke
tfn %tfw

%2. log line = +1 line
tfg
-2. MULX mul 10. log MULY mul moveto
 0. MULX mul 10. log MULY mul lineto stroke
tfn %tfw
} bind def %cero line

%***************************************************
/Yi_Yxyi {%BEG Yi_Yxyi definition for 4 functions: L*i, dYi, dYi/Yi, Yi/dYi
%i=0 not available
%no log, 1 1 100 allways possible
i1 1 100 {/i exch def %i=1,100

         X00i i Yi  i get put

         Y00i i L*i i get put
         Y0ui i L*i i get L*u div put
         Y10i i dYi i get put
         Y1ui i dYi i get dYu div put
         Y20i i dYi i get Yi i get div put
         Y2ui i dYi i get Yi i get div
                dYu       Yu       div div put
         Y30i i Yi i get dYi i get div put
         Y3ui i Yi i get dYi i get div
                Yu       dYu       div div put

         xchartl 0 eq {Yx0i i Y00i i get put
                       Yxui i Y0ui i get put} if
         xchartl 1 eq {Yx0i i Y10i i get put
                       Yxui i Y1ui i get put} if
         xchartl 2 eq {Yx0i i Y20i i get put
                       Yxui i Y2ui i get put} if
         xchartl 3 eq {Yx0i i Y30i i get put
                       Yxui i Y3ui i get put} if
        } for %i=1,100

} bind def %END Yi_Yxyi

%***************************************************
/Y_curve {%BEG Y_curve
50 setlinewidth
0 1 1 {/je exch def %je=0,1
je 0 eq {1 1 1 setrgbcolor} if
je 1 eq {0 0 0 setrgbcolor [100] 0 setdash} if
i1 1 i2 {/i exch def %i=i1,i2
         X00i i get log MULX mul
         Yxui i get jlog 1 eq {log} if MULY mul
         i i1 eq {moveto} if
         i i1 1 add ge
         i i2 1 sub le and {lineto} if
         i i2 eq {stroke} if
        } for %i=i1,i2
} for %je=0,1

50 setlinewidth
0 1 1 {/ki1i2 exch def %ki1i2=0,1 
       ki1i2 0 eq {/i i1 def} if
       ki1i2 1 eq {/i i2 def} if
       tfb
       X00i i get log MULX mul
       Yxui i get jlog 1 eq {log} if MULY mul 060 0 360 arc fill
       newpath
       X00i i get log i 1 eq {0.80 sub}{0.10 add} ifelse MULX mul
       Yxui i get jlog 1 eq {log} if 0.10 sub MULY mul moveto
       TBK
       Yxui i get jlog 1 eq {log} if cvsshow3x
       0 setgray
       newpath
      } for %ki1i2=0,1

[ ] 0 setdash
0 setgray

} bind def %END Y_curve

%***************************************************
/proc_appli {%proc_appli
/x00a 4300 def
%y00a is to be defined in main program
tfn %tfw
x00a y00a moveto
1000 0 rlineto stroke
tfg
[100] 0 setdash
x00a y00a moveto
1000 0 rlineto stroke
[ ] 0 setdash

tfn %tfw
x00a y00a 0.8 ydel mul sub moveto
TBL (application) showen
(Anwendungs\255) showde
x00a y00a 1.5 ydel mul sub moveto
TBL (range) showen
(bereich) showde
} bind def %proc_appli

%***************************************************
/proc_toptext {%proc_toptext
%jlog 0:without log, 1:with log in main program
tfn %tfw
/ytr0t0 3750 ytr0 sub def
/ytr0t1 ytr0t0 250 sub def
/ytr0t2 ytr0t0 500 sub def
 
050 xtr0 sub 3725 ytr0 sub moveto
xchartl 00 eq {%xchartl=01
               jlog 1 eq {(log \050) bKs} if 
               (L*) jKs  ifunc 0 eq {(80) ibKs} if
               (/L*) jKs ifunc 0 eq {(80,) ibKs} if (u) ibKs
               jlog 1 eq {(\051 ) bKs} if

               2100 xtr0 sub ytr0t0 moveto
               TBK ifunc 0 eq {(HAULAB) show} if
                   ifunc 1 eq {(CIELAB) show} if
                   ifunc 2 eq {(IECsRGB) show} if
                   ifunc 3 eq {(TUBsRGB) show} if
               TBK ( lightness ) showen
               (\255Helligkeit ) showde
               (L*) jKs ifunc 0 eq {(80) ibKs} if
               TBK ( normalized) showen
               ( normiert) showde

               2100 xtr0 sub ytr0t1 moveto
               (to the background lightness ) showen
               (f\374r die UmgebungsHelligkeit ) showde
               (L*) jKs ifunc 0 eq {(80,) ibKs} if 
               (u) ibKs TBK
              } if %xchartl=00

xchartl 01 eq {%xchartl=01
               jlog 1 eq {(log \050) bKs} if
               (D) sMs 20 0 rmoveto
               (Y) jKs (/) bKs (D) sMs
               (Y) jKs (u) ibKs
               jlog 1 eq {(\051 ) bKs} if

               1800 xtr0 sub ytr0t0 moveto
               TBK ifunc 0 eq {(HAULAB) show} if
                   ifunc 1 eq {(CIELAB) show} if
                   ifunc 2 eq {(IECsRGB) show} if
                   ifunc 3 eq {(TUBsRGB) show} if
               TBK ( tristimulus value difference) showen
               (\255Normfarbwertdifferenz) showde

               1800 xtr0 sub ytr0t1 moveto
               (D) sMs 20rm (Y) jKs TBK
               ( normalized to ) showen
               ( normiert f\374r ) showde
               (D) sMs 20rm (Y) jKs (u) ibKs TBK

              } if %xchartl=01

xchartl 02 eq {%xchartl=02
               jlog 1 eq {(log [) bKs} if
               20 0 rmoveto (\050) bKs
               (D) sMs (Y/Y) jKs (\051 / \050) bKs
               (D) sMs (Y/Y) jKs (\051) bKs (u) ibKs 
               jlog 1 eq {(]) bKs} if               

               2800 xtr0 sub ytr0t0 moveto
               TBK ifunc 0 eq {(HAULAB) show} if
                   ifunc 1 eq {(CIELAB) show} if
                   ifunc 2 eq {(IECsRGB) show} if
                   ifunc 3 eq {(TUBsRGB) show} if
               (\255Y) jKs TBK
               ( sensitivity) showen
               (\255Empfindlichkeit) showde

               2800 xtr0 sub ytr0t1 moveto
               (normalized to ) showen
               (normiert f\374r ) showde
               (\050) bKs (D) sMs (Y/Y) jKs (\051) bKs (u) ibKs
              } if %xchartl=02


xchartl 03 eq {%xchartl=03
               jlog 1 eq {(log[) bKs} if
               20 0 rmoveto (\050) bKs
               (Y/) jKs (D) sMs (Y) jKs (\051 / \050) show
               (Y/) bKs (D) sMs (Y) jKs (\051) bKs (u) ibKs
               jlog 1 eq {(\051]) bKs} if

               2800 xtr0 sub ytr0t0 moveto
               TBK ifunc 0 eq {(HAULAB) show} if
                   ifunc 1 eq {(CIELAB) show} if
                   ifunc 2 eq {(IECsRGB) show} if
                   ifunc 3 eq {(TUBsRGB) show} if
               (\255Y) jKs TBK
               ( contrast) showen
               (\255Kontrast) showde

               2800 xtr0 sub ytr0t1 moveto
               (normalized to ) showen
               (normiert f\374r ) showde
               (\050) bKs (Y/) jKs (D) sMs (Y) jKs (\051) bKs (u) ibKs

              } if %xchartl=02
tfn %tfw
} bind def %proc_toptext

%***************************************************
/prog_L*top {%BEG prog_L*top ifunc=0,1,2,3
/20rm {20 0 rmoveto} def

%top-eq. BEG HAULAB, CIELAB, IECsRGB, TUBsRGB_proc_C02_C08
%for ifunc=0 to 3

ifunc 0 eq {/s1 134.6 def /n1 0.31 def  /d1 34.6 def} if
ifunc 1 eq {/s1 116 def /n1 1 3     div def /d1 16 def} if
ifunc 2 eq {/s1 100 def /n1 1 2.4   div def /d1  0 def} if
ifunc 3 eq {/s1 100 def /n1 1 10 ln div def /d1  0 def} if
/r1 s1 0.18 n1 exp mul def
/g1 r1 r1 d1 sub div def
/h1 d1 r1 d1 sub div def

x00t y00t moveto
(L*) jLs 20rm
(=) bLs 20rm (s) show 20rm
(\050Y/Y) jLs (n) ibLs (\051) bLs
(n) ebLs
(-) bLs 20rm (d) bLs

x01t y00t moveto
(\050) nLs (Y) kLs (n) iLs (=100,) nLs 20rm
           (Y) kLs (u) iLs (=) nLs TL Yu cvishow (,) nLs 20rm
TL (s=) show s1 cvsshow1x
ifunc 0 eq {(, n=0,31)     show} if
ifunc 1 eq {(, n=1/3)      show} if
ifunc 2 eq {(, n=1/2,4)    show} if
ifunc 3 eq {(, n=1/ln(10)) show} if
TL (, d=) show d1 cvsshow1x (\051) nLs

x00e y00t moveto
TL ([1a]) show

x00t y00t 250 sub moveto
(L*) jLs 20rm
(=) bLs 20rm (r \050) bLs
(Y/Y) jLs (u) ibLs (\051) bLs
(n) ebLs
(-) bLs 20rm (d) bLs

x01t y00t 250 sub moveto
(\050r = s ) nLs (\050Y) kLs (u) iLs
(/Y) kLs (n) iLs (\051) nLs (n) ebLs
(=) nLs TL r1 cvsshow2x (, ) nLs
(L*) kLs (u) iLs (= r-d =) nLs
TL r1 d1 sub cvsshow1x
(\051) nLs

x00e y00t 250 sub moveto
TL ([1b]) show

} bind def %END prog_L*top ifunc=1,2,3

%***************************************************
/prog_L*DL*u_C02 {%BEG prog_L*DL*u_C02
%for ifunc=0,1,2,3 HAULAB, CIELAB, IECsRGB, TUBsRGB
%ioutC02 0/1:no or with output of equations in main program

tfr %C02
/y00t 2990 def
/ydel 0270 def

prog_L*top %ifunc=0,1,2,3 ALL

%********************************************
ifunc 0 eq {%ifunc=0 special BEG HAULAB_C02
x00t y00t ydel 2 mul sub moveto
(L*/L*) jLs (u) ibLs 
TBL (=) show 20rm (g) show 20rm
(\050Y/Y) jLs (u) ibLs (\051) bLs
(n) ebLs
TBL (-h) show

x01t 450 add y00t ydel 2 mul sub moveto
TL (\050g=r/(r-d)=) show r1 r1 d1 sub div cvsshow2x
TL (, h=d/(r-d)=) show d1 r1 d1 sub div cvsshow2x
TL (\051) show

x00e y00t ydel 2 mul sub moveto
TL ([1c]) show

jlog 1 eq {%jlog=1
x00t y00t ydel 3 mul sub moveto
(log [\050) bLs
(L*/L*) jLs (u) ibLs
20rm (+) bLs 20rm (h) bLs 20rm
(\051 / g ] = n log \050) bLs
(Y/Y) jLs (u) ibLs (\051) bLs
( = ) bLs n1 cvsshow2x 20rm (log) bLs
(\050) show (Y/) jLs Yu cvishow (\051) bLs

x00e y00t ydel 3 mul sub moveto
TL ([1d]) show

x00t y00t ydel 3.9 mul sub moveto
(ln [\050) bLs
(L*/L*) jLs (u) ibLs
( + h) bLs
(\051 / g] = n) bLs 20rm (ln(10)) bLs 20rm (log\050) bLs 20rm
(Y/Y) jLs (u) ibLs (\051) bLs
( = ) bLs n1 10 ln mul cvsshow2x 20rm (log \050) bLs
(Y/) jLs Yu cvishow (\051) bLs

x00e y00t ydel 4.0 mul sub moveto
TL ([1e]) show

x00t y00t ydel 5 mul sub moveto
(\050) bLs (L*/L*) jLs (u) ibLs ( + h) bLs
(\051 / g ]) bLs
( = e) bLs 0 90 rmoveto
(n ln(10) log \050) bSs (Y/Y) jSs (u) ibSs (\051) bSs
0 -90 rmoveto
( = e) bLs 0 90 rmoveto
TBS n1 10 ln mul cvsshow2x 20rm %80 smaller
(log \050) bSs (Y/) jSs TBS Yu cvishow (\051) bSs
0 -90 rmoveto

x00e y00t ydel 5 mul sub moveto
TL ([1f]) show
} if %jlog=1
} if %ifunc=0 special END HAULAB_C02

%********************************************
ifunc 1 eq {%ifunc=1 special BEG CIELAB_C02
x00t y00t ydel 2 mul sub moveto
(L*/L*) jLs (u) ibLs 
TBL (=) show 20rm (g) show 20rm
(\050Y/Y) jLs (u) ibLs (\051) bLs
(n) ebLs
TBL (-h) show

x01t 450 add y00t ydel 2 mul sub moveto
TL (\050g=r/(r-d)=) show r1 r1 d1 sub div cvsshow2x
TL (, h=d/(r-d)=) show d1 r1 d1 sub div cvsshow2x
TL (\051) show

x00e y00t ydel 2 mul sub moveto
TL ([1c]) show


jlog 1 eq {%jlog=1
x00t y00t ydel 3 mul sub moveto
TBL (log [\050) show
(L*/L*) jLs (u) ibLs
TBL 20rm (+) show 20rm (h) show 20rm
TBL (\051 / g ] = n log \050) show
(Y/Y) jLs (u) ibLs (\051) bLs
x00e y00t ydel 3 mul sub moveto
TL ([1d]) show

x00t y00t ydel 3.9 mul sub moveto
TBL (ln [\050) show
(L*/L*) jLs (u) ibLs
TBL ( + h) show
TBL (\051 / g ] = ln(10) n log \050) show
(Y/Y) jLs (u) ibLs (\051) bLs
x00e y00t ydel 4.0 mul sub moveto
TL ([1e]) show

x00t y00t ydel 5 mul sub moveto
TBL (\050) show (L*/L*) jLs (u) ibLs
TBL ( + h) show
TBL (\051 / g ]) show
TBL ( = e) show (ln(10) n log \050) ebLs (Y/Y) jbLs (u) ebSs
(\051) ebLs
x00e y00t ydel 5 mul sub moveto
TL ([1f]) show
} if %jlog=1
} if %ifunc=1 special END CIELAB_C02

%********************************************
ifunc 2 eq {%ifunc=2 special BEG IECsRGB_C02
x00t y00t ydel 2 mul sub moveto
(L*/L*) jLs (u) ibLs 
TBL (=) show
(\050Y/Y) jLs (u) ibLs (\051) bLs
(n) ebLs
x00e y00t ydel 2 mul sub moveto
TL ([1c]) show

jlog 1 eq {%jlog=1
x00t y00t ydel 3 mul sub moveto
TBL (log \050) show
(L*/L*) jLs (u) ibLs
TBL (\051 = n log \050) show
(Y/Y) jLs (u) ibLs (\051) bLs
x00e y00t ydel 3 mul sub moveto
TL ([1d]) show

x00t y00t ydel 3.9 mul sub moveto
TBL (ln \050) show
(L*/L*) jLs (u) ibLs
TBL (\051 = ln(10) n log \050) show
(Y/Y) jLs (u) ibLs (\051) bLs
x00e y00t ydel 4.0 mul sub moveto
TL ([1e]) show

x00t y00t ydel 5 mul sub moveto
(L*/L*) jLs (u) ibLs
TBL ( = e) show (ln(10) n log \050) ebLs (Y/Y) jbLs (u) ebSs
(\051) ebLs
x00e y00t ydel 5 mul sub moveto
TL ([1f]) show
} if %jlog=1
} if %ifunc=2 special END IECsRGB_C02

%********************************************
ifunc 3 eq {%ifunc=3 special BEG TUBsRGB_C02
x00t y00t ydel 2 mul sub moveto
(L*/L*) jLs (u) ibLs 
TBL (=) show
(\050Y/Y) jLs (u) ibLs (\051) bLs
(1/ln(10)) ebLs
TBL ( \050ln(x)=ln(10) log(x)\051) show
x00e y00t ydel 2 mul sub moveto
TL ([1c]) show

jlog 1 eq {%jlog=1
x00t y00t ydel 3 mul sub moveto
TBL (log\050) show
(L*/L*) jLs (u) ibLs
TBL (\051=(1/ln(10)) log\050) show
(Y/Y) jLs (u) ibLs (\051) bLs
x00e y00t ydel 3 mul sub moveto
TL ([1d]) show

x00t y00t ydel 3.9 mul sub moveto
TBL (ln\050) show
(L*/L*) jLs (u) ibLs
TBL (\051=log\050) show 
(Y/Y) jLs (u) ibLs (\051) bLs
x00e y00t ydel 3.9 mul sub moveto
TL ([1e]) show

x00t y00t ydel 5 mul sub moveto
(L*/L*) jLs (u) ibLs
TBL (= e) show (log\050) ebLs (Y/Y) jbLs (u) ebSs
(\051) ebLs
x00e y00t ydel 5 mul sub moveto
TL ([1f]) show
} if %jlog=1
} if %ifunc=3 special END TUBsRGB_C02

} bind def %END prog_L*DL*u_C02

%***************************************************
/proc_YDYu_C04 {%BEG proc_YDYu_C04
%for ifunc=1,2,3 CIELAB, IECsRGB, TUBsRGB
%ioutC04 0/1:no or with output of equations in main program

tfr %C04
/y00t 2990 def
/ydel 0270 def

prog_L*top %for ifunc=1,2,3

ifunc 0 eq {%ifunc=0 BEG special HAULAB_C04
%2
x00t y00t ydel 2 mul sub moveto
(dY) jLs TBL ( = [) show
(Y) jLs (n) ibLs (/ \050 n s \051] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
x00e y00t ydel 2 mul sub moveto
TL ([2c]) show

%3
x00t y00t ydel 3 mul sub moveto
(dY) jLs (u) ibLs
TBL ( = [) show
(Y) jLs (n) ibLs TBL (/ \050 n s \051] ) show
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
TBL ( = ) show
100 n1 s1 mul div 18 100 div n1 exp mul cvsshow4x
x00e y00t ydel 3 mul sub moveto
TL ([2d]) show

%4
x00t y00t ydel 4 mul sub moveto
(dY / dY) jLs (u) ibLs 
TBL ( = ) show
(\050Y / Y) jLs (u) ibLs
(\051) bLs (1-n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([2e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
TBL (log\050) show (dY / dY) jLs (u) ibLs
TBL (\051 = (1-n) log) show
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([2f]) show
} if %jlog=1
} if %ifunc=0 special END HAULAB_C04

ifunc 1 eq {%ifunc=1 special CIELAB_C04
%2
x00t y00t ydel 2 mul sub moveto
(dY) jLs TBL ( = [) show
(Y) jLs (n) ibLs (/ \050 n s \051] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
x00e y00t ydel 2 mul sub moveto
TL ([2c]) show

%3
x00t y00t ydel 3 mul sub moveto
(dY) jLs (u) ibLs
TBL ( = [) show
(Y) jLs (n) ibLs TBL (/ \050 n s \051] ) show
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
TBL ( = ) show
100 n1 s1 mul div 18 100 div n1 exp mul cvsshow4x
x00e y00t ydel 3 mul sub moveto
TL ([2d]) show

%4
x00t y00t ydel 4 mul sub moveto
(dY / dY) jLs (u) ibLs 
TBL ( = ) show
(\050Y / Y) jLs (u) ibLs
(\051) bLs (1-n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([2e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
TBL (log\050) show (dY / dY) jLs (u) ibLs
TBL (\051 = (1-n) log) show
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([2f]) show
} if %jlog=1
} if %ifunc=1 special END CIELAB_C04

ifunc 2 eq 
ifunc 3 eq or {%ifunc=2,3 special IECsRGB/TUBsRGB_C04
%2
x00t y00t ydel 2 mul sub moveto
(dY) jLs TBL ( = [) show
(Y) jLs (n) ibLs (/ \050 n s \051] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
x00e y00t ydel 2 mul sub moveto
TL ([2c]) show

%3
x00t y00t ydel 3 mul sub moveto
(dY) jLs (u) ibLs
TBL ( = [) show
(Y) jLs (n) ibLs TBL (/ \050 n s \051] ) show
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
TBL ( = ) show
100 n1 s1 mul div 18 100 div n1 exp mul cvsshow4x
x00e y00t ydel 3 mul sub moveto
TL ([2d]) show

%4
x00t y00t ydel 4 mul sub moveto
(dY / dY) jLs (u) ibLs 
TBL ( = ) show
(\050Y / Y) jLs (u) ibLs
(\051) bLs (1-n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([2e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
TBL (log\050) show (dY / dY) jLs (u) ibLs
TBL (\051 = (1-n) log) show
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([2f]) show
} if %jlog=1
} if %ifunc=2,3 END IECsRGB/TUBsRGB_C04

} bind def %END proc_YDYu_C04

%***************************************************
/proc_dYDY_C06 {%BEG proc_dYDY_C06
%for ifunc=1,2,3 CIELAB, IECsRGB, TUBsRGB
%ioutC06 0/1:no or with output of equations in main program

tfr %C06

/y00t 2990 def
/ydel 0270 def

prog_L*top %for ifunc=0,1,2,3

ifunc 0 eq {%ifunc=0 special BEG HAULAB_C06
%2
x00t y00t ydel 2 mul sub moveto
(dY / Y) jLs TBL ( = [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
( / ) bLs (Y) jLs
x00e y00t ydel 2 mul sub moveto
TL ([3c]) show

%3
x00t y00t ydel 3 mul sub moveto
((dY / Y)) jLs (u) ibLs
( = [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
( / ) bLs (Y) jLs (u) ibLs
x00e y00t ydel 3 mul sub moveto
TL ([3d]) show

%4
x00t y00t ydel 4 mul sub moveto
((dY / Y)) jLs ( / ) bLs ((dY / Y)) jLs (u) ibLs
( = ) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs (-n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([3e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
(log [) bLs
((dY / Y)) jLs ( / ) bLs ((dY / Y)) jLs (u) ibLs
(] = (-n) log) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([3f]) show
} if %jlog=1
} if %ifunc=0 special END HAULAB_C06

ifunc 1 eq {%ifunc=1 special BEG CIELAB_C06
%2
x00t y00t ydel 2 mul sub moveto
(dY / Y) jLs TBL ( = [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
( / ) bLs (Y) jLs
x00e y00t ydel 2 mul sub moveto
TL ([3c]) show

%3
x00t y00t ydel 3 mul sub moveto
((dY / Y)) jLs (u) ibLs
( = [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
( / ) bLs (Y) jLs (u) ibLs
x00e y00t ydel 3 mul sub moveto
TL ([3d]) show

%4
x00t y00t ydel 4 mul sub moveto
((dY / Y)) jLs ( / ) bLs ((dY / Y)) jLs (u) ibLs
( = ) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs (-n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([3e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
(log [) bLs
((dY / Y)) jLs ( / ) bLs ((dY / Y)) jLs (u) ibLs
(] = (-n) log) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([3f]) show
} if %jlog=1
} if %ifunc=1 special END CIELAB_C06

ifunc 2 eq 
ifunc 3 eq or {%ifunc=2,3 special BEG IECsRGB/TUBsRGB_C06
%2
x00t y00t ydel 2 mul sub moveto
(dY / Y) jLs TBL ( = [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
( / ) bLs (Y) jLs
x00e y00t ydel 2 mul sub moveto
TL ([3c]) show

%3
x00t y00t ydel 3 mul sub moveto
((dY / Y)) jLs (u) ibLs
( = [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs
( / ) bLs (Y) jLs (u) ibLs
x00e y00t ydel 3 mul sub moveto
TL ([3d]) show

%4
x00t y00t ydel 4 mul sub moveto
((dY / Y)) jLs ( / ) bLs ((dY / Y)) jLs (u) ibLs
( = ) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs (-n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([3e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
(log [) bLs
((dY / Y)) jLs ( / ) bLs ((dY / Y)) jLs (u) ibLs
(] = (-n) log) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([3f]) show
} if %jlog=1
} if %ifunc=2,3 special END IECsRGB/TUBsRGB_C06

} bind def %END proc_dYDY_C06

%***************************************************
/proc_YDdY_C08 {%BEG proc_YDdY_C08
%for ifunc=1,2,3 CIELAB, IECsRGB, TUBsRGB
%ioutC08 0/1:no or with output of equations in main program

tfr %C08

/y00t 2990 def
/ydel 0270 def

prog_L*top %for ifunc=1,2,3

ifunc 0 eq {%ifunc=0 special BEG HAULAB_C08
%2
x00t y00t ydel 2 mul sub moveto
(Y / dY) jLs ( = ) bLs (Y) jLs ( / { [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs ( }) bLs
x00e y00t ydel 2 mul sub moveto
TL ([4c]) show

%3
x00t y00t ydel 3 mul sub moveto
((Y / Y)) jLs (u) ibLs
( = ) bLs (Y) jLs (u) ibLs ( / { [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs ( }) bLs
x00e y00t ydel 3 mul sub moveto
TL ([4d]) show

%4
x00t y00t ydel 4 mul sub moveto
((Y / dY)) jLs ( / ) bLs ((Y / dY)) jLs (u) ibLs
( = ) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs (n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([4e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
(log [) bLs
((Y / dY)) jLs ( / ) bLs ((Y / dY)) jLs (u) ibLs
(] = (n) log) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([4f]) show
} if %jlog=1
} if %ifunc=0 special END HAULAB_C08

ifunc 1 eq {%ifunc=1 special BEG CIELAB_C08
%2
x00t y00t ydel 2 mul sub moveto
(Y / dY) jLs ( = ) bLs (Y) jLs ( / { [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs ( }) bLs
x00e y00t ydel 2 mul sub moveto
TL ([4c]) show

%3
x00t y00t ydel 3 mul sub moveto
((Y / Y)) jLs (u) ibLs
( = ) bLs (Y) jLs (u) ibLs ( / { [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs ( }) bLs
x00e y00t ydel 3 mul sub moveto
TL ([4d]) show

%4
x00t y00t ydel 4 mul sub moveto
((Y / dY)) jLs ( / ) bLs ((Y / dY)) jLs (u) ibLs
( = ) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs (n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([4e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
(log [) bLs
((Y / dY)) jLs ( / ) bLs ((Y / dY)) jLs (u) ibLs
(] = (n) log) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([4f]) show
} if %jlog=1
} if %ifunc=1 special END CIELAB_C08

ifunc 2 eq 
ifunc 3 eq or {%ifunc=2,3 special BEG IECsRGB/TUBsRGB_C08
%2
x00t y00t ydel 2 mul sub moveto
(Y / dY) jLs ( = ) bLs (Y) jLs ( / { [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs ( }) bLs
x00e y00t ydel 2 mul sub moveto
TL ([4c]) show

%3
x00t y00t ydel 3 mul sub moveto
((Y / Y)) jLs (u) ibLs
( = ) bLs (Y) jLs (u) ibLs ( / { [ ) bLs
(\050) bLs (Y) jLs (n) ibLs (/ \050 n s \051 ] ) bLs
(\050Y) jLs (u) ibLs ( / Y) jLs (n) ibLs
(\051) bLs (1-n) ebLs ( }) bLs
x00e y00t ydel 3 mul sub moveto
TL ([4d]) show

%4
x00t y00t ydel 4 mul sub moveto
((Y / dY)) jLs ( / ) bLs ((Y / dY)) jLs (u) ibLs
( = ) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs (n) ebLs
x00e y00t ydel 4 mul sub moveto
TL ([4e]) show

%5
jlog 1 eq {%jlog=1
x00t y00t ydel 5 mul sub moveto
(log [) bLs
((Y / dY)) jLs ( / ) bLs ((Y / dY)) jLs (u) ibLs
(] = (n) log) bLs
(\050Y / Y) jLs (u) ibLs
(\051) bLs
x00e y00t ydel 5 mul sub moveto
TL ([4f]) show
} if %jlog=1

} if %ifunc=2,3 special END IECsRGB/TUBsRGB_C08

} bind def %END proc_YDdY_C08

%***************************************************
/proc_mdu {%BEG proc_mdu
%for ifunc=0,1,2,3 HAULAB, CIELAB, IECsRGB, TUBsRGB
%for C02, C04, C06, C08

xtr neg ytr neg translate %new cero point

%available Yx0i and Yxui from Yi_Yxyi
%for x: xchartl=0_C02, xchartl=1_C04, xchartl=2_C06, xchartl=3_C08

/n090 090 def
/n004 004 def
/t090 (90) def
/t004 (4) def
/x00x x00t 200 sub def

        x00x 1 MULY mul 120 add moveto
        jlog 0 eq {%jlog=0
                   (m) jLs (u) ibLs t090 ibLs (_) ibLs t004 ibLs TBL ( = ) show
                   Yx0i n090 get Yx0i n004 get sub
                   Yi   n090 get Yi   n004 get sub div cvsshow3x (, ) show
                   (f) jLs t090 ibLs (=) bLs TBL Yx0i n090 get cvishow (, ) bLs
                   (f) jLs t004 ibLs (=) bLs TBL Yx0i n004 get cvishow
                  } %jlog=0
                  { %jlog=1
                   (m) jLs (nu) ibLs TBL 
                   xchartl 0 eq {( = n = ) show n cvsshow3x} if
                   xchartl 1 eq {( = 1-n = ) show 1 n sub cvsshow3x} if
                   xchartl 2 eq {( = -n = ) show n neg cvsshow3x} if
                   xchartl 3 eq {( = n = ) show n cvsshow3x} if 
                  } ifelse %jlog=0,1
        x00x 1 MULY mul 200 sub moveto
        (m) jLs (u) ibLs TBL ( = ) show
        Yx0i 20 get jlog 1 eq {log} if
        Yx0i iu get jlog 1 eq {log} if sub
        Yi 20 get jlog 1 eq {log} if
        Yi iu get jlog 1 eq {log} if sub div cvsshow3x

xtr ytr translate %new cero point

} bind def %proc_mdu

%*********************************************************************
/proc_C02_ALOG_L*DL*u {%BEG proc_C02_ALOG_L*DTu*

jlog 0 eq {%jlog=0
/Fx0log -2.0 def
/Fy0lin  1.0 def
/xtr 2000 def
/ytr 0000 def
/iys 1000 def
/MULY iys 1 div def %scale=1
/ytrl ytr 1000 div def %lin shift
           } %jlog=0
           { %jlog=1
/Fx0log -2.0 def
/Fy0log -1.0 def
/xtr 2000 def
/ytr 1000 def
/ytrl ytr 1000 div def %log shift
           } ifelse %jlog=1

xtr ytr translate %new cero point
%0 not available

Y_curve

[ ] 0 setdash
50 setlinewidth %50 setlinewidth
0 setgray

%C02 xchartl=0
proc_mdu

/i18 18 def
tfz %tfz
[100] 0 setdash
jlog 0 eq {%jlog=0,1
Fx0log MULX mul L*i i18 get L*u div MULY mul moveto
i18 log MULX mul L*i i18 get L*u div MULY mul lineto
iu log MULX mul 0                  MULY mul lineto stroke
          } %jlog=0
          { %jlog=1
Fx0log MULX mul L*i i18 get L*u div log MULY mul moveto
i18 log MULX mul L*i i18 get L*u div log MULY mul lineto
i18 log MULX mul Fy0log                 MULY mul lineto stroke
          } ifelse %jlog=0,1
[ ] 0 setdash
tfn

ioute 1 eq {cero_line} if

xtr neg ytr neg translate %old cero point

tfn %tfw
/Y001DYu L*i  i1 get 
         L*i  i18 get div def
/Y100DYu L*i 100 get 
         L*i  i18 get div def
tfg
[100] 0 setdash
0 Y001DYu jlog 1 eq {log} if ytrl add MULY mul moveto
0 Y100DYu jlog 1 eq {log} if ytrl add MULY mul lineto stroke
[ ] 0 setdash
tfn %tfw

/y0 2900 def
/x1  800 def
/y1 2700 def
/xdel 800 def
/ydel 300 def

%for ifunc=0,1,2,3 HAULAB, CIELAB, IECsRGB, TUBsRGB
/ioutC02 1 def %0/1:no or with output of equations in main program
ioutC02 1 eq {prog_L*DL*u_C02} if

/y00a 0850 def
proc_appli

} def %END proc_C02_ALOG_L*DL*u

%*********************************************************************
/proc_C04_ALOG_DLn {%BEG proc_C04_ALOG_DLn
jlog 0 eq {%jlog=0
/Fx0log -2.0 def
/Fy0lin  1.0 def
/xtr 2000 def
/ytr 0000 def
/iys 1000 def
/MULY iys 2.0 div def %scale=2.0
/ytrl ytr 1000 div def %lin shift
           } %jlog=0
           { %jlog=1
%/Fx0log -2.0 def
%/Fy0log -1.0 def
/xtr 2000 def
/ytr 1000 def
/ytrl ytr 1000 div def %log shift
           } ifelse %jlog=1

xtr ytr translate %new cero point
%0 not available

Y_curve

%C04 xchartl=1
proc_mdu

tfz %tfz
[100] 0 setdash
jlog 0 eq {%jlog=0,1
/rYdY 1 def
-2.    MULX mul rYdY MULY mul moveto
iu log MULX mul rYdY MULY mul lineto
iu log MULX mul 0    MULY mul lineto stroke
          } %jlog=0
          { %jlog=1
/rYdY 1. def
-2.    MULX mul rYdY log MULY mul moveto
iu log MULX mul rYdY log MULY mul lineto
iu log MULX mul -1.      MULY mul lineto stroke
          } ifelse %jlog=0,1
[ ] 0 setdash
tfn

xtr neg ytr neg translate %new cero point

tfn %tfw
/Y001DYu dYi  i1 get dYu div def
/Y100DYu dYi 100 get dYu div def
tfg
[100] 0 setdash
0 Y001DYu jlog 1 eq {log} if ytrl add MULY mul moveto
0 Y100DYu jlog 1 eq {log} if ytrl add MULY mul lineto stroke
[ ] 0 setdash

tfn %tfw
/y0 2400 y0del sub def
/ydel 300 def

%C04
/x02t 0900 def
/x03t 2300 def
/x04t 3300 def

%for ifunc=0,1,2,3 HAULAB, CIELAB, IECsRGB, TUBsRGB
/ioutC04 1 def %0/1:no or with output of equations in main program
ioutC04 1 eq {proc_YDYu_C04} if

/y00a 0750 def
proc_appli

} def %END proc_C04_ALOG_DLn

%*********************************************************************
/proc_C06_ALOG_DL_Ln {%BEG proc_C06_ALOG_DL_Ln
jlog 0 eq {%jlog=0
/Fx0log -2.0 def
/Fy0lin  1.0 def
/xtr 2000 def
/ytr 0000 def
/iys 1000 def
/MULY iys 2.0 div def %scale=2.0
/ytrl ytr 1000 div def %lin shift
           } %jlog=0
           { %jlog=1
%/Fx0log -2.0 def
%/Fy0lin  1.0 def
/xtr 2000 def
/ytr 1000 def
/ytrl ytr 1000 div def %log shift
           } ifelse %jlog=1

xtr ytr translate %new cero point
%0 not available

Y_curve

%C06 xchartl=2
proc_mdu


tfz %tfz
[100] 0 setdash
jlog 0 eq {%jlog=0,1
[100] 0 setdash
/rYdY 1 def
-2.    MULX mul rYdY MULY mul moveto
iu log MULX mul rYdY MULY mul lineto
iu log MULX mul 0    MULY mul lineto stroke
          } %jlog=0
          { %jlog=1
/rYdY 1 def
-2.    MULX mul rYdY log MULY mul moveto
iu log MULX mul rYdY log MULY mul lineto
iu log MULX mul -1.      MULY mul lineto stroke
          } ifelse %jlog=0,1
[ ] 0 setdash
tfn

%/iu 18 def
%/rYdY Yi i18 get dYi i18 get div
%      Yu        dYu        div div e10D30 exp def
%-2.    MULX mul rYdY log MULY mul moveto
%iu log MULX mul rYdY log MULY mul lineto
%iu log MULX mul -1.      MULY mul lineto stroke

[ ] 0 setdash
tfn

ioute 1 eq {cero_line} if

xtr neg ytr neg translate %old cero point

tfn %tfw
/Y001DYu dYi  i1 get Yi  i1 get div 
         dYu         Yu         div div def
/Y100DYu dYi 100 get Yi 100 get div 
         dYu         Yu         div div def
tfg
[100] 0 setdash
0 Y001DYu jlog 1 eq {log} if ytrl add MULY mul moveto
0 Y100DYu jlog 1 eq {log} if ytrl add MULY mul lineto stroke
[ ] 0 setdash
tfn %tfw

/y0 2400 y0del sub def
/ydel 300 def

%C06
/x02t 0900 def
/x03t 2300 def
/x04t 3300 def

%for ifunc=0,1,2,3 HAULAB, CIELAB, IECsRGB, TUBsRGB
/ioutC06 1 def %0/1:no or with output of equations in main program

ioutC06 1 eq {proc_dYDY_C06} if

/y00a 01400 def
proc_appli

} def %END proc_C06_ALOG_DL_Ln

%*********************************************************************
/proc_C08_ALOG_L_DLn {%BEG proc_C08_ALOG_L_DLn
jlog 0 eq {%jlog=0
/Fx0log -2.0 def
/Fy0lin  1.0 def
/xtr 2000 def
/ytr 0000 def
/iys 1000 def
/MULY iys 1 div def %scale=1
/ytrl ytr 1000 div def %lin shift
           } %jlog=0
           { %jlog=1
%/Fx0log -2.0 def
%/Fy0log -1.0 def
/xtr 2000 def
/ytr 1000 def
/ytrl ytr 1000 div def %log shift
           } ifelse %jlog=1

xtr ytr translate %new cero point

Y_curve

%C08 xchartl=3
proc_mdu

tfz %tfz
[100] 0 setdash
jlog 0 eq {%jlog=0,1
[100] 0 setdash
/rYdY 1 def
-2.    MULX mul rYdY MULY mul moveto
iu log MULX mul rYdY MULY mul lineto
iu log MULX mul 0    MULY mul lineto stroke
          } %jlog=0
          { %jlog=1
/rYdY 1 def
-2.    MULX mul rYdY log MULY mul moveto
iu log MULX mul rYdY log MULY mul lineto
iu log MULX mul -1.      MULY mul lineto stroke
          } ifelse %jlog=0,1
[ ] 0 setdash
tfn

ioute 1 eq {cero_line} if

xtr neg ytr neg translate %old cero point

tfn %tfw
/Y001DYu Yi  i1 get dYi  i1 get div 
         Yu         dYu         div div def
/Y100DYu Yi 100 get dYi 100 get div 
         Yu         dYu         div div def
tfg
[100] 0 setdash
0 Y001DYu jlog 1 eq {log} if ytrl add MULY mul moveto
0 Y100DYu jlog 1 eq {log} if ytrl add MULY mul lineto stroke
[ ] 0 setdash

%C08
/y0 2400 y0del sub def
/ydel 300 def

%for ifunc=0,1,2,3 HAULAB, CIELAB, IECsRGB, TUBsRGB
/ioutC08 1 def %0/1:no or with output of equations in main program
ioutC08 1 eq {proc_YDdY_C08} if

/y00a 0850 def
proc_appli

} def %END proc_C08_ALOG_L_DL
%%EndProlog
gsave

/lanind 1 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 3 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 /Fxx (d) def} if %colorm=0, deintp=0
colorm 0 eq deintp 1 eq and {/Txx (e) def /Fxx (e) def} if %colorm=0, deintp=1
colorm 1 eq deintp 0 eq and {/Txx (dd) def /Fxx (d) def} if %colorm=1, deintp=0
colorm 1 eq deintp 1 eq and {/Txx (de) def /Fxx (e) def} if %colorm=1, deintp=1
xchart 0 eq {/Txx (-) def /Fxx (-) def} if %always independent of intended output

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

gsave
%XCHA01.PS END

%STOPB
%ifunc, i1, i2 defined near line 292

ifunc 0 eq {Yi_L*i_dYi_H_0} if
ifunc 1 eq {Yi_L*i_dYi_C_0} if
ifunc 2 eq {Yi_L*i_dYi_I_0} if
ifunc 3 eq {Yi_L*i_dYi_T_0} if

/cvishow0 {cvi 6 string cvs show} def

/kchartl 0 def %0:left page, 1:right page
/pchartl 0 def %0:top  page, 4:down  page %not used
/jlog 1 def %0,1 without/with log

72 90 translate

0.010 MM dup scale

/xbtex0 1 def  %xbtex=0 for files Y1(0/1)-(3/7)n.EPS
xbtex0 1 eq {%xbtex0=1

40 setlinewidth
/ymax1 08550 def
/xmax1 12250 def
1.0 setgray
0 0  moveto xmax1 0 rlineto 0 ymax1 rlineto xmax1 neg 0 rlineto
closepath fill
0 setgray
0 0  moveto xmax1 0 rlineto 0 ymax1 rlineto xmax1 neg 0 rlineto
closepath stroke

TK
0 setgray
150  /Times-ISOL1 FS
150 -140 moveto
(hep80-3n) show
} if %xbtex0=1

/xpos [00100 06150 00100 06150] def
/ypos [04480 04480 00220 00220] def

/xchartl 0 def
0 1 3 {/xchartl exch def %xchartl=0,3

gsave

Yi_Yxyi

%STOPC

20 setlinewidth

xpos xchartl get ypos xchartl get translate

0 setgray
150  /Times-ISOL1 FS
150 -140 moveto
(hep8) show kchartl cvishow0 (-) show
xchartl 1 add pchartl add cvishow0
(a) show %a

/xwidth 6000 def
/ywidth 4000 def
25 setlinewidth
1 1 1 setrgbcolor
0 0 moveto xwidth     0 rlineto 0 ywidth rlineto
           xwidth neg 0 rlineto closepath fill

0 setlinewidth
0 0 moveto xwidth     0 rlineto 0 ywidth rlineto
           xwidth neg 0 rlineto closepath clip

20 setlinewidth
0 setgray
0 0 moveto xwidth     0 rlineto 0 ywidth rlineto
           xwidth neg 0 rlineto closepath stroke


tfn %tfw
/xtr0 380 def
/ytr0 280 def

xtr0 ytr0 translate

%*********************************************************************
50 setlinewidth %50 setlinewidth
tfn
 0 0 moveto 5000 0 rlineto stroke
 0 0 moveto 0 3100 rlineto stroke

5000 100 add 0 moveto
-100 50 rlineto 0 -100 rlineto closepath fill
0 3100 100 add moveto
-50 -100 rlineto 100 0 rlineto closepath fill

TBL
/tx [(-2) (-1) (  0) (  1) (  2)] def
/txl [( ) (0,1) (  1) ( 10) (100)] def

%!x-Achse: 100 Einheiten = 0600 Skalen-Einheiten
0 1 4 {/i exch def
       /ixt {-150 i 1000 mul add} def
       /ixl { 000 i 1000 mul add} def
        ixt -230 moveto tx i get exec show
        tfb
        i 1 ge {ixt 100 moveto txl i get exec show} if
        tfn
        ixl   60 moveto 0 -120 rlineto stroke
     } for

tfn %tfw
3300  100 moveto (Y) jLs (u) ibLs TBL (=18) show
tfn

/ixtt 4.5 1000 mul def
/iytt -200 def
ixtt 200 add iytt moveto
(log) bLs 20rm (Y) jLs
/iytt 050 def
ixtt 200 add iytt moveto
tfb (Y) jLs tfn

tfn %tfw
2000 0 moveto 1900 0 rlineto stroke
tfg
[100] 0 setdash
2000 0 moveto 1900 0 rlineto stroke
[ ] 0 setdash

jlog 0 eq {%jlog=0,1
tfn
%!y-Achse: 100 S-Einheiten = 1000 Skalen-Einheiten
/j1y 2 def %default
/j2y 4 def
%ifunc 0 eq {%ifunc=0 HAULAB
%xchartl  0 eq {/j1y 4 def /j2y 6 def} if
%xchartl  2 eq {/j1y 1 def /j2y 3 def} if
%xchartl  4 eq {/j1y 0 def /j2y 2 def} if
%xchartl  6 eq {/j1y 4 def /j2y 6 def} if
%           } if %ifunc=0 HAULAB
%ifunc 1 eq {%ifunc=0 CIELAB
%xchartl  0 eq {/j1y 3 def /j2y 5 def} if
%xchartl  2 eq {/j1y 2 def /j2y 4 def} if
%xchartl  4 eq {/j1y 1 def /j2y 3 def} if
%xchartl  6 eq {/j1y 3 def /j2y 5 def} if
%           } if %ifunc=1 CIELAB
TBL
%yshift -1000  0   1000    2000    3000    4000
%        0     1      2      3       4       5       6       7       8
/ty [(  -3)(  -2)(  -1) (    0) (    1) (    2) (    3) (    4) (    5)] def
/tyl[(0,001)(0,01) (0,1)     (1)    (10)   (100)  (1000) (10000) (100000)] def

ifunc 0 eq
ifunc 2 eq or 
ifunc 3 eq or {%ifunc=0,2,3 HAULAB, IECsRGB, TUBsRGB 
/ty0[(    0)(  500)(1000)(1500)] def
/ty1[(  0)(  1)(  2)(  3)] def
/ty2[(0,0)(0,2)(0,4)(0,6)] def
/ty3[(  0)(  2)(  4)(  6)] def
/ty4[(0,00)(0,01)(0,02)(0,03)] def
/ty5[(  0)(  2)(  4)(  6)] def
/ty6[(    0)(200)(400)(600)] def
/ty7[(  0)(  1)(  2)(  3)] def
      } if %ifunc=0,2,3 HAULAB, IECsRGB, TUBsRGB
ifunc 1 eq {%ifunc=1 CIELAB
/ty0[(    0)(  50)(100)(150)] def
/ty1[(  0)(  1)(  2)(  3)] def
/ty2[(  0)(  2)(  4)(  6)] def
/ty3[(  0)(  2)(  4)(  6)] def
/ty4[(0,0)(0,1)(0,2)(0,3)] def
/ty5[(  0)(  2)(  4)(  6)] def
/ty6[(  0) (20) (40)(60) ] def
/ty7[(  0)(  1)(  2)(  3)] def
      } if %ifunc=1 CIELAB

/j1y 0 def
/j2y 3 def
j1y 1 j2y {/j exch def
           /jyt {-50 j j1y sub 1000 mul add} def
           /jyl {000 j j1y sub 1000 mul add} def
           -400
           jyt moveto
           xchartl 0 eq {ty1 j get show} if
           xchartl 1 eq {ty3 j get show} if
           xchartl 2 eq {ty5 j get show} if
           xchartl 3 eq {ty7 j get show} if
           tfn
           -60 jyl moveto 120 0 rlineto stroke
          } for
        } %jlog=0
        { %jlog=1
tfn
%!y-Achse: 100 S-Einheiten = 1000 Skalen-Einheiten
/j1y 2 def %default
/j2y 5 def
%ifunc 0 eq {%ifunc=0 HAULAB
%xchartl  0 eq {/j1y 4 def /j2y 6 def} if
%xchartl  2 eq {/j1y 1 def /j2y 4 def} if
%xchartl  4 eq {/j1y 0 def /j2y 3 def} if
%xchartl  6 eq {/j1y 4 def /j2y 6 def} if
%           } if %ifunc=0 HAULAB
%ifunc 1 eq {%ifunc=0 CIELAB
%xchartl  0 eq {/j1y 3 def /j2y 5 def} if
%xchartl  2 eq {/j1y 2 def /j2y 5 def} if
%xchartl  4 eq {/j1y 1 def /j2y 4 def} if
%xchartl  6 eq {/j1y 3 def /j2y 5 def} if
%           } if %ifunc=1 CIELAB
TBL
%yshift -1000  0   1000    2000    3000    4000
%        0     1      2      3       4       5       6       7       8
/ty [(  -3)(  -2)(  -1) (    0) (    1) (    2) (    3) (    4) (    5)] def
/tyl[(0,001)(0,01) (0,1)     (1)    (10)   (100)  (1000) (10000) (100000)] def
j1y 1 j2y {/j exch def
           /jyt {-50 j j1y sub 1000 mul add} def
           /jyl {000 j j1y sub 1000 mul add} def
           -400 jyt moveto ty j get show
           tfb
           j j1y 1 add ge {100 jyt moveto tyl j get show} if
           tfn
           -60 jyl moveto 120 0 rlineto stroke
          } for
        } ifelse %jlog=0,1

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

%BEG C02, C04, C06, C08 ********************************************
%jlog 0:without log, 1:with log in main program

ifunc 0 eq {/n 0.3100 def} if %HAULAB
ifunc 1 eq {/n 1 3.0 div def} if
ifunc 2 eq {/n 1 2.4 div def} if
ifunc 3 eq {/n 1 10 ln div def} if

proc_toptext

/y0del 100 def
50 setlinewidth

%****************************************************************
tfb
/xtfb 00 def
xchartl 00 eq {100 3200 xtfb sub moveto
               (L*/L*) jLs ifunc 0 eq {(80,) ibLs} if (u) ibLs
              } if

xchartl 01 eq {100 3200 xtfb sub moveto
               (D) sKs (Y) jLs (/) bLs
               (D) sKs (Y) jLs (u) ibLs
              } if

xchartl 02 eq {100 3200 xtfb sub moveto
               (S) jLs (r)  ibLs (/) bLs (S) jLs (ru) ibLs
               (=\050) bLs (D) sKs (Y/Y) jLs (\051/\050) bLs
               (D) sKs (Y/Y) jLs (\051) bLs (u) ibLs
              } if

xchartl 03 eq {100 3200 xtfb sub moveto
               (C) jLs (r)  ibLs (/) bLs
               (C) jLs (ru) ibLs
               (=\050) bLs (Y/) jLs (D) sKs (Y) jLs (\051/\050) bLs
               (Y) jLs (/) bLs (D) sKs (Y) jLs (\051) bLs (u) ibLs
              } if
tfn

%**************************************************************
xchartl 00 eq {proc_C02_ALOG_L*DL*u} if
xchartl 01 eq {proc_C04_ALOG_DLn} if
xchartl 02 eq {proc_C06_ALOG_DL_Ln} if
xchartl 03 eq {proc_C08_ALOG_L_DLn} if

%END C01_C08**********************************************

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

xtr0 neg ytr0 neg translate
%**************************************************************

xpos xchartl get neg ypos xchartl get neg translate

grestore

} for %xchartl=0,3

showpage
grestore

%%Trailer