p = input('Bitte das p :> ') ;
q = input('und nun das q :> ' ) ;

p_halbe = p/2 ;
radikand = p_halbe^2 - q ;

if radikand > 0
   fprintf('Der Radikand ist positiv\n')
   fprintf('Die Nullstellen von x^2+%6.2fx+%6.2f sind \n',p,q)
   x_1 = -p_halbe + sqrt(radikand) ;
   x_2 = -p_halbe - sqrt(radikand) ;
   fprintf('x_1 = %10.4f \n', x_1) ;
   fprintf('x_2 = %10.4f \n', x_2) ;
   % und jetzt die Probe 
   p_1 = x_1 * x_2 - q
   p_2 = - (x_1+x_2) - p
   fprintf('Probe nach Vieta: %3.2f %3.2f\n',p_1, p_2)
elseif radikand == 0
    % eine Nullstelle
else
   fprintf('Der Radikand ist negativ\n')
   neg_wurzel = sqrt(-radikand) ;
   x_1_real = -p_halbe ;
   x_1_imag = neg_wurzel ;
   x_2_real = -p_halbe  ;
   x_2_imag = -neg_wurzel;
   fprintf('x_1 = %10.4f + %10.4f i \n', x_1_real, x_1_imag)
   fprintf('x_2 = %10.4f + %10.4f i \n', x_2_real, x_2_imag)
end