%
%  pqFormel2: 
%    p =  0, q = -4
%    p = -2, q = 1
%    p =  0, q = 9
%    p =  'W', q = 9
% 
function pqFormel2()
  fprintf( 'Loesung der quadr. Gl. mit pq-Formel: \n' );
  p = input( '> p = ' );
  if( ~ isnumeric(p) )   % Check, ob Zahl eingegeben wurde
    fprintf( 'Eingabefehler: Zahl erwartet! \n' );
    return;
  end
  q = input( '> q = ' ); 
  if( ~ isnumeric(q) )
    fprintf( 'Eingabefehler: Zahl erwartet! \n' );
    return;
  end
  
  D = (p/2)^2 - q; 
  % 3 Fälle: 
  if( D > 0 ) 
    x1 = - p/2 + sqrt( D ); 
    x2 = - p/2 - sqrt( D ); 
    fprintf( '-> 2 reelle Lsgen.: x1 = %g, x2 = %g \n', x1, x2 );
  elseif( D < 0 )
    Rex = - p/2;
    Imx = sqrt( -D );
    fprintf( '-> 2 komplexe Lsgen.: x1 = %g+%gi, x2 = %g-%gi \n', ... 
             Rex, Imx, Rex, Imx );    
  else  % D == 0
    x = -p/2;
    fprintf( '-> 1 Doppel-Lsg.: x = %g \n', x );
  end
    
end