function answer = discrim( A, B, C )
% discriminant of a quadratic equation

answer = B.^2 - 4 .* A .* C;