function [ xp xm ] = SolveEq( p,q )
% Naive Loesung einer quadratischen Gleichung der 
% Form x^2 + p*x + q = 0 mit pq-Formel
%   Input: p,q
%   Output: xp, xm = -p/2 +/- sqrt(p^2/4-q)
xp = -p/2 + sqrt(p^2/4-q);
xm = -p/2 - sqrt(p^2/4-q);
end

