% FIGURE 5.4.2
      function OMEGA = OMEGA542(IDER,X)
      global N
%                              EVALUATE OMEGA (IDER=0) OR ITS IDER-TH
%                              DERIVATIVE AT X
      OMEGA = 1.1*HERM('S',N,IDER,X);
