% FIGURE 5.8.1
      function OMEGA = OMEGA581(IDER,X)
      global TKP1
%                              EVALUATE OMEGA (IDER=0) OR ITS IDER-TH
%                              DERIVATIVE AT X
      OMEGA = exp(TKP1)*HERM('H',0,IDER,X);

