%
% Serie 03, Aufgabe 2
%
% m.file um die quadratische Gleichung x*x + p*x + q = 0 zu loesen

% Eingabe der beiden Parameter
clear
format long g
p = input('geben Sie den Parameter p ein: ');
q = input('geben Sie den Parameter q ein: ');

% Berechnung und Ausgabe der Loesung

a = -p/2;
D = a*a-q;

% Algorithmus 3.2 aus dem Skript

disp('Algorithmus 3.2')

if (D>=0)
   x1 = a+sqrt(D)
   x2 = a-sqrt(D)   
else
   sprintf('keine reelle Loesung !')
end

% Algorithmus 3.3 aus dem Skript

disp('Algorithmus 3.3')
if (D<0)
   sprintf('keine reelle Loesung !')   
else
   if (a<0)
      x1 = a-sqrt(D)
      x2 = q/x1
   elseif (a>0)
      x1 = a+sqrt(D)
      x2 = q/x1
   else
      x1 = sqrt(-q)
      x2 = -x1
   end   
end
