  function  [ discrim ] = discriminant( A, B, C )    

  % DISCRIMINANT : Compute discriminant in quadratic equation.
  %

  discrim = B*B - 4*A*C;
