function testOPO()

  opo = Data( 'OPO_Daten.xlsx' );
  opo.getData();
  
  plot( opo );
  title( 'Laser-Spektrum' )
  xlabel( 'Lambda/my' )
  ylabel( 'Emission' )
  
  % Integral
  f = trapez( opo )
  % Mittelwert 
  mw = mean( opo )
  % Standardabweichung 
  sa = std( opo )
  
  x = linspace( min(opo.x), max(opo.x), 200 );
  
  y = f * exp( -(x-mw).^2 / (2*sa^2) ) / sqrt(2*pi*sa^2);
    
  hold on 
  plot( x, y, 'm-' )
  
%   max_y = max( opo.y );
%   plot( mw, max_y/100, 'sr' )
%   plot( mw+sa, max_y/100, '+g' )
%   plot( mw-sa, max_y/100, '+g' )

  hold off

end
