%
%  [x1,x2] = pqFormel( 1, -6 )
%  [x1,x2] = pqFormel( 0, -4 )
% 
function [x1,x2] = pqFormel( p, q )
  x1 = - p/2 + sqrt( (p/2)^2 - q ); 
  x2 = - p/2 - sqrt( (p/2)^2 - q ); 
end