Função para Cálculo pelo Método Quadrático para Gretl

function [x_neg,x_pos]=quadratica_1(b,c)

  // [x_neg,x_pos]=quadratica_1(b,c)

  // Resolve a equa��o quadr�tica x^2+b*x+c=0 com controle num�rico do erro.

  // (C) Rudnei Dias da Cunha 2009

  if (b<0)

      x_pos = (-b+sqrt(b^2-4*c))/2

      x_neg = c/x_pos

  else

      x_neg = (-b-sqrt(b^2-4*c))/2

      x_pos = c/x_neg

  end

endfunction