% FIGURE 5.4.1
      function OMEGA = OMEGA541(IDER,X)
%                              EVALUATE OMEGA (IDER=0) OR ITS IDER-TH
%                              DERIVATIVE AT X
      OMEGA = 0.0;