Title ; FOOTNOTE ' ';
data anno;  f=1/sqrt(2*3.14159*.36);
  xsys='2'; ysys='2';
  length color function $ 8;

  type = 21; function = 'label'; color='orange';
  x =-1.5; y= 1.5; text='Lognormal Densities';
  angle=90; output;
  x = -2.7;  y=exp(2+.5*.36)/50;  text = '1';  output;
  x = -1.1;  y=exp(3+.5*.36)/50;  text = '2';  output;
  x = -0.4;  y=exp(4+.5*.36)/50 ; text = '3';  output;

  type = 21; function = 'label'; color='orange';
  x = 3; y=-.1; text='Normal Densities'; angle=0;    output;
  x = 2; y=.55; text = '1';     output;
  x = 3; y=.55; text = '2';     output;
  x = 4; y=.55; text = '3';     output;

            *  (first pair of densities)  ;
  x=-0.4; function = 'move'; color='blue';
  y=f*exp(-.5*((x-2)/.6)**2); type=.5; output;
  y=exp(x)/50; x = -probnorm((x-2)/.6); type=1.5; output;
 function='draw'; do  x=-.4 to 4.4 by .1;
  y=f*exp(-.5*((x-2)/.6)**2); type=1; output;
  if d=. then d=0; p = .1*y; y = exp(x)/50; int=y-d;
  d=y; x1 = -p/int; type=2;   output;
           end;

            *  (second pair of densities)  ;
  x=0.6; function = 'move'; color='cyan'; d=0;
  y=f*exp(-.5*((x-3)/.6)**2); type=2.5; output;
  y=exp(x)/50; x = -probnorm((x-3)/.6); type=3.5; output;
 function='draw'; do  x=0.6 to 5.4 by .1;
  y=f*exp(-.5*((x-3)/.6)**2); type=3; output;
  p = .1*y; y = exp(x)/50; int=y-d;
  d=y; x1 = -p/int; type=4;output;
  end;

            *  (third pair of densities)  ;
  x=1.6; function = 'move'; color='green'; d=0;
  y=f*exp(-.5*((x-4)/.6)**2); type=4.5; output;
  y=exp(x)/50; x = -probnorm((x-3)/.6); type=5.5; output;
 function='draw'; do  x=1.6 to 5.5 by .1;
  y=f*exp(-.5*((x-4)/.6)**2); type=5;  output;
  p = .1*y; y = exp(x)/50; int=y-d;
  d=y; x1 = -p/int; type=6;   output;

               *   put in the reference lines;

        if x>3.69 and x < 3.71 then do;  type=-1; put x type;
        function = 'move'; y=0;                  output;
        function = 'draw'; line=1;
        y=f*exp(-.5*((x-4)/.6)**2);              output;
            color='red'; y=exp(x)/50; line=2;    output;
            hold = x; x=0;                       output;
            color='green';         x=x1;line=1;  output; x=hold;
                 end;
        if x>3.59 and x < 3.61 then do; type=-1;  put x x1 hold;
        function = 'move'; y=0;  output;
        function = 'draw'; line=1;
        y=f*exp(-.5*((x-4)/.6)**2);              output;
            color='red'; y=exp(x)/50; line=2;    output;
            hold = x; x=0;                       output;
            color='green';         x=x1;line=1;  output; x=hold;
                 end;
        if x>4.29 and x < 4.31 then do; type=-1;  put x;
        function = 'move'; y=0;  output;
        function = 'draw'; line=1;
        y=f*exp(-.5*((x-4)/.6)**2);              output;
            color='red'; y=exp(x)/50; line=2;    output;
            hold = x; x=0;                       output;
            color='green';         x=x1;line=1;  output; x=hold;
                 end;
        if x>4.39 and x < 4.41 then do; type=-1;  put x;
        function = 'move'; y=0;  output;
        function = 'draw'; line=1;
        y=f*exp(-.5*((x-4)/.6)**2);              output;
            color='red'; y=exp(x)/50; line=2;    output;
            hold = x; x=0;                       output;
            color='green';         x=x1;line=1;  output; x=hold;
                 end;
        end;
proc sort; by type;
data anno; set anno; if y < 3;
if type=2 or type=4 or type=6 then x=x1;
drop x1 p int f;
data next; do x = 0 to 5.5 by .1;
   y = exp(x)/50; output; end;
   data next; set next; if y < 3;
proc gplot anno=anno;
   plot y*x/ vref=0 href=0 chref=black cvref=black haxis=-6.5 to 6
   vaxis = -.2 to 3 noaxes nolegend;
  symbol1 v=diamond i=none c=black  h=1;
run;