% Wende pq-Formel aus SolveEq() auf unterschiedliche 
% q Werte an, studiere Qualitaet der Loesung
% Warum genuegt es q zu variieren?
clc
clear
p=1;
disp('Loese quadratische Gleichung fuer p = 1 und unterschiedliche q')
fprintf('q            xm        xm^2+p*xm+q     xp         xp^2+p*xp + q\n')
fprintf('----------------------------------------------------------\n')
for i=1:16
    q=10^(-i);
[xp,xm] = SolveEq(p,q);
fprintf('%5.0e %20.15e %10.1e %20.15e %10.1e\n',q, xm, xm^2+p*xm+q, xp,...
    xp^2+p*xp+q)
end
